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information related to metastability, 


1. Introduction 


cavitation 
inte 


The prediction of 
ment has been of 
turers, and users for many Essentially all 
of the work in this area consists of testing which 
pro\ ides information that is applicable only to the 
fluid and equipment used in the Investi- 
gations concerned with basic problems of cavitation 
have not been pursued as vigorously. The purpose 
of this paper is twofold: (1) To present an analysis 
which indicates some of the basic problems in cavi- 
tation and which may permit 
cavitation characteristics, and 
experimental results which have 
compare these 

Cavitation « 


in hydraulic 
to designers, 


equip- 
rest manufac- 


years. 


tests. 


To discuss some 
been obtained and 
with the results of the analysis. 
occur in any type of hydraulic 
equipment in the local static can 
drop below the local vapor pressure, with a conse- 
quent formation of a vapor phase. It can occur in 
simple equipment such as a pipe and in the most 
complicate 
turbine. 

Cavitation may be said to occur whenever a cavity 
forms in a liquid; in hydraulic equipment this cavity 
is usually filled with vapor from the surrounding 
liquid, or with a mixture of this vapor and gaseous 
impurities. Cavitation is frequently associated with 
certuin symptoms. To some, cavitation exists as 
soon one, very few, small vapor bubbles 
appear in the system; to others, cavitation exists if 
damage results; to a third group, cavitation exists 
when audible sounds created. Clearly, the 
symptoms of interest must be specified if the concept 
of the “existence of cavitation” is to be significant. 

This study will be concerned with cavitation that 
is sufficiently developed to cause a detectable alter- 
ation in the performance of equipment. Usually, 
an appreciable amount of vapor will form and there 
may be damage or noise. The symptoms of cavita- 
tion to be predicted are therefore changes in per- 
formance characteristics such as the head-capacity 
characteristics of a pump, and the calibration curve 
of a flow meter. 

Examples of the empirical work 
tation include that of 
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problems in cavitation 


and which may 
of hydraulic 


equipment is presented. 
results of the analysis. 
prediction of symptoms of 
), but that 
and vapor-phase dynamics is required. 


Weltmer [3], Hartman [4], and Wood [5]. In these 
investigations the application of the results is limited 
because an insight into the basic phenomena is not 
a primary objective of the experiments. <A limited 
number of studies Stahl and Stepanoff [6]) 
have attempted to gain more widely applicable infor- 
mation through more detailed considerations of the 
phenomena in systems where appreciable cavitation 
occurs. Here the influence of fluid properties on 
cavitation characteristics will be emphasized. (Fluids 
of significantly different properties are considered 
as different fluids, even though they 
chemical composition. ) 


to the 


(oe 
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have the same 


2. Analysis 


Consider the device in figure 
ric equipment such as a 
Regardless of the liquid, 
be the same at the 


1, a piece of volumet- 
pump, venturi, or valve. 
the volume flow rates will 
same non-cavitating operating 
point. Viscous and compressibility effects are 
assumed negligible (experiments by Richards [7], 
Purcell [8], and Martin [9] show that this assumption 
is acceptable Under the same assumptions it fol- 
lows (and experience confirms) that the static pres- 
_ gradients within the device (expressed in height 

f liquid) are essentially independent of the liquid. 
F or example, in a non-cavitating pump operating at 
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FIGURE 1 Pertinent lecations and head drops in a cavitaling 


device. 
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a given speed and capacity, the static head-drop 
from the inlet to any point within the pump is inde- 
pendent of the liquid. 

Consider the pressure level in the device to be 
lowered until a change in performance, due to 
cavitation, is incipient. This does not mean how- 
ever that the vapor phase may not be well developed. 
Immediately before the symptoms of cavitation 
become evident the pressure and velocity distribu- 
tions in the device (to location 2 in fig. 1 where the 
vapor phase is starting to form) do not depend upon 
the fluid. Thus, the differences in cavitation 
characteristics of different fluids, in the same 
device, are due to the differences in what occurs in 
those regions of the device that are downstream of 
the point of vapor inception. (Note that differences 
in nucleation characteristics of liquids can alter the 
location of the point of vapor inception for different 
fluids.) Attention therefore focused on the 
differences in behavior of the fluids after the static 
pressure has dropped to, or perhaps somewhat below, 
the vapor pressure. The concern here is with what 
occurs between locations 2 and 3 in the device in 
figure 1; if stable equilibrium exists, locations 1 and 
2 coincide. 

The work of Stahl and Stepanoff [6] has the same 
purpose as this analysis, but the basic assumptions 
and therefore the final results of the two analyses 
differ. The major difference is that: 
the head depression, below the pressure at which 
vaporization starts (not necessarily the 
pressure), required to generate a ‘“‘cavitating volume’”’ 
of vapor is calculated, and it is lesa that differ- 
ences in this head account for the variations in 
cavitation characteristics of different fluids. Stahl 
and Stepanoff calculate the vapor-to-liquid-volume 
ratio formed because of a given head depression, and 
assume that this ratio is an indication of cavitation 
characteristics. In addition, this analysis accounts 
for characteristics of the hydraulic device. Sale- 
mann [10] has discussed some aspects of the anal- 
ysis presented here, as well as that of Stahl and 
Stepanoff. 


is 


2.1. Nomenclature 


C'=defined as t ‘ott P|: tem- 


perature coefficient for determining 

changes in the enthalpy of the liquid 

during the vaporization process. 
(’,=isobaric specific heat of liquid. 


v7; (l— T 181) ( dp 


h=specific enthalpy of fluid in region of 
cavitation. 

N=pump speed. 

p=pressure. 

(/=pump capacity. 


7T=temperature. 


u=mean fluid velocity in the region of cavita- 
tion. 
V=volume of vapor per unit volume of mixture 


in region of cavitation when symptoms 
of cavitation are incipient. 
v=specific volume. 





In this analysis | 


vapor | 





xz=quality of fluid in region of cavitation. 


z-=quality in the region of cavitation when 
symptoms of cavitation are incipient. 

Ah,.=drop in head, measured in height of liquid, 
which is necessary to cause incipient 
symptoms of cavitation. 

Ah,=head drop experienced by a particle of 
fluid as it travels from the pump inlet to 
the region where cavitation is occurring. 

Ap-=pressure drop, after vaporization begins, 
necessary to cause sufficient flashing to 
produce incipient symptoms of cavita- 
tion. 

sg ; ;, 1 / ov 
8=coefficient of expansion, (= :) 
v\or/, 

5=a function of ¢ and fluid properties, defined 

by eq (13). 
A\=latent heat of vaporization. 
p=density. 
@=defined as V/(1—V), and is generally a func- 

tion of design, operating point, and fluid. 
x=a function defined by eqs (7) and (8). 


Subscripts: f refers to saturated liquid. 

g refers to saturated vapor. 

/ refers to liquid. 

v refers to vapor. 

a eB 
fluids. 


nN used to identify 


2.2 


Assuming that steady flow exists, heat transfer and 
shaft work are negligible during the flashing process, 


Derivation of Equations 


and that effects of force fields are negligible, the 
energy equation may be written 
uu? 
dh+d(5)=0. 
Assume further that h=h,+-r\, \ is constant, and 


that temperature is a function of pressure only. 
The energy equation becomes: 


a(S L ) dp +-Ada+d ( =) =(), 

The temperature coefficient, C, is defined as 
dl} ; athe 2 ; 

bi ( [ Pp) w1—T18) (J) | 


dT 
The derivative (|) is the slope of the 


ip (1) 


‘dp 


vapor- 


pressure curve only if the liquid and vapor are in 
stable equilibrium; generally, this is not the case, and 
the function 7(p) will depend upon those properties 
of the fluid which determine its ability to exist in 
metastable states. Also \ may not be exactly equal 
to the usual (stable equilibrium) latent heat of 
vaporization. 

Assuming that steady state exists, and that viscous 
drag, changes in cross-sectional area, and changes 
in flow direction can be neglected during the flashing 
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process, the momentum equation can be written 


dp t d(pu?) Q. (2) 
The limitations introduced by the preceding assump- 
tions may be relaxed in the application of the results 
of this analysis because the characteristics of differ- 
ent fluids passing through the same device are to be 
compared; it is anticipated that the fluid properties 
will be important and that the geometrical effects 
will tend to cancel out. 

Using the previous assumptions, 
equation may be written 


the continuity 


d(pu)=0. (3) 


and = that 
, and combining equations 1, 2, and 


Assuming that = 1/p=[(1—.)v,+ z2,| 
z lore v 
3, we have 


] dr 
ap — 
(dl 
I Jl ( ( dp | 
Note that 7, and v, are the specific volumes of 


saturated liquid and vapor only if the fluid is in stable 


é qu ilibr ium ; 


dT 
Assume that 7, v,, C, and ( ] ) are not 
( pP 


functions of x. The preceding expression can then 
be integrated. Evaluating the constant of integra- 
tion from the condition that Ap=0O when «s=0, we 
obtain 


dT 
(i) 
It is important to note that the pressure at which 
Ap=0 is not the local vapor pressure, but depends 
upon the metastability and nucleation characteristics 
of the liquid. As characteristics vary from 
liquid liquid, it follows that the initial pressure, 
and hence the region in the device at which eq (4) 
to be applied (i.e., between locations 2 and 3) depends 
upon the fluid. Equation (4) gives the change in 
pressure (Ap) which must occur if a fraction x of 
pure liquid is to vaporize by flashing. 

Recall the premise that the symptoms of 
cavitation will occur when a ‘‘cavitating volume”’ of 
vapor forms in some region of the device. The loca- 
tion of the region and the volume of vapor depend 
upon the design of the device, the operating point, 
and probably upon the fluid; the location can be in 
the eye of an impeller, in a small localized low pres- 
sure region, or across a flow passage. If 2, is the 
quality in the cavitating region, then z,”, is the vol- 
ume of vapor per unit mass of mixture in the region. 
Therefore 


these 


is 


basic 


~ 








where Ap, is the pressure drop, after vaporization 
begins, necessary to cause sufficient flashing to pro- 
duce symptoms of cavitation. 

Note that 7,7, is the volume of vapor formed per 
unit mass of mixture. As equal volumes, not equal 
masses, of fluid flow through the device per unit 
time, we should be concerned with cavitating volumes 
of vapor per unit volume of mixture. Because of 
previous assumptions, we may write (z.v,)~v[V/ 
(1—V)], where V is the volume of vapor per unit 
volume of mixture in the cavitating region. (This 
V is the same as V’ mentioned by Salemann [10].) 

In eq (5), Ape is expressed in pressure units; Le., 
force/area. "This quantity could be directly com- 
pared from fluid to fluid if it were expressed in height 
of the flowing liquid. Substituting for (2.v,) and 
Ap, in eq (5), we obtain 


(x 


Vy 
(i—V) 


“ i : dT | 


dp 
Equation (6) gives the drop in head Ah,, 
height of liquid, which would accompany the forma- 
tion of a volume, V, vapor per unit volume of 
mixture in the cavitating region, and which is there- 
fore necessary to initiate symptoms of cavitation. 

One is tempted to Say that, for a given device and 
operating point for a given speed, N, and 
capacity, QY, in a pump) the vapor volume ratio, V, 
may be the same for all fluids. In view of what has 
been said previously concerning metastability, nucle- 
ation, and vapor-phase dynamics, V is probably also 
a function of the fluid. Defining a function 


(6) 


measured in 


of 


(@.2., 


@ (operating point, <9 


fluid) = ) 


we have 


‘) 
Av _ 
Ah, -In 1+- ? — ? (7) 
ee ) 
v; dp 
where the function ¢@ is to be determined either 


experimentally or by further analysis. The analvsis 
presented by Wu [11] is an example of a theoretical 
approach which can predict the vapor volume ratio, 
V. Wu determines the effect of an attached region 
of vapor upon the lift and drag coefficients of a body 
in the flow passage. When the volume of this region 
is great enough to appreciably effect the lift and drag 
of the body, the performance of the device (e.g., 
pump) will be affected; this computed volume will 
then be the “vapor volume” in the cavitating region. 


If no information is available regarding metasta- 
bility, nucleation, and vapor- -phase dynamics, the 
functions T(p) and C cannot be predicted. In 


addition, if the fluids to be compared do not cavitate 
similarly (i.e., @ depends upon the fluid) eq (7) may 
be written in the form 


149 








\» 
Ah,= “7 x (8) 


0 


(operating point, fluid), 


where the function x must be obtained by cavitation 
tests on the device at different operating points and 
with different fluids. 


2.3. Application to Pumps 


the results obtained in 
can be used, and to indicate 


In order to indicate how 
the preceding section 


some of the unresolved problems, the application of | 


the results to the prediction of symptoms of cavita- 
tion in pumps will be discussed. 

First it-is necessary to relate the head drop, Ah,, 
with some measurable quantity which is related to 

cavitation characteristics. The relations formulated 

vas are primarily for illustrative purposes; their 
predictions compare poorly with the available experi- 
mental data. Possible reasons for this will be indi- 
cated as the formulation develops. Because of its 
wide use in pumping, the net positive suction head 
(NPSH) (i.e., the difference between the total 
pressure and the vapor pressure at the pump suction) 
will be used as the measurable quantity. During 
the aonuiader of this discussion, the term (NPSH) 
will mean the minimum NPSH required to suppress 
symptoms of cavitation. 

One way to set up the required relationship is as 
follows: Consider the drop in head, Ah,, experienced 
by a particle of liquid as it travels from the region 
where the NPSH is measured to the region where the 
cavitation is occurring. (See fig. 1.) For a given 
pump and given operating point, Ah, is assumed to 
be about the same for all fluids if the pump is not 
cavitating. (It has been implied that the location 
of the cavitation does not depend upon the fluid, 
which is not generally true.) Assume that 


Ah N PSH)—AA,. (9) 
(Both Ah, and (NPSH) are positive, while Ah, is 
negative. ) Noting that NPSH represents the head 


drop from the inlet to the vapor pressure, and that 


Ah. is reckoned, not from the vapor pressure, but 
from the pressure at which vaporization begins 
(refer to fig. 1), it is apparent that eq (9) is not exact. 


In addition, the alteration of the flow pattern within 
the pump due to the presence of the cavitating fluid 
is ignored. If, in spite of these limitations, Ah, is 
assumed independent of the fluid being pumped, it 
follows that 
(NPSH), 


(Ah): = (NPSH).— (Ah) 2 


(NPSH),,—(Ah,),=const., (10) 


where the subscripts 1, 2,...m... designate 
different fluids. 

If the required fluid properties (A, %, v,, T'(p), ete. 
are known, the prediction of the N PSH requireme ot 
for a given pump with any fluid by means of the 
results of this analysis requires the determination of 





two functions, Ah, and @ (or x). If both functions 
can be predicted theoretically no experimental work 
is required; this will probably not be the case. 

A second relationship between NPSH and Ah, 
is presented here because it approximately predicts 
some NPSH’s which have been measured with liquid 
hydrogen, liquid nitrogen and liquid oxygen; its use 
has not yet been justified on any other basis. Be- 
cause (N PSH) Ah.) increases, 
we assume that 


should decrease as ( 


(NPSH)(Ah.) =a constant. (11) 
Thus, 
a RES Ah, 
{ > 9 
NPSH),—(NPSH) 4 ‘} (12) 


The comments made above in connection with eqs 


(9 and 10), concerning the determination of the 
unknown functions, apply here also. Equations 
7 or 8) and either (9, 10, 11, or 12) constitute pairs 
of relations which are examples of results of this 


analysis. 
Due to the form of eq (7 
function @ from eq (10 or 12 


, the computation of the 
) involves the solution of 


a transcendental equation. In order to simplify 
computations an approximation to eq (7), based 
upon the assumption that 
C (dT 
o<fi-S(4) 
v) dp 
will be derived. The error introduced by the 
assumption is not more than a few percent, the 


magnitude of the error depending upon the magni- 
tude of the function ¢. Define: 


(a 
l 
EE i (ip) | 


(13) 


)] 


, we can write 
IY) l l l 
Ah.=— if =6+58—... ;3 | (14) 
7 2 » i 


The limitations of this restriction (i.e., 61) must be 
determined through numerical computation. As 
there is no information available describing a fluid 
which deviates from stable equilibrium, stable 
equilibrium properties must be used in the computa- 


a \. ; 
tions; ( j -) is then the slope of the vapor pressure 
ap 


curve, and the other properties are for saturated 
liquid and vapor. ‘The data required for the compu- 
tation and the results for three important cryogenic 
fluids at atmospheric pressure are given in table 1. 
The computations show that the upper limit cand 
upon the volume fraction of vapor in the cavitating 
region, V, to ensure the validity of eq (14) is very 
high. It is higher than any values for V computed 
from the cryogenic data obtained by the author, and 
much higher than any values computed by Salemann 
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[10]. ‘This limit permits more than 99 percent of the 
cavitation region, with nitrogen or oxygen, to be 
filled with vapor, and more than 95 percent of the 
region to be filled with vapor in the case of hydrogen. 
Therefore this limitation on V (or ¢) is probably 
unimportant. 

Replacing eq (7) by eq (14) does not simplify com- 
putations because the transcendental function is 
merely replaced by its infinite series representation. 
A further development, which answers the following 
question, must be made: How many terms in eq (14) 
can be neglected without introducing unacceptable 
errors? The simplest computation results if only the 
first term is used. The maximum error introduced 
in the computed value of Ah, by neglecting all terms 
except the first is (1/2)6; this error is plotted as a 
function of V in figure 2. Curves are shown for 
hydrogen, nitrogen, and oxygen at their normal boil- 
ing points, and for water at 70 °F. The error with 
water is negligible, even when 95 percent of the 
cavitation region is filled with vapor. With nitrogen 
and oxygen the error is only 6 percent when 95 per- 
cent of the cavitation region is filled with vapor, and 
is less than 1 percent when there is less than 75 per- 
cent vapor. ‘The largest error, with hydrogen, is less 
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Figure 2: Percent error in Ah, due to use of only the first 
term in eq (14) versus the volume fraction of vapor (V) in 


the region of cavitation. 


than 10 percent when the volume of vapor is less than 
82 percent. The values of V computed from our 
| data and those computed by Salemann [10] indicate 
that the use of only the first term in eq (14) will 


yield reasonably accurate values of Ah,. There- 
fore Ah, will be computed from 

hv, dv o 

Ah, 6 = aed 15) 

’ ' C/di (15) 


bed bes) 


_ Substitution of eq (15) in eqs (10 and 12) gives the 
simplified formulations for the prediction of NPSH: 


(NPSH), sane 
| 2, [ 
l 


v1 


dp 








| 
| 
| 
NY, QD 
— - oe ( b) 
v C' /d7I 
iF 
v dp 1 
and 
NY, .9) 
v, & dT 
fara by 
‘ ; . . 0; \ap 1 - 
(NPSH).=(NPSH), - (17) 
NY, re 
v c; dT 
ae 
v, \dp . 
TasLe 1. Cryogenic fluid properties and data 
Quantity H ydro- Nitrogen Oxygen 
gen 
Boiling point, °K 20.4 77.32 90. 13 
Latent heat of vaporization, A, cal/g 106. 5 47.6 50. 8 
Isobaric specific heat, Cp, cal/g °K -. 2 50 0. 489 0. 400 
Specific volume of liquid, », cm4/g 14,11 1, 237 87 
Thermal expansion coefficient, B;, 1/°K 0.0175 0. 00580 . 00738 
Ratio of specific volume of saturated vapor to 
that of saturated liquid, v,/ry__- 5O 183 265 
iT 
Slope of vapor pressure curve ( dp ), P 
K cm‘/cal 132 383 336 
Upper limit on V for eq (14) to be valid 0. 958 0. 994 0. 994 
Minimum NPSH_ required to suppress 
symptoms of cavitation, ft 0. 167 7.0 10. 6 to 


Equation (16 or 17) would permit the prediction 
of the NPSH required to prevent symptoms of 
cavitation with fluid 2 if the NPSH requirements 
for fluid 1, the required fluid properties, the 
¢-functions were known, and if the assumed relations 
between NPSH and Ah, (eqs 9 or 11) were valid. 
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If it is further assumed that ¢ is not a function 
of the fluid properties, eqs (16 and 17) simplify to 


m 
it 
(NPSH),—(NPSH),+4(N,Q) 4 | — ( ap 
~ oXdp/ |, 
2! 
v. : 
ee as i (18) 
C/dT 
tnd ip) 
and 
+ : 
fe) 
(NPSH),=(NPSH), =>—"““P/s" = (NPSH), (22), 
rn T) 
0, 
dT 
a -( 7) 








(19) 


“ec 


respectively, where 7 has been called the “cavitation 
tendency”. The greater the cavitation tendency, 
the greater the required NPSH. The assumptions 
involved in the derivation of eqs (18 and 19) must 
be emphasized: the location of the cavitation region 


in a pump is independent of the fluid; the volume of 


vapor in this region, when cavitation symptoms 
arise, is independent of the fluid. The NPSH 


requirements which are to be compared are at the 
same operating points (i.e., the speed (N) and 
flow (Q) are the same for both fluids). In applying 
eq (18), the function, ¢(NV,Q) must be known, while 
application of eq (19) does not require a knowledge 
of the function. 


3. Experimental Comparison 
3.1. Centrifugal Pumping 


In this section some experimental results which 
are pertinent to the analysis presented in the preced- 


ing section will be discussed. These indicate that 
the effects of the existence of metastable states, 


nucleation characteristics, and vapor-phase dynamics, 
are significant. 

Only a meager amount of data is available to 
check the validity of (18 and 19); some data 
have been reported by Salemann [10], and some have 
been obtained by the author. As no information is 
available concerning the deviation of the fluids from 
states of stable equilibrium during the tests, prop- 


eqs 


erties of saturated fluids are used in making 
comparisons. 
One of the more drastic assumptions in the 


derivation of eq (18) is that @ is independent of the 
fluid being pumped. A straightforward way to 
check the equation is to substitute measured values 
of the minimum required NPSH and known values 
of fluid properties into the equation, and to determine 








if o(N,Q)[=V/(1—V)] is constant. Salemann [10] 
did this for trichloromonofluoromethane, water, 
butane, and benzene, and found that ¢ varied from 
0.25 to 1. 


In order to compare eq (19) with the data, we 
shall rewrite it as 
» (Led) (D=aGy)I 
(NPSH), _T »\dp / | (20) 
(NPSH), 1 . 


AEN] 


The validity of this expression can be checked by 
comparing the left-hand side (the experimental 
NPSH ratios) with the right-hand side (the theoreti- 
cal cavitation-tendency ratios computed from fluid 
properties). The results of this comparison, using 
some of Salemann’s data, are given in table 2. It 
apparent that Salemann’s data indicate that 
eqs (18 and 20) are not quantitatively valid, although 
the second equation may predict relative NPSH 
requirements for different fluids. 

During the development of a hydrogen pumping 


is 


system at this laboratory some crude NPSH data 
for hydrogen, nitrogen, and oxygen were obtained 


in a centrifugal pump designed for water. These 
data will be compared with eqs (18 and 20). Figure 3 
is a schematic of the test apparatus used to obtain 
the data. The pump was a 10-stage submersible 


water pump rated at 3,450 rpm, 9.7 gpm, and 100 ft 
NPSH 


above 


discharge head. The 
adding the liquid level 


was 
the 


determined by 
pump suction, 
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measured with a hydrostatic device, to the pressure 
in the vapor space above the liquid. Readings were 
taken immediately after pressurization to minimize 
the rise in temperature of the fluid. The pressure 
in the vapor space was reduced to atmospheric 
between each set of readings so that the temperature 
would drop to the saturation temperature. (In a 
later apparatus the NPSH was determined by a 
vapor-pressure device [12] and the state of the fluid 
at the pump suction was controlled by a _ heat 
exchanger and a standpipe.) As the liquid in the 
test Dewar was very agitated, its temperature was 
assumed to be uniform at the saturation temperature 
corresponding to atmospheric pressure. 

TABLE 2. Comparison of NPSH ratios to cavitation-tendency 
ratios for some of Salemann’s [10] data 


Experi- (Theoretical 
Fluids mental | cavitation 
(NPSH) tendency 
ratio ratios 
Freon-11 at 85 °F A 1.21 2 6 
Freon-11 at 120 °1 
Water at 250 °1 1. 28 4.10 
Water at 300 °F 
Water at 250 °I paeda 1.31 1. 80 
Freon-11 at 120 ° 





The minimum NP SH’ s require dt to suppress symp- 
toms of cavitation were obtained by plotting suction- 
head depression curves: for a given operating point 
the capacity je discharge head were plotted as 
functions of NPSH. Above the minimum required 
NPSH both the capacity and discharge head are 
independent of NPSH, while below this value both 
quantities decrease very rapidly (with decreasing 
NPSH). The breaks in both curves occurred at the 
same NPSH and were readily located on the graphs. 
The symptom of cavitation was taken to be the break 
in the curves; the minimum NPSH required to sup- 
press symptoms of cavitation defined the 
NPSH adjacent to, and above, thie breaks in the suc- 
tion-head depression curves. (Figure 4 is included 
to indicate the type of data obtained.) 

The measured values of minimum NPSH required 
to suppress symptoms of cavitation are given in 
table The low value for hydrogen has since been 
verified on other pumps by the author and by other 
investigators [13]. The minimum NPSH required 
for oxygen a only be bracketed between 10.6 ft 
and 15.7 ft because the tests were terminated when 
it was discovered that the liquid oxygen was eroding 
the pump. 

Using data in table 1, we can compute values of 
@ from eq (18), and then values of V from the 
definition of ¢ Equation (18) satisfies our hydrogen 
and nitrogen data if the volume fraction of vapor, 
V, is 0.053. It satisfies our nitrogen and oxygen 
data if the volume fraction of vapor is between 0.848 
and 0.930. Thus, these data lead to the same con- 
clusion as the data of Salemann: that the formulation 
expressed by eq (18) (in which @ is assumed to be 
independent of the fluid being pumped) will not 
predict NPSH requirements. 
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Rated flow. 


The comparison between the cryogenic data and 
eq (20) is given in table 3. The theoretical predic- 
tions are in qualitative agreement with the experi- 
mental results. It is not possible to say whether 
more refined experiments would produce better or 
worse agreement. 

Thus, both Salemann’s data and the cryogenic data 
indicate that eq (18) is invalid, while both sets of 
data indicate that eq (20) may have some value in 
predicting NPSH requirements. It should be 
pointed out that Salemann’s experimental (N PSH) 
requirements were those which permitted a 3 percent 
drop in discharge head, not those required to prevent 
symptoms of cavitation entirely. The less satisfac- 
tory correlation of Salemann’s data with eq (20) may 
be due, in part, to the possibility that a 3 percent 
decrease in performance appreciably affects the flow 
throughout the whole pump, and thus seriously vio- 
lates assumptions made in the derivation of the 
equation. It should be noted that the assumptions 
that Ah, and @ are independent of fluid character- 
istics are very drastic, and may be responsible for 
the poor correlations presented above. 


3.2. Flow Measurement 


It has long been known (e.g., see Benjamin and 
Miller [14]) that even though a liquid passes through 
a region in which the static pressure falls below the 
vapor pressure, symptoms of cavitation need not 
occur. This means that even if the pressure at the 
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vena contracta of an orifice or at the throat of a 
venturi is below the vapor pressure, such devices 
may still be accurate flow meters. As this situation 
exists in many applications it is desirable to be able 
to predict when cavitation will affect the performance 
of a flow meter. 


TaBLeE 3. Comparison of NPSH ratios to cavitation tendency 
ratios for cryogenic liquids 
Experi- Theoretical 
Fluids mental cavitation 
NPSH tendency 
ratio ratios 
Nitrogen: hydrogen 42 5 
Oxygen: nitrogen 1.51 to 2.24 1. 37 
Oxygen: hydrogen__- 64 to 94 75.5 


During an investigation to determine the be- 
haviors of sharp-edged orifices with water, liquid 
nitrogen, and liquid hydrogen [7], it was not possible 
to produce cavitation symptoms as long as pure liquid 
entered the orifices. With liquid nitrogen, the pres- 
sures at the venae contractae were as much as 170 
in. of liquid (4.9 psi) below the vapor pressure, while 
with liquid hydrogen the pressures at the venae con- 
tractae were as much as 192 in. of liquid (0.49 psi) 
below the vapor pressure. (These were the lowest 
pressures attainable with the apparatus.) The only 
way that symptoms of cavitation could be produced 
was to have two-phase flow entering the orifices; in 
many tests, even when two-phase flow entered the 
orifices, cavitation symptoms were not evident. 
Figure 5 is included here to show the results of some 
of the tests. Notice that only about one half of the 
nitrogen test points for which the upstream static 
pressure was below the vapor pressure show evidence 
of the presence of two-phase flow. Neither of the 
hydrogen points with upstream pressures below the 
vapor pressure appears to be affected by the two- 
phase flow. It is also interesting that orifice cali- 
bration curves obtained with water may be quite 
satisfactory for use with liquid nitrogen and liquid 
hydrogen. 


The absence of cavitation symptoms when the 
static pressure is below the vapor pressure is due to 
the fact that appreciable time is required to form a 
significant vapor phase and that the fluid elements 
are not in the low pressure region for the required 
length of time. The time depends upon the ability 
of the fluid to persist in metastable states, on nu- 
cleation characteristics, and on bubble dynamics 
including heat and mass transfer. Thus the same 
basic problems are involved here as with cavitation 
in pumps, and the analysis presented in section 2 
may be applicable. 


3.3. Two-Phase Flow 


Except in the region where the vapor phase is 
initiated, two-phase single-component fluid flow is 
not normally associated with cavitation. However, 
as metastability is one of the basic phenomena asso- 
ciated with cavitation and as it has been observed 
that metastable states can exist throughout two- 
phase, single-component fluid flows, some observa- 
tions of these flows are presented here. The purpose 
in this seetion is to show that fluids can exist in 
metastable states for relatively long times; therefore, 
information of the types discussed in preceding sec- 
tions would also be useful in this area. 

The observations reported here were made in an 
apparatus designed study steady, two-phase, 
single-component, fluid flow. This apparatus, de- 
scribed by Hatch [15], employs trichloromonofluoro- 
methane. The curves in figure 6 show the existence 
of metastable states in the test section. The ‘meas- 
ured temperature’? curves were obtained from actual 
temperature measurements, while the ‘saturation 
temperature’ curves were obtained from the satura- 
tion temperatures taken from tables of thermody- 
namic data) corresponding to the measured pressures. 
It is obvious that the measured temperature lagged 
behind the saturation temperature, resulting in a 
superheated fluid. Here again, it can be concluded 
that vaporization does not occur rapidly enough to 
permit the attainment of stable equilibrium, 
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Calibration curve. 
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Measurements indicating differences between 
measured and saturation temperatures of more than 
3 °F were common. If the two-phase pressure drops 
are computed from the Martinelli and Nelson [16] 
correlation, the existence of these metastable states 
can cause errors up to 15 percent. . These observa- 
tions stress once again the desirability of research 
into metastable equilibria. (It should be pointed 
out that the spatial oscillations of the measured 
temperatures in figure 6 are believed to actually occur; 
no fault could be found with the thermometry. 
They may be due to the nozzle at the inlet to the test 
section.) 


4. Conclusions 


Observations of hydraulic equipment in which 
‘avitation occurs show that the existence of metasta- 
ble equilibrium states, nucleation characteristics, 
and vapor-phase dynamics can be important in 
the prediction of cavitation characteristics. In- 
vestigations of fluid characteristics which influence 
such behavior are therefore necessary if cavitation 
symptoms are to be predicted. 

The meager data available indicate that an 
approach of the type presented in section 2 may be 
fruitful in predicting svmptoms of cavitation. 
In addition to the types of investigations suggested 
in the preceding paragraph, theoretical and experi- 
mental studies concerned with the determination 
of the ¢-function (or volume fraction of vapor in 
the cavitating region, V) should be continued. 
These studies should include work concerned with 
detailed pressure distribution within the equip- 
ment and the influence of the vapor region upon 
the detailed fluid mechanics of the equipment 
involved. 
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Tu o-phase flow data, 


The qualitatively accurate predictions of the 
“Cavitation Tendency” concept suggest a desira- 
bility for further evaluation and refinement of 
this concept. However, the assumption that the 
¢-function is not dependent upon fluid properties 
appears to be too drastic and eqs (16) and (17) should 
be used instead of eqs (18) and (19) whenever the 
¢-function is known. As eq (16) was _ logically 
developed, while eq (17) was not, the former formula- 
tion is to be preferred; its poor correlation with 
available data may be due to the drastic assumption 
(i.e., @ independent of fluid) mentioned above. 
Finally it is suggested that, when sufficient informa- 
tion is available, eq (9) (and eqs (10) and (16)) should 
be modified to account for the delay in vapor forma- 
tion which occurs between locations 1 
figure 1. 


and 2 in 
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Fresh or outside air 
through shafts or tunnels. For a shaft or 
from the air to the rock in summer 


heating svstem of the installation cannot be 
means of estimating the influence 
computing the 
assuming a sinusoidal variation of the 
substantial accordance 
model tunnel. 


1. Introduction 


When outside air for ventilating an underground 
chamber, or for other purposes, is drawn continuously 
through a long tunnel or airway in rock, heat inter- 
change between the rock and the flowing air con- 
siderably reduces, by the time the air reaches the 
far end, the fluctuations of temperature to which 
the entering outside air was subject as a result of 
diurnal and seasonal changes. The tempering thus 
effected by the rock in cooling summer air and heating 
winter air may have an important benefit in reducing 
the design size and operating cost of the air cooling 
and heating systems of an underground installation. 
Mine and tunnel air ventilation studies [1, 2]! have 
shown qualitatively the reduction of the fluctuations 
of temperature of outside air flowing into mines, 
but have ventured no analysis of this process for 
prediction of temperatures and heat flow quantities. 

This paper presents the results of a mathematical 
analysis of the heat transfer between the air and the 
rock, assuming that the temperature of the entering 
air varies periodically and sinusoidally. For the 
case of long tunnels in rock, diurnal and short-time 
climatic fluctuations are quickly attenuated and 
damped out in the first part of the tunnel. Con- 
sequently, the significant fluctuation is the seasonal 
or annual periodic variation in outdoor air tempera- 
ture, which is known to approximate a sinusoidal 
oscillation [3]. 

The mathematical solutions obtained are for the 
ease of dry air. Efforts to take account of heat 
exchanges between the air and rock due to con- 
densation of moisture from the air, or its evaporation 
from the rock, introduce mathematical complexities 
that preclude analytic solutions. However, the 
results obtained for dry air are conservative in 
regard to the amounts of heat transferred between 
the air and rock, as the effect of latent heat transfer 

*This work was supported by the Protective Structures Branch, Office of the 


Chief of Engineers, U.S. Army 
Figures in brackets indicate the literature references at the end of this paper. 


realized in the design stage without 
of the tunnel. 
tempering effect are developed by mathematical analysis of the 
outside 
with experimental results obtained in tests made in 
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needed for ventilating underground installations is introduced 
tunnel in continuous use, heat is transferred 
and from 
reducing the seasonal cooling and heating loads of the underground installation. 
benefit of this tempering effect in reducing the size 


the rock to the air in the winter, thereby 
The possible 
air cooling and 
a reliable 
In this paper, equations and functions for 
problem, 
shown to be in 
a small-scale 


and operating cost of the 
air temperature, and are 


is to increase the total heat exchanges to values 
greater than those given by the equations for dry air 
conditions. 

To determine the validity of the analytical 
treatment for design purposes, experiments were 
conducted on a small-scale model underground 
tunnel, and the data thus obtained were compared 


with the results predicted by the analytical method. 


2. Analysis 


To obtain a mathematical expression for the heat- 
transfer characteristics of the rock surrounding an 
air-intake tunnel, let us assume that the tunnel 
takes the form of a hollow cylinder. Let us further 
assume that the temperature 7, of the air entering 
the tunnel varies periodically about its mean annual 
value 7; in accordance with the following relation- 
ship: 

(T= T) =(Toeisc T;) cos wt 
where the angle wt increases by 27 radians (360°) in 
l year. 

For convenience in the development, we will use a 

subrogated temperature 6, defined as (7—T,), 


where 7, is the steady temperature of the rock 
remotely surrounding the tunnel. On this basis, 


the subrogated annual mean temperature 6; of the 
entering air is (7';—T7,); the amplitude of entering 
air temperature variation (To max — T;) is designated 
as A. Accordingly, the equation above becomes: 


T)— T,=%=6:+A Cos at. (1) 


The hollow cylinder has a radius, a, but is other- 
wise unbounded. The surrounding solid has a 
conductivity K, specific heat c, and thermal dif- 
fusivity a. The air flows in the cylinder at a velocity 
V in the direction of the positive z-axis. It is 
assumed that in the surrounding rock, heat is 
conducted radially but not parallel to the axis of 
the cylinder, and the thermal properties are in- 
dependent of temperature. 
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Let M’ be the mass of air in contact with unit 
surface area of the cylinder, c’ the specific heat of 
the air, 6, the air temperature at a distance r and 
time t, and A the coefficient of heat transfer between 
the air and the cylinder surface. The heat balance 
equation satisfied by the temperature of the air at 
distance z is then: 


ro dé, 16, 
oui E fe 


where 6, is the temperature at the surface of the 


+h(0,—0, ,)=0 (2) 


cylinder at distance x. The heat transfer at a 
distance x along the cylinder surface is: 
of 
—K( ne ) h(0,—6,, 7). (3) 
dr ° 


Equations (1), (2), and (3) provide the boundary 
conditions for the differential equation for conduction 
of heat in the solid: 


0,2, 108,., 126, 


- 120), 
or Tv? or a Ot 


A closely approximate solution for the temperature of 


the air at distance x, assuming steady periodic heat 
flow, is: 
T.—T,=0,=¢ 


[A cos (wt—AB’)+6,; cos AB’) (4) 
where A involves z, and A’ and B’ involve h. 

The derivation of eq (4) is given in section 8.1. 
The symbols are defined in section 8.3. The functions 
A’ and B’ are plotted on figure 1 for a useful range of 


the variables 2(=ayw/a) and n(=ah/K). 


3. Model Underground Tunnel 

The model underground tunnel (see fig. 2) was 
cylindrical, 6 in. in diameter and 33.7 ft long, sur- 
rounded by 18 in. of 1:4 cement-sand mortar. 
Several thicknesses of !4-in. hair-felt insulation were 
placed on the outside of the mortar. Air was drawn 
through the tunnel by a blower and its rate of flow 
was measured by means of an orifice in the sheet 
metal duct connecting the tunnel and the blower. 
At the entrance to the tunnel was a sheet metal 
duct containing chilled water coils and electric strip 
heaters. The chilled water coils were used to cool 
the air to a temperature below the lowest tempera- 
ture to be encountered at the entrance to the tunnel. 
The electric strip heaters were governed by a pro- 
gram controller that modulated the heat input to 
the air so that the entering air temperature varied 
sinusoidally as a function of time. Thermocouples 
were located at the entrance and exit and at various 
locations in the tunnel, and in the surrounding 
cement mortar, as shown in figure 2. The thermal 
properties of the cement-sand mortar were deter- 
mined from samples obtained during pouring 
operations. 
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Figure 1. Values of A’ and B’ as functions of 2(=avyw/a) 


4). 


and hi ah K) for use in eq | 

Each test consisted of four cycles of temperature 
variation with a period of 16 hr and an amplitude 
of about 20 °F. Because the transients due to the 
initial nonperiodic condition were virtually elimi- 
nated in the first two cycles, the last two cycles were 
used for test results. 


4. Results 


Three tests, in which weight rate of air flow, W’, 
was varied, were carried out on the model under- 
ground tunnel. The parameters for these tests are 


given in table 1. 


TABLE 1 Test parameters 


Angular velocity, w 0.3927 radian/hr 


runne!l radius, a 25 ft 

Tunnel length, Z 33.7 ft 

Thermal conductivity (mortar), K 0.747 Btu/hr-ft-°F 
Thermal diffusivity (mortar), @ O31 ft2/hr 

Specific heat (air), c’ 24 Btu/lb-°F 
Z=dyw/a. SY 


are given in table 2. The 
temperatures were substan- 


Pertinent test data 
observed entering air 
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tially sinusoidal and conformed very closely to the PARLt | pplied entering and predicted eaving , air 
equations (within +0.25 °F), of the form of eq (1 uti ea 
set forth in table 3. The observed values of air 
. . 907 ¢ k I ng I 7 (eq 
temperature at the exit of the tunnel (L=33.7 ft 
were found to lag in phase those of the entering air ies eS: a 
by about 0.14 radian. All experimental quantities | 1 T=8).7-+ T.=8 
aie . 80.4 & 10.2 1457 
involved in eq (4 


were determined directly, 
the coefficient of heat transfer, h. 
maximum when (at—As’ 
for computing A’, from which the values of h were 
obtained. The values of h thus computed are pre- 
sented in table 2 and are compared with values of 
the forced convection heat transfer coefficient cal- 
culated from an equation given by MeAdams [4 
which, with appropriate open for the model 
tunnel and air substituted, 


except 
Since @, is a 
equals zero, eq (4) serves 





Substituting the derived values of h, and observed 
values of A, 6;, a the parameters from table 1 in 
eq (4), for L=33.7 ft, we obtain the time variation 
of the leaving wie tonipernters as given by the equa- 
tions in table 3. When these equations were plotted, 
the curves followed the experimental data within 

0.25 °F. 

The tempering effect of the tunnel on the air is 
shown by referring to test No. 3, table 3, where the 








h—0.016(W’)°* -) | amplitude of the air which was 20.4 °F at the en- 
‘ | trance, was reduced to 10.2 °F at the exit. 
. Although the experimental data on air temperature 
ly z 7 al ane results tron fests on node ( rqre ne 
Pasi 2. Data and , mn model underground | satisfied the analytical solution of the problem, it 
rive 
was felt that further substantiating information was 
Ni desirable. Equations for the time-variation of tem- 
peratures at various radii in the surrounding cement- 
ecisring air tems. Te * sand mortar were therefore computed from eq (12) 
poate “=== ms | es | for test No. 3, as shown in table 4. Plots of these 
Mean, T : te 87.3 80.7 s0.4 equations with associated experimental data appear 
d slitude, A = 22. 1 20.7 0. : . ‘ 
a in figure 3. 
Leaving air temps., Tx, °F -_ = 
Maximum 95.0 90. 2 91.5 | TABLE 4. Temperatures in the cement-sand mortar for test No. 
may . vb “ | he 81.35 3 at cross-section x=22.7 ft, as predicted by eq (12) of section 
Mean remote temp. of id (mortar), 7., °F 83. 2 | 80.3 | 82. 2 | 8.1 
6:= T:—T., °F hei einpeadanaed ae 1.2 | 0.4 | —18 | = ———— = = ——— 
Air flow rate, W’, lb/hr Scape 336.4. | 367.9 | 451.2 | 
| | At radius r=0.33 ft T'z.r=82.0+-3.1 cos (wt—0.733) 
Coefficient of heat transfer h, Btu/hr-ft?-°F: | At radius r=0.50 ft jt 2.1+1.7 cos (wt—1.169) 
From eq (4)..-.. pica 1.76 | 1.85 2. 22 At radius r=0.75 ft ‘ 0.7 cos (wi—1.815) 
From eq (5)..-- SE ere ‘ 169 | 1. 82 2.14 | At radius r=1.25 ft T: 2.2+0.1 
a ee —— = — | 
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FIGURE 3 Air and mortar temperature history for model underground tunnel (test No. 8, table 
Experimental data for temperatures in mortar at various radii from center of tunnel are denoted by symbo 
5. Discussion leaving air temperatures given in table 3. These 


For the model tunnel, values of the coefficient of 
heat transfer were computed from eq (4) and the 
temperature amplitudes of the leaving air by the 
method outlined in section 4. A comparison of 
these coefficients with those determined from eq (5) 
shows agreement within 5 percent, which is within 
the probable accuracy of the thermal properties of 
cement-sand mortar as determined by tests on a few 
samples. 

For tunnels or shafts probable in undergrouid 
installations, the diameters are likely to be greater 
than 3 ft. Values of the coefficient of heat transfer 
are not found in the literature for airways of this 
size and would be particularly difficult to predeter- 
mine for the rough and irregular surfaces usually 
resulting from the excavating process. However, 
reasonable values may be estimated from compro- 
mise considerations of values determined for air flow 
in small cylinders and air flow parallel to planes, 
and with an allowance for roughness of surfaces. 

Substitution of the parameters and the derived 
values of A into eq (4) gave the equations for the 


| 


equations showed differences less than 0.3 
plotted with the experimental data. Differences 
less than 0.5 °F were found between observed tem- 
peratures in the surrounding solid (mortar), and those 
computed by eq (12), as shown in figure 3. 

Because the experimental data agree closely with 
values computed by eqs (4) and (12), the analytic 
approach developed in this paper should be appli- 
cable for determining temperatures in tunnels of 
larger size. 


IF when 


6. Application to a Practical Problem 


To design heating and air conditioning equipment 
for an underground installation, it is necessary to 
predict the maximum and minimum temperatures 
of the air when it is delivered at the installation. 
From eq (4), the maximum and minimum tempera- 
tures occur when wf=ABP’ and wt=—7r+ABP’, respec- 
tively. Using eq (4), these extreme values may be 
written 

(0,) r==e 44’ (4+A4 


6; cos AB’). (6) 
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Figure 1 is a plot of A’ and B’ as functions of n and 
z computed for a range of values expected in a 
practical problem. 

For a sample hypothetical problem, consider a 
square tunnel with nominal dimensions 7 ft by 7 
ft by 2000 ft long, continuously delivering air to an 
underground installation at a rate of 50,000 cubic 
feet per minute, or 228,000 pounds per hour. For 
noncircular cross sections, the equivalent cylindrical 
radius, a, is equal to 2 (cross sectional area) /(per- 
imeter). Assumed parameters for substitution in 
eq (6) are given in table 5. 

TaBLe 5. Pa 


rameters fo 


/ spotl elical proble m 











Thermal conductivity of surrounding rock K=1.2 Btu/hr-ft-°l 
rhermal diffusivity of surrounding rock a=(0),035 {t2/hr 
Angular velocity (period one year w=0.000717 radian/hr 
Coetlicient of heat tra h=3.43 Btu/hr-ft’-° 
Equivalent radius of tunne a=3.5 ft 
Amplitude of outside air A=45 °F 
Mean ten ure ke Ti=52 ° 
Remote ro mip T K 
Specific heat of air c 24 Btu/lb-° 
Weight rate of air flow Ww’ 238.000 Ib hr 
7] T;—T 3 

ah/ kk 10 

Lyw/a=0.5 
1’ (funct of l figure 4.48 
B funct ! i figure | 215 
1 Kl i 0.0439 


The extreme values predicted by eq (6) are 


(0,), 34 and —39 °F. 


Because @ is a subrogated temperature based on a 
zero Value equal to the remote temperature of the 
surrounding rock, the maximum and minimum air 
temperatures at the end of the tunnel (22000 ft) 
are 554+-34=89 °F and 55—39=16 °F, respectively. 
The maximum and minimum outside air tempera- 
tures (eq (1 are 97 °F and 7 °F, respectively. 
The maximum and minimum air temperatures at the 
end of the tunnel lag the entrance by an 
angle equal to AB’, which in this case equals 0.094 
radian, equivalent to a time-lag of 5.5 days. 

Although the rates of heat loss and gain of the air 
in a tunnel were not in the experimental 
portion of this paper, they are of considerable interest 
in the application to practical problems. For the 
hypothetical problem, the average rate of cooling 
by the rock during the cooling season (eq (14), 
tion 8) is 251,000 Btu/hr, or an equivalent of about 
21 tons of refrigeration. The average rate of heating 
by the rock in the heating season is 309,000 Btu/hr, 
equivalent to the heating effected by burning about 
3.35 gallons of fuel oil per hour. The maximum 
cooling and heating rates are substantially greater 
than the average, as indicated by eq (13) 


those at 


cliscussed 


sec- 
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7. Conclusions 


Mathematical equations are derived for computing 
the tempering effect of a long tunnel in reducing the 
amplitude of climatic temperature variation of out- 
side air flowing through a tunnel. Results of experi- 
ments support the application of the developed 
equations for design purposes 











The functions necessary for design use of the 
equations have been evaluated for a useful range of 
the parameters and are presented graphically. 

In the design for heating and air conditioning sys- 
tems, fresh or outside air necessary for ventilation of 
an underground installation represents a portion of 
the heating and cooling load for the installation. 
Without a means for determining the temperature 
- the air at the end of “ passage through tunnels 

shafts, the designer is forced to rely upon the 
ean air temperatures for design conditions. The 
hypothetical problem of section 6 shows that the 
heating and cooling loads for design purposes are 
considerably reduced when the temperature of air at 
the end of a tunnel is predicted by the derived equa- 
tions. These equations are included in a design 
manual for military installations [5]. 

An assumption has been made which develops a 
degree of uncertainty toward the use of the equations. 
An effective coefficient of heat transfer value ‘h’’ 
— for a particular application may, in most 

rases, only be determined as an educated guess, be- 
cause no predictable values for large size tunnel 
airways are available from literature. 


8. Appendix 
8.1. Solution of Heat Transfer Problem 
With the boundary conditions as set forth in eqs 
(1), (2), and (3) of section 2, eq (4) may be obtained 
by assuming that the temperature in the surrounding 
solid at a cross section at a distance z from the en- 
trance 1s 
Coie Fe'*', where F f(r) (7) 
and also that the temperature of the flowing alr is: 


6, Ge ade where G fi | (8) 


By substituting the appropriate derivatives of ed (7) 
in the differential equation for conduction heat 
radially in a eylindrical system having peo 
svmimetry: 


O°A ! Oo¢ 1086 
or ror adt 


the latter becomes 
Pr 4 ) ages -—=(). 


This is the modified Bessel equation of zero order, 
for which the general solution given in terms of the 
modified Bessel functions I>, Ko, is 


om (- \ ‘s) 4 K,(r/*) 


For F to be finite as 7 approaches infinity, D=0, 
and a particular solution is: 
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Pax Ko(ry/ ~\ 


and from eq (7), 


6, 2=C'e Ko (ry = ): (9) 
Qa 


Differentiating, 


d 
(ir), 


Substituting eqs (9) and (10) in eq (3), and solving 
for C, eq (10) becomes: 


: is lw > lw 
te** \ : Ki(ay - ): (10) 


a 


h /w,, ie 
i) * wik,. K\ a Ki (ay a 


j : j : (10a) 
ar/);=a 2 es) (® (W 
2 a ri 'K ) 
Va kilay st RK ays 
Simplify by letting 
avw/a 
Ky (2v1 N(z)¢ 
"s = f 27) 
yi Ky (ey? Ni (2) eX? Tas 
ah K 
Ni 
P=z=> 
Not 
+ v { 
re] Q) | 4 rea) 
Where values of the functions of N,(< N(2), Dy 2 


and @;(2) are civen by MeLachlen and Mevers [6], eq 
(10a) becomes 


oe) 6,) nPetet*/2 


dr a n+Pe press 


Rationalizing and further simplifying, eq (3) be- 


comes 


K( =) Kké, P—n sin ut in COS mM 
my ; 
di a Vi P ° 
pt —2 sin yp 
i] 

K@ , 12D 

a (A’+1B’)=h(0,—9,, 2) 

aT 


where 


2r(P nmsin pw 


A’ » (fig. 1) 
’ o. 
aw 
" + P «- SID yw 
9 
, 271i COS Ut . 
B P ? fig. | 
| “ -2 sin u 
n 


| Substituting the above in eq (2) 


v do, . de, Ke, 


a *__ (4’+2B’)=—0. 
dx‘ dt T2oeM ca 


| Substituting from eq (8), 


GK 


VG’ +i0G-4 (4 
2nrM ca 


i’+7B’)=0. 


Kx Kx 


Letting A SaV Me? We? 


the solution of the differential equation is: 


Substituting in eq (9), 


0,=JS | it = A(A’ iB’) | 


Equation (1) states 
ther 


Making this substitution, 


— | ae yp AB 4+ 6. (+40) 


Considering only the real part, 


— “| acos(o = AB’ )+6 cos (7+-4B") | 


For frequencies w under consideration, the term wr/V 
is negligible and a closely approximate solution is 


6,=e-“44'[A cos (wt—AB’)+6; cos AB’). (4) 


Without giving the derivation from eq (9), the 
temperature at radius, 7, 
the tunnel is: 


in the solid surrounding 


‘a 
: cos (wt+S—r)-4 cos (S—r) 
1A N,(R) 
Ae = 
No! P\2 » 
( ) +1—2—sinug 
\ n n 
(12) 
where 
if Py w/a rea 
S rea) (R Do\ } AB’ 
COS UL 
7=— tan 
fi ° 
7 SIN pb 
/ 
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8.2. Air Heat Exchange 


At any time, ¢, the rate at which heat is lost by 
the flowing air is: 


q We’ (4 0,). 


The maximum rates of cooling and heating of the 
air are given by 


cos AB’ 


6,(1 e~ 44’ eos AB’)). 


J 


If it is assumed that cos Ab’~1.0, then 
(1 hen A-+6;). (13) 


The average rate of heat 


siven cooling 
or heating evele Is: 


flow for a 


(2>++.): 


8.3. Symbols 


(14) 


Unless otherwise labeled, quantities 


dimen- 
sionless when consistent units are 


are 
used: 


A=—K2/W'e’ 


27 pP “silk uw 


A! =- » plotted on figure 1 as a function 
/ / P | 
! 9 sin 5 of n anc 
os 
a—radius of hollow evlinder, ft, 
2(cross section)/(perimeter), for noncircular 
evlinders, 
- Irn COS LU : ’ F 
b P » plotted on figure las a function 
dt 9 sin u of n and 
aoe ; 


c=specific heat of surrounding solid, Btu/lb-°F, 
c’ specific heat of air flowime in 
stu lb I, 
e=base of Naperian (natural) logarithms, 
h—coefficient of heat transfer between 
surface of evlinder, Btu/hr-ft?-°F, 
K=thermal conductivity of 
Btu/hr-ft-°F, 
M’=weight of fluid (air 
ot evlinder, Ib /ft?, 
n ah K, 
P=2N,(z2)/N ; 
q=heat flow rate from air to tunnel wall, Btu/hr, 
r=radius, {t,! 


evlinder, 


air and 
surrounding = solid, 


per unit of surface area 


T=temperature, °F. 7) is temperature of air 
entering tunnel; 7; is annual mean value of 
T,. T, is the steady temperature of the deep 
underground rock remotely surrounding the 
tunnel, 

t=time, hr, 

V= velocity in tunnel, ft/hr, 

W’=weight rate of air flow, lb/hr, 
x=distance from entrance of tunnel, ft, 


(ly W@/ a, 
a=thermal diffusivity of surrounding - solid, 
ft?/hr, 


A=amplitude of annual variation of 7) about 


the mean value 7,, °F, 


\ T) max T’;) 1 /2( T max T min 
MQ (2 w/4 DMo\Z), 
§—subrogated temperature, equal to (7— 7,), °F, 


#,—subrogated annual mean entering air temper- 
ature, equal to (7,;— 7,), °F, 


w=angular velocity, radian/hr. 


Temperatures at various positions are indicated by 
the subscripts: 


Q—at tunnel inlet (r=0, ra), 
gat distance z from tunnel inlet (r<a 
L=at tunnel outlet (x L,? a) 


r.e—at radius 7 at distance x from tunnel inlet 
(roa 

a—at tunnel surface (r=a) at distance « from 
inlet. 


The subscript m denotes a maximum or 


minimum 


value 
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3oth the incubation period prior to stress-corrosion cracking and its initiation in the 
AZ31B magnesium alloy were studied using electrochemical solution potential-time and 


extension-time oscillographic curves 


obtained 


simultaneously. From these data it is 


postulated that the crystals most favorably oriented for slip were deformed, on initial stress- 
ing of the specimen, to such an extent that they were resistant to further extension and the 


protective films on their surfaces, destroyed in the early deformation, were repaired. 


Sub- 


sequent extension occurred as the result of stresses concentrated in relatively more resistant 
crystals with strain rates as high as 1 percent/see at limited sites rupturing the protective 


surface film. 


This exposed metal was anodic to very much larger areas on which the film 


had been repaired and cracks developed by electrochemical processes. 
Stress-corrosion cracks were initiated on planes that made large angles with the basal 
crystallographic plane of the magnesium alloy and in crystals of a polycrystalline aggregate 


that were unfavorably oriented for basal slip 


1. Introduction 


The stress-corrosion cracking of magnesium-base 
allovs has been the subject of several studies at 
NBs. A study of the relative susceptibility of the 
various wrought magnesium alloys then available 
was made during World War II and the results were 
subsequently published [/].' A service failure that 
led to the initiation of that study proved to be an 
isolated case. Stress-corrosion cracking of magne- 
service not been a problem. 
However, magnesium alloys lend themselves par- 
ticularly well to ceneral studies of stress-corrosion 
cracking in that complete failures, for example, of 
the FS-1 alloy, can be obtained in the laboratory in 
less than 5 min. Studies since the war vears have 
been made with the objective of increasing our 
understanding of the mechanism of stress-corrosion 
cracking. 

A paper published in 1958 [2] gave the results of an 
investigation initiated to answer certain questions 
that had been raised but never answered during the 
war work. In that investigation [2] it was shown 
that if the specimen was strained in the corrodent at 
a rate in 500 107°/min, measured 1 min 
after the stress * was applied, failure could be ex- 
pected in less than 5 min. On the other hand if the 
strain rate was less than 500 10~°/min the specimen 
would most probably not fail even in 1000 min. 
Further, it was found on examining many specimens 
that there was no evidence of cracking for the first 
40 sec after the stress was applied even though 
failures occurred in 2's min. A short hypothesis was 
presented to explain this incubation period. In the 


NN 


sium alloys in has 


eXCeSS of 


Italicized figures in brackets indicate the literature references at the end of 
this paper 
2 A stress in excess of 95.5 percent of the yield strength of the annealed material 
was required to produce stress-corrosion cracking of this alloy in the corrodent 
used 
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present paper, experiments designed to verify some 
of the ramifications of that hypothesis are described 
and the hypothesis is expanded. The crystallo- 
graphic planes on which stress-corrosion cracks 
originated were also determined. 


2. Experimental Procedure 


Specimens for this investigation were prepared 
from sheet and extruded AZ31B magnesium alloy 
stock. The chemical compositions of these materials 
were well within the tolerances of ASTM Specifica- 
tion B-90-58 for this alloy. 

Sheet tensile specimens, having a '4-in. reduced 
section, were machined from ‘\¢ and 's-in. thick sheet 
material originally obtained for the earlier work [1]. 
Most of the specimens used in the present investiga- 
tion were in the annealed condition although a few 
were from “hard-rolled”’ sheet. 

A very few notched specimens were machined from 
-in. diam extruded rod. The specimens contained 
a 90° notch 'sin. deep. The diameter at the root of 
the notch was 0.250 in. and the root radius 0.005 in. 

Specimens used in determining the crystallo- 
graphic planes on which the stress-corrosion cracks 
originated were also machined from the extruded 
rod. The reduced sections of these specimens were 
about 2 in. long with a square cross section 4% in. on 
aside. They were similar in design to the low carbon 
steel specimens used by the author in another in- 
vestigation [3]. After machining, the specimens were 
annealed for 30 min at 400 °C, cooled to room 
temperature, strained 0.3 to 0.5 percent, and sub- 
sequently heated slowly to 520 °C. The resulting 
specimen was multigrained, with individual grains 

’ The nominal composition of the AZ31B alloy is Al 3.0 percent, Zn 1 0 percent, 


Mn 0.20 min, with maximum values for Si, Cu, Fe, Ni, and Ca of 0.10 percent, 0.05 
percent, 0.005 percent, 0.005 percent, and 0.04 percent re spectively. 








sufficiently large that their crystal orientations could 
be determined by the conventional back reflection 
X-ray technique. 

The sheet specimens were insulated from the 
corrodent, except for a 2 in. gage length, by means of 
plastic tape. All of the notched specimens were 
also insulated with tape except over the notched 
areas. One set of the notched specimens was exposed 
to the corrodent with the entire notched area exposed. 
In the other set the notch was covered with lacquer 
and after the lacquer had dried it was cut through 
with a razor blade to expose the metal at the root of 
the notch. 

All specimens were contained in glass cells of the 
type described earlier [/]. The corrodent was the 


conventional NaCl, K,CrQ, solution (35 ¢ NaCl 
and 20 g K,CrO,/liter) at approximately 30°C. The 


specimens were stressed in tension by dead loading 
with a lever system. In all cases care was exercised 
to avoid impact in applying the stress. If the spec- 
imen strain was to be followed, the lever was sup- 
ported by a hydraulic jack until the load, etc., were 
adjusted. The jack was then released, applying the 
stress at a rapid rate. 

The strain or extension of the specimen was fol- 
lowed by measuring the movement of the end of the 
loading lever (lever ratio 25 to 1). This was done by 
reading a dial gage (smallest division 1/10,000 in.) 
or by recording the unbalance of a bridge circuit, 
one arm of which contained wire strain gages actuated 
by the motion of the end of the lever. In a few 
instances wire strain gages cemented onto the spec- 
imen itself (but insulated from the solution) were 
incorporated in the bridge circuit. 

The electrochemical solution potentials of the sheet 
and notched specimens were determined against a 
calomel electrode of the saturated KCI type. Con- 
nection between the electrode and the corrodent was 
made by means of a salt bridge. The electrochemical 
solution potential and the extension of the specimen 
(as determined by the unbalance of the bridge 
circuit) were recorded simultaneously on a_ two- 
channel oscillograph. 

No extensions nor electrochemical solution poten- 
tials were determined for the large-grained specimens 
on which the crystallographic orientation of the plane 
of fracture was determined. The was re- 
moved from some of these specimens after cracks 
had developed but before the final fracture had 
occurred. If fracture was complete, no attempt was 
made to determine the plane on which this fracture 
occurred, but analyses were made to determine the 
planes of origin of other cracks that formed but had 
not led to complete failure. 

Visual observations showed that cracks generally 
originated at a corner of the specimen where adjacent 
surfaces met. Hence, the crystallographic planes 
on which the cracks originated were determined from 
the angles that the cracks made on the adjacent faces 
with the corners of the specimen [4]. 

The orientation of the tension axis of the specimen 
with respect to the crystallographic axes was deter- 
mined for the grains containing cracks and also for 
a number of grains free from cracks. 


stress 
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| 
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3. Results and Discussion 


3.1. Extension Time and Electrochemical Solution 
Potential Curves 


The idealized curves shown in figure 1 were based 
on a study of simultaneous oscillographic records of 
the strain and accompanying changes in electro- 
chemical solution potentials for a number of speci- 
mens. The extension increased very rapidly upon 
initial loading of the specimen. It continued to in- 
crease but at a diminishing rate for 40 to 50 sec. 
The rate then increased again, with the development 
of stress-corrosion cracks, to final failure of the speci- 
men after about 150 sec exposure. The electrochem- 
ical solution potential changed about 90 my in the 
anodic direction on the initial loading of the 
specimen. This was attributed to the rupturing of 
the protective film on the surface. The further 
marked change of potential in the anodic direction 
with the exposure of bare metal, on the failure of the 
specunen, substantiated this idea. The total change 
in potential from the unstressed condition amounted 
to approximately 200 my, a value that is in good 
agreement with previously published data [4]. It 
should be noted that once failure had occurred and 
the metal was ho longer being strained the potential 
changed very abruptly, probably within 1 sec, to a 
value only about 20 mv more active than that of the 
specimen in the original condition. The extremely 
rapid rate of change in potential after the fracture of 
the specimen indicated that a protective film was 
being formed on the surface of what bad been bare 
metal. 

The yield strength of the annealed allov can be 
increased about 50 percent by cold work. From the 
extension curve, figure 1, it is to be expected that the 
alloy is being cold worked. It is postulated that 
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with the application of stress, those grains most 
favorably oriented with respect to the stress pattern 
were deformed by basal slip, kinking, etc., and were 
hardened to the point where further deformation 
occurred only very slowly. Specifically, dislocations 
moved through the material until they were blocked 
at barriers such as grain boundaries, subgrain bound- 
aries, inclusions, ete. Stress concentrations were 
built up as a result of this dislocation pileup and slip 
then occurred in some of these grains that had pre- 
viously resisted deformation. This process would be 
repe ‘ated as stresses concentrated in suc cessively 
more resistant grains and produced deformation. 
The delayed strain in these relatively resistant grains 
thus accounted for the slow increase in extension 
after the initial deformation. While the deforma- 
tion of individual grains was postulated to be a 
discontinuous process, the average effect over many 
grains was recorded as a continuous process. 

The increase in the potential (in the active direc- 
tion) on the initial loading of the specimen is postu- 
lated, was indicated above, to be due to the 
rupturing of the protective film on those grains that 
were most readily plastically deformed. 

The data obtained at the final fracture of the 
specimen would indicate that as a given grain was 
work hardened to the point where it became resistant 
to further strain, where the strain rate was low, 
the protective film was repaired on its surface. As 
the ratio of grains being strained at a rate sufficient 
to rupture the film decreases with respect to those on 
which the film had been repaired, the potential would 
change in the noble direction. This would account 
for the change in potential, after the initial extension 
of the specimen. If the average extension rate of the 
specimen, determined 40 or 60 sec after the stress was 
applied, was below some critical value, the potential 
would continue to change in the noble direction until 
it approached the value obtained prior to stressing 
the specimen. 

If the average strain rate exceeds this critical value, 
stress concentrations at a few sites will deform the 
specimen surface at a sufficient rate to prevent the 
film from reforming. If these sites are localized for 
a sufficient time, cracks will be initiated by an elec- 
trochemical process. As was indicated earlier, the 
bare metal is at least 200 mv anodic to the film- 
covered surface in the corrodent used. If a relatively 
large area is film-free, attack will be of the pitting 
type. If the film-free area is very small and is sur- 
rounded by a large film-covered area, the condition 
of a small anode and large cathode will exist; this is 
an ideal condition for rapid corrosion of the anode. 
Such corrosion on a small area will form a sharp pit 

r trench that will further concentrate the stresses, 
elaine the film from reforming at its apex, and 
will lead to the initiation of stress-corrosion cracks. 


as 


1.€., 


Once cracking is initiated, the continuous exposure 
of bare metal at the tips of the cracks will change 
the potential of the specimen in the anodic direction. 
A very marked change occurred, as was pointed out 
earlier, with the exposure of a considerable surface of 
bare metal at the failure of the specimen. 


3.2. Strain-Rate Associated With Initiation of Cracks 


The extension of many specimens was followed by 


readings, at stated intervals, of the dial gage in 
contact with the end of the lever system. This 
method was highly sensitive and although it was 


not possible to measure the extension of the specimen 
during the first 10 to 20 see of exposure, data was 
obtained subsequently from which extension rate 

time curves could be drawn. Three such curves 
are shown in figure 2. Breaks in curves A and B 
resulted from the development of stress-corrosion 
cracks in those specimens. No cracks developed 
in spec ~ n C which extended at a slower rate than 
either A or B. Cracks began to develop in speci- 
mens fees at the approximate rate shown for 
B after about 40 sec exposure in the corrodent [2]. 
If the strain at sites where cracks will develop is 
postulated to account for the difference in the strain 
rates of specimens B and C after 40 sec, if certain 
assumptions are made as to the width of the regions 
in which cracks will develop, and if the number of 
cracks that develop per inch of specimen length is 


determined, an estimate can be made of the strain 


rate at these crack sites. 


It was found that approxi- 
mately 


25 crack sites develop for each inch of gage 
length. Assume that the strain is confined to : 
region 1000 times the ‘“‘a’’ parameter of the mag- 
nesium lattice; then the average strain rate in the 
areas where cracks were being initiated can be 
calculated to have been approximately 1X 10~2/see 
(1%/sec). This value seems reasonable within one 
order of magnitude. It is comparable with the value 
of 5 percent/see reported by Hoar and West [6] 
the strain rate necessary to prevent polariza- 
tion that would have stifled development of cracks in 
austenitic stainless steel specimens. 
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FIGURE 2. Relations between time and strain rate for specimens 


stressed above and below the critical strain rate at which stress- 
corrosion cracks develop. 


Breaks in curves A and B indicate the development of stress-corrosion cracks 
No cracks developed in the specimen from which curve C was obtained 
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3.3. Effect of Ratio of Cathode to Anode Areas on 
Incubation Period 


Several notched specimens were exposed to the 
corrodent at a stress of approximately 46,000 lb/in.’; 
this stress * could produce failures by stress-corrosion 
cracking in 60 sec or less. As was indicated above, 
all of the metal in the notched areas of some of these 
specimens was exposed to the corrodent while only 
the metal at the root of the notch was exposed in 
other specimens. 

All specimens failed at the stress level used. How- 
ever, the incubation period (exposure period before 
there were indications that cracks were initiated) 
was about six times as long for specimens that were 
completely insulated except at the root of the notch, 
as for those specimens where the whole notched area 
was exposed to the corrodent. 

The concentrated at the roots of the 
notches of these latter specimens limited all strains 
of any appreciable magnitude to those regions. 
The protective film will be ruptured only at the root 
of the notch. The sides of the notch (when exposed 
to the corrodent) will be film-covered, cathodic to the 
root of the notch, and very much larger than the 
anodic (film-free) area at the root of the notch. 
Thus there exists the condition of a very small anode 
and large cathode that is ideal for 
by electrochemical processes. 

In specimens on which only the roots of the notches 
are exposed to the corrodent there will be some small 
areas where the strain rate is higher than others, 
that is, on which the protective film will be ruptured. 
Conversely, there will be some areas where there will 
be more opportunity for the film to repair than on 
others. The resulting anodic and cathodic areas 
may be nearly the same size and while electrochemi- 


stresses 


rapid corrosion 


initiated 
eal corrosion may occur and cracking eventually 
i ‘velop, the rate will be much slower than if a large | Taste 1. Components of stress normal to the plane of stress- 
cathode-small anode relationship exists. corrosion cracks, and probability of ip ote ages gle 
The evidence [7, 8] indicated that once cracks 19, wh «ei hi ap : “eile ae a nt ee 
have penetrated into the metal, the sides of the ipiabistnsualbis cient e em ey 
cracks become the cathodic areas; whether or not ft Sie aad the slip pune and lip direction, respectively. 
the sides of the notches are insulated is then of little 
further importance. Angie between | Component of Probability of 
Crystal ion axis and the tensile stress) slip at con- 
sion erack |corrosionerack | 
3.4. Crystallographic Data on Crystals Containing, | - 
and Free From, Stress-Corrosion Cracks deg sin d 008 d 
on 12 0. Y7s 0.000 
C2 5 909 OOS 
The poles of the planes on which stress-corrosion | hy 4 a oa 
cracks were initiated are shown in a single stereo- | CE 10 ws 043 
graphic triangle in figure 3. In every case these 4 16 961 vr 
planes make large angles with the basal plane. If C8 ‘ - rd: 
the same condition exists in this alloy as has been C 10 W 5 078 
reported for relatively pure magnesium [9], namely, | C12 ti 966 M4 
that the basal planes of the cry ‘stals in an extrusion | a _ 
are approximately parallel to the direction of ex- N2 oo 
trusion, then the fact that the poles of the stress- | a 130 
corrosion cracks make high angles with the pole of | N 5 | 
the basal plane is to be expected. N7 - 342 
ye “408 
‘ Stresses of this magnitude will produce a true strain of approximately 0.03; sn wi ~_ 
this included a deviation from the modulus line of about 0.015. 
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With the exception of C-2, C-8, and possibly C-11, 
the planes of fracture do not lie close to any low-index 
crystallographic planes. These results agree in a 
general way with those obtained by Meller and 
Metzger [10] on single crystals of an AZ61 allov and 
do not agree with those reported by Priest [7] for 
the same alloy. 

It is well recognized that the planes of stress- 
corrosion cracks are generally approximately normal 
to the stress axis. The data obtained in this in- 
vestigation further confirm that observation. The 
resolved tensile stress normal to the plane of the 
stress-corrosion crack is proportional to the cosine of 
the angle between the tension axis and the normal to 
the plane of the stress-corrosion crack. These angles 
for the various stress-corrosion cracks are given in col- 
umn 2 of table 1. The corresponding cosines are list- 
edin column 3. The component of the stress normal 
to the stress-corrosion crack is in no case smaller than 


0.934 times the tensile stress. 


(O110) 








(0001) (1120) 


FIGURE 3. Poles of planes on which stress-corrosion cracks were 











(O00!) C7. 


FIGURI }. Tension axes for crystals “C” in 


j 


which crac k ing 
occurred and for crystals ‘*N”’ in 


which no cracks were 


found. 


The ease with which slip will occur in a given 
crystal Is dependent, in part al least, on the magni- 
tude of the product sin @ cos A, where ¢ is the angle 
between the stress axis and the slip plane and 
is the angle between the stress axis and the slip 
direction. The slip plane is (0001) and the slip 
direction is [1120] for polyerystalline magnesium [9] 
and it is to be expected that they are the same for the 
AZ31 alloy. The values of sin @ cos \ have been 
calculated for each of the 12 crystals in which cracks 
were analyzed and for 10 crystals, 
same specimens, in which no cracks were found. 
These data are given in column 4 of table The 
tension axis, with respect to each crystal, is plotted 
in a single stereographic triangle, in figure 4. 

Data in table 1 and figure 4 indicate that crystals 
in which cracking occurred were generally less 
favorably oriented for slip than those crystals ana- 
lyzed in which no cracks were found. — It is recognized 
that in polverystalline specimens the neighboring 
crystals may be the determining factor in whether or 
not slip occurs in a particular crystal. However, 
the fact that the stress-corrosion cracks occurred 1 
the crystals least favorably oriented for slip is in 
agreement with the mechanism postulated earlier for 
the initiation of cracking in the 
magnesium alloy. 


from some of the 


stress-corrosion 


4. Summary 


1. Idealized curves were prepared showing the 
simultaneous changes in extension and_ electro- 
chemical solution potential for AZ31B magnesium 
alloy specimens that were exposed under constant 
tensile load in an aqueous solution of NaCl and 
K,CrO, until failures by cracking 
occurred. 


stress-corrosion 








2. Initiation of stress-corrosion cracking, after an 
incubation period, is postulated to occur in grains 
that were deformed as the result of stress-concentra- 
tions built up by the initial deformation and the 


hardening of those grains that were favorably 
oriented for initial plastic extension. 
3. Making several assumptions, the strain rate 


(in microscopic sections) associated with, and most 

probably required for, the initiation of stress-cor- 

rosion cracking was calculated to be of the order of 
percent/sec. 

The incubation period (period prior to initiation 
of cracking) was much longer for specimens on which 
the passive (film covered) and active (film free) areas 
were approximately the same than for those on which 
the passive area was many times larger than the 
active area. 

5. Stress-corrosion cracks made large angles with 
the basal crystallographic plane in large-grained 
specimens but were not generally associated, at their 
initiation, with any low-index crystallographic plane. 

6. Stress-corrosion cracks were generally approxi- 
mately normal to the stress axis and occurred in 
grains unfavorably oriented for slip. 
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Corrosion Research Council and the National Bureau 
of Standards. The author appreciates the assistance 
of Joan Calvert, June Chung Fugh, James G. Early, 
Jr., Joseph E. Himes, and Howard T. Yolken in the 
experimental work in this investigation. 
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Three steel specimens were continuously exposed in the laboratory for 
by 


in city water to which 
specimens were under 


added 3 percent 
continuous cathodic 


was 


and the other by current from a rectifier through 
left to corrode freely. 
As a result of the cathodic protection, 


1961) 


almost 5 vears 
weight of sodium chloride. Two of these 
protection, one by current from a zine anode, 


a carbonanode. The third specimen was 


carbonates and silicates formed protective 


coatings which eventually reduced the current density required for protection from about 


5 to 0.02 ma/ft \ coated specimen, 
including 12 days out of the salt water), 


while without protective currents were 


zation curves The currents 
be related to AV/A/J values from the curves 
eorrosion rates as measured by weight loss, 
1. Introduction 
When cathodic protection is mentioned in con- 


junction with coatings for protecting steel structures 
against corrosion, the coatings referred to are in- 
variably organic in nature. Relatively little is said 
of the benefits of the inorganic or natural coatings 
such as carbonates or silicates which can be deposited 
result of cathodic protection and which the 
authors have found to be quite corrosion-resistant in 
a salt water environment. Although the protective 
character not necessarily in con- 
nection with cathodic protection, has been emphasized 
by other investigators it is believed that the present 
discussion will be of additional value and interest 
because the very nature of cathodic protection when 
properly applied is such as to automatically keep 
these coatings in a state of re pair, 

Speller [1] ' states that silicious and carbonaceous 
coatings are sometimes formed on metal surfaces 
by adsorption in the atmosphere, underwater, and 
in soils, and often act as the main natural agency 
opposing corrosion. 

Evans [2] in talking about the behavior of zinc- 
coated steel continually immersed in water has 
suggested that the carbonate electrochemically 
deposited on the steel, when the zine into 
solution, may persist after all the zine has 
disappeared. 

In some of the early NBS soil corrosion tests, 
Logan [3] attributed the protective effect of the zine 
coating on galvanized pipe specimens to a protective 


as a 


of such coatings, 


goes 


even 


nr 
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indicate the literature references at the end of this paper. 


after 
required only 0.3 ma/ft 
rhe instantaneous corrosion rates on the coated specimens (ser: 


measured 
at which breaks occurred in the 
which 


eing without protective current for 32 days 
for initial protection. 

atched and unscratched) 
(breaks) in polari- 
curves were found to 
bore a relation to the 


by changes-in-slope 


cathodic 
values in turn 


film found on the surface of the zine rather than in 
cathodic protection by the zine. Later, Denison 
and Romanoff, [4] in reporting on the behavior of 
galvanized steel specimens exposed to corrosive soils 


for 13 yr, tentatively attributed their remarkable 
corrosion resistance in alkaline soils to a film or 
coating, probably silicious in nature, that was 


deposited cathodically by galvanic action between 
the remaining zine coating and the alloy layer on the 
underlying steel. The protection of the steel was 
obvious even though the outer zine coating had been 
virtually removed by corrosion. 

Hoar [5] in a private communication suggested 


that the protection produced by the cathodic 
deposition of calcium carbonate on a metal surface 
can persist in many cases long after the cathodic 


current is stopped. 

Pearson [6] associates the useful ampere-hour 
effect of cathodic protection on ferrous materials 
underground with a film formed by the precipits ation 


of insoluble carbonates. He describes the “film” 
formed as being effective even when barely visible 
and building up with time. 

This paper describes the effects observed by 


cathodically protecting cold-rolled steel rods exposed 
to city water to which was added 3 percent by weight 
of sodium chloride. It is shown how the effective- 
ness of the coatings formed on the steel surfaces can 
be evaluated by the means of polarization measure- 
ments. The applied currents associated with 
changes-in-slope of the current-potential curve during 
cathodic polarization are shown to be related to 
the corrosion caused by physically disturbing the 
uniformity of a coating. The study extended over 
a period of almost 5 yr. 
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2. Experimental Procedure 


The experiments to be described started out as 
a simple demonstration of cathodic protection for 
students spending their vacations working on various 
projects at NBS. A }4-in. diam rod of cold-rolled 
steel was coupled to a zine rod through a milli- 
ammeter and exposed to Washington City water 
to which had been added 3 percent by weight of 
sodium chloride. The steel rod was protected by 
tape at the water line and had approximately 0.1 ft? 
of surface (9 in. of its length) exposed to about 17 
liters of the salt water contained in an open Pyrex 
jar. Henceforth, all currents will be expressed as 
apparent current densities in ma/ft? or ma/dm?’. 
During initial exposure, the galvanic current rapidly 
diminished from 200 ma/ft? to 10 ma/ft? and after 
3 weeks of exposure was down to 6.0 ma/ft? at which 
time the rod was removed from the water, washed 
under running water, and examined. A_ white 
chalklike coating which could readily be scratched 
off was noticeable on the lower end of the rod. The 
rod was put back into the water and again coupled 
to the zinc. As time went on, the galvanic current 
became smaller until, after 1 yr of exposure, it was 
0.5 ma/ft?. During the exposure period it was 
observed that the potential of the steel, with reference 
to a saturated calomel half cell (used throughout this 
investigation), was about —1.0 v when read immedi- 
ately after breaking the galvanic circuit. After the 
circuit had been broken for 1 min, the potential 
of the steel was still between —0.85 and —1.0 v. At 
the end of the year, the rod was uniformly covered 
with the chalky deposit (2 to 3 mils thick) containing 
zinc silicate as shown by the X-ray diffraction powder 
pattern analyses. Figure 1 is a photograph of the 
rod showing the white coating. On this occasion, 
no rust was in evidence, anywhere, even under the 
tape laps. 

Upon the conclusion of the foregoing experiment 
3 additional weighed steel specimens, prepared as 
the one 5 ‘viously described, were individually 
exposed in 3 Pyrex jars each cont: aining 17 liters of 
city sates to which had been added 3 percent by 
weight of sodium chloride. One specimen, No. 1, 
was left to corrode freely and act as a control. 
Specimen No. 2 was connected galvanically to zine 
and specimen No. 3 received protective current 
from a }s-in. diam carbon rod fed by a copper oxide 
rectifier. Here again, the specimens were set up 
primarily for demonstration purposes, but protective 
currents were measured from time to time as were 
the potentials of specimens Nos. 2 and 3 under 
cathodic protection. The water level in the open 
jars was never permitted to drop below the tape 
line, makeup city water being added weekly to 
compensate for evaporatior: loss. 

After 12 months the 3 specimens were removed 
from the jars, cleaned and reweighed. The specimens 
were put back into the same jars and electrolyte 
from which they had been removed. The specimens 
were not removed from the electrolyte again, except 
for brief periods, until the 56th month and finally 

















Inorganic 


FIGuRE l. coating formed on steel cathodically 
protected by zinc for one year in a 3 percent sodium chloride 
solution. 

for chemical analysis of the coatings formed on 

the specimens under cathodic protection, and for 

cleaning and reweighing after 57 months of exposure. 

Until the 55th month, specimens Nos. 2 and 3 
were continually under cathodic protection except 
for relatively brief periods which will be discussed 
later. The potential of specimen No. 2 (connected 
to the zinc) was always well within the protective 
range because of the good conductivity of the 

electrolyte, while the potential of specimen No. 3 

was always at least at —0.77 v and usually more 

electronegative. After 34 months exposure, a po- 
tentiometer was shunted across the secondary of the 
rectifier supplying specimen No. 3 and the applied 
voltage adjusted to a value consistent with main- 
taining a minimum protective potential of approxi- 
mately —0.77 v. By this time, a uniform natural 
coating was in evidence on specimen No. 3 and it 
was found on further reductions in the current applied 
to No. 3, after 52 months, that it was easily possible 

to maintain a polarization potential of —0.77 v 

with 0.1 ma/ft®. In fact, even with this low- 

current density the potential of specimen No. 3 

was usually more electronegative than —0.77 v. 
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During the last 3 months of exposure cathodic 
polarization curves were obtained on specimens 
Nos. 2 and 3. Curves were also obtained on speci- 
men No. 1 (freely corroding) before and after the 
removal of the heavy accumulation of corrosion 
products. After cleaning, curves were obtained 
during a 14-day exposure “period to provide data on 
the initial corrosion rate of the unprotected rod and 
for comparison with the curves obtained on speci- 
mens 2 and 3. 

Potential checks were made with a 200,000 
ohms/v voltmeter. Polarization data were obtained 
with a two-pen (current and potential) strip chart 
electronic recorder. The polarizing current was 
varied by linear increments in the applied voltage 
which were obtained with a synchronously driven 
10-turn potentiometer. The Holler null circuit 
described in a previous paper [7] was used in con- 
junction with the recorder but IR-drop compensation 


was unnecessary in the low-resistivity electrolyte. 
3. Electrical Measurements and 
Observations 
The more important information in this study 


pertains to the significance of the polarization data 
and the discussion later pertaining to the formation 
of the natural coatings and their importance in 
cathodic protection. Before discussing the polar- 
ization data, — ‘r, some of the prior ‘observations 
on specimens Nos. 1, 2, and 3 or in order. 

The potentials, wet and other observations 
pertaining to these specimens during the first 12 
months of exposure are of interest. Until the end 
of the first month specimens Nos. 2 and 3 each 
received 5.0 ma/ft? of protective current, with the 
potential of No. 2 around —0.97 v and No. 3 at 
-0.90 v. After the first month, the protective 
galvanic current on specimen No. 2 rapidly dimin- 
ished to 0.6 ma/ft? for at least 6 months of the 
vear and never exceeded 1.2 ma/ft?. On the other 
hand, the current applied to specimen No. 3 usually 
varied between 5 and 13 ma/ft?. The potentials 
for specimen No. 2 varied during the year between 

0.97 v and —1.0 v while the potential for specimen 
No. 3 varied between —0.80 v and —0.97 v, averag- 
ing about —0.85 v. All potentials were read while 
the anode was momentarily disconnected. During 
the year, the voltage applied between the carbon 
anode and specimen No. 3 was very close to 2 v, 
so it can be assumed safely that the driving voltage 
was greater than the driving voltage between the 
zine anode and specimen No. 2, even allowing for 
considerable anodic polarization of the carbon anode. 
At the end of the vear, the 3 specimens were removed 
for examination, cleaning and reweighing. Specimen 
No. 1 was cathodically cleaned in salt water at a 
current density of 1 amp/ft? for 3 hr and then wire- 
brushed by hand to the bare metal. The specimen 
lost 4185 mg during the first 12 months which is 
equivalent from Faraday’s Law (assuming 100% 
corrosion efficiency) to an average corrosion current 
density of about 4.6 ma/ft?. Specimens Nos. 2 and 
3 were washed under hot running water and examined. 
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Number 2 had a fairly adherent coating, similar to 
that formed on the preliminary demonstration speci- 
men exposed earlier. Although no rust was visible, 
the specimen was subjected to a cathodic current 
density of 1 amp/ft? for 1 hr. After considerable 
wire-brushing and reweighing, specimen No. 2 
weighed 110 mg more than before exposure. The 
increase in weight was attributable to the coating 
which was still visible after the wire-brushing. 
Similar treatment of specimen No. 3, also without 
evidence of corrosion, revealed an increase in weight 
of only 10 mg. Before cleaning, specimen No. 3 
had a dark-green appearance. Thus, during the 
first year, the coating formed on specimen No. 2 
was apparently superior to that on specimen No. 3. 
This explains the lower current density required to 
protect specimen No. 2. 

After the 3 specimens were cleaned and reweighed 
they were put back into their respective jars with 
no other changes. Currents and potentials were 
about the same during the next year as they had 
been previously. After 24 months of exposure, 
specimen No. 3 was taken out of the salt water for 
a few hours. It was brushed under hot running 
water with a brass-bristle brush, permitted to dry 
and then weighed. The coating now weighed about 
25 mg and its appearance was the same as after 12 
months. The specimen was put back into the salt 
water with cathodic protection. After 33 months of 
total exposure time the protective current, which 
was 3 ma/ft? at the time, was temporarily removed 
from specimen No. 3. Immediately after opening 
the circuit, the potential was around —1 v. After 
5 min, the potential was —0.82 v and after 35 min 
with no protective current the potential was —0.77 v. 
The protective current was then restored at an arbi- 
trarily reduced value of 1 ma/ft? which soon polar- 
ized the steel to —0.8 v. The next day the potential 
was —0.84 v and the current 0.5 ma/ft®, a condition 
that persisted for the next 7 days. The applied 
voltage between the carbon anode and the steel 
specimen was further reduced to 1.45 v and it was 
found that for several days 0.25 ma/ft? was sufficient 
to hold the potential at —0.77 v. During the next 
18 months the current density varied between 0.1 
and 0.5 ma/ft? except for one occasion when it went 
to 1 ma/ft? when the power failed over a weekend. 
This resulted in a pinhole through the natural coating 
as evidenced by a spot of rust. The coating soon 
re-formed over this spot and the current again grad- 
ually lowered until, between the 52d and 5: oth months 


of exposure it varied between 0.1 and 0.25 ma/ft? 
with the potential around —1.0 v. At c ‘se rela- 
tively low-current densities no attempt was then 


made to reduce the current further but it was 
obviously possible to do so, as will be shown later, 
in view of the potential. 

An experience with specimen No. 2 was also rather 
enlightening. After 42 months exposure, the zine 
anode was disconnected and the potential of the 
specimen measured at intervals during a period of 
8 hr. Prior to opening the circuit the protective 
current was 0.5 ma/ft®?. Immediately after breaking 
the anode circuit, the potential of the specimen was 








—0.96 v and after 30 min it was —0.83 v. After 
1 hr without protective current, the potential was 
still —0.78 v and after 64 hr it was —0.62 v with 
no visible evidence of rusting. The protective cur- 
rent was left off over the weekend. After 64 hr the 
potential was —0.665 v and two small spots of rust 
were evident on the coated surface, otherwise there 
was no evidence of corrosion. The rust was brushed 
off and the specimen put back into the solution 
under protection by the zinc. Over the weekend, 
although no measurements were made, the potential 
probably rose to a value somewhat more noble than 
—0.62 v and then when breakdown of the points on 
the surface occurred, the potential gradually drifted 
in the more electro-negative direction. A similar 
potential drift in the noble direction was noted on 
steel specimens in soils [8] after protective currents 
were removed and then a reversed trend in potential 
drift as corrosion again commenced. 

A repetition of the foregoing experiment followed 
by polarization measurements was conducted on 
specimen No. 2 after 56 months total exposure time. 
At this time, the galvanic current had decreased to 
0.1 ma/ft? from a value of 0.5 ma/ft® on the previous 
occasion. Within a few sec after opening the gal- 
vanic circuit the potential of the steel was —0.92 v 
at which time a cathodic polarization curve was 
obtained, The polarization was rather rapid, AV/A/ 
being 1000 mv/ma from 0 to 1 ma/ft?. After running 
the curve, the potential of the specimen which had 
been left on open circuit drifted to —0.77 v in 2'4 hr, 
and after 314 hr the potential was —0.75 v. Twenty- 
one hours later, the potential was —0.64 v and there 
was no visible evidence of corrosion. Cathodic po- 
larization curve A, figure 2, was obtained at this 
time. This was the first of a series of curves taken 
on this specimen during the following 36 days to 
show how damage of the coating affected the polari- 
zation, and the effect on polarization of removing 
the specimen from the salt water, esposing it to the 
indoor atmosphere for 12 days and then putting it 
back into the same water again. 

Curve B, figure 2, shows the polarization of speci- 
men No. 2 after 45 hr without protective current 
when the first visible evidence of rust appeared about 
2.5 in. below the waterline, or around 1 in. below 
the lower edge of the tape. Curve C, was obtained 
after 5 days without protective current and then 
the specimen was removed from the water in order 
to scrape off some of the coating about 3 in. below 
the lower edge of the tape. The coating was removed 
in the form of a 0.2 in. wide band encircling the rod. 
Curve D, obtained the next day, shows the effect 
of the coating removal on the polarization. Curve 
E, was obtained after 12 days without protective 
current, no additional corrosion being visible other 
than on the two areas just mentioned. On the 14th 
day without protective current, this specimen was 
removed from its corrosive medium and exposed to 
the air for 12 days with the other two specimens. 
Just prior to putting specimen No. 2 back into the 
salt water, a little dilute hydrochloric acid was placed 


on a small area of the coating. Gas evolution, indic- 
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Figure 2. Cathodic polarization of an inorganically coated 


steel specimen, No c. mas perce nt sodium chloride solution, 
as affecte d by progressively damaging the coating during a 
pe riod of 36 days without cathodic protection, 

Polarization rates A\V/AI, mv/m 


X0.1) indicated by the 
504: F, 320 and G, 353 


i, through the range of applied current (density 
. 7 


irrow ire as follows: A, 940; B, 808; (¢ ; D, 520; E, 


ative of a carbonate coating, was observed. The 
specimen was then rinsed under running water and 
returned to the salt water from which it had been 
removed. The next day, curve F, figure 2, was ob- 
tained. There was some visible ey idence of rust on 
the area touched by the acid the previous day. 
During the next 10 days there was no evidence of 
increased corrosion. On the 86th day without 
cathodic protection the final cathodic polarization 
curve, G, was obtained. After stabilization of the 
specimen potential, an anodic polarization curve was 
run. The specimen was then removed from the 
water for a spectrographic analysis of the natural 
coating formed during the 57 months of exposure 
and for reweighing of the specimen after removal of 
the coating. The anodic polarization was purposely 
left until the end of the experiment to avoid the 
possibility of disturbing the coating. The signifi- 
cance of the polarization curves will be discussed 
later. 

Similar polarization studies were made on speci- 
men No. 3 except that no attempt was made to 
intentionally disturb the natural coating. After 55 
months exposure, the protective current on specimen 
No. 3 was 0.1 ma/ft? and a quick reading of its poten- 
tial upon removing the current gave —1.05v. After 
the circuit had remained open for 5 min, the potential 
was still —0.98 v. A potential-time study was then 


carried out by use of the potentiometer-recorder. | ft eas 
After 35 min without protective current, the potential | 
was —0.97 v. Note the relatively slow decay in : - 
potential as compared with that observed in a similar mh 
experiment, previously described, after 33 months ex- o—e 
posure. After 24, 48, 72, and 96 hr, the potentials 8 TS 
were —0.89 v, —0.81 v, —0.81 v, and —0.76 v, 
respectively. The potential continued to drift in 
the electropositive direction until on the 9th day 
without protective current, the potential was —0.65 : oe 
v with no visible corrosion. Cathodic polarization 
curves obtained on specimen No. 3 during this initial 
9-day period without protective current showed that 
the polarization rate (designated AV/MS/) was rela- 
tively high in comparison with that observed on the 
corroding specimen No. 1. 7 eg 
The protective current from the rectifier was again os 
applied to specimen No. 3 and left on for 12 days. 
Within 4 hr, the potential was —0.84 v and the i 
current density 0.3 ma/ft?. After 7 days, the current din 
density was 0.03 ma/ft®, the potential —0.90 v. On 
the 10th day, the applied voltage from the rectifier ¢ e 
was reduced to 1.33 v between anode and cathode. 
On the 12th day, the potential of the specimen was 





0.85 v and the protective current density was 
0.02 ma/ft Cathodic polarization from this poten- 
tial revealed AV/A/ as being 3400 mv/ma between 0 : mae 
and 0.5 ma/ft? of applied current density. On the 
12th day the protective current was removed and the 
specimen left without protection. After 24, 48, and ganically coated steel specimen, No. 3, in a 3 percent sodium 
96 hr, the specimen potential was, 0.83 v, 0.30 v, chloride solution for 29 days without cathodic protection. 
and 0.76 v, re spectively. \fter the second 9-day Polarization rates, my/ma, are as follows: A, 2100; B, 2400; C ,2670; D, 2570; 1 
period without protection specimen No. 3 was re- | ~ — 
moved from the salt water and exposed to the air in 
the laboratory for 12 days. This was followed by 
11 additional days of exposure without protection to 
the same salt water from which it had been removed. 
Cathodic polarization curves were obtained during 
the two 9-day periods preceding exposure to the air 
and during the following 11-day period. The polari- 
zation runs were made at times when some corrosion a 
could have been taking place, that is, when the speci- =~ 6 
men potential was electropositive to 0.77 Vv as sane 
shown by a previous study [9]. These cathodic 
curves are shown in figure 3. Curve A, was obtained 
on the 4th day, curve B, on the 9th day of the initial 
9-day period, and curves C, D, and E during the _- 
second 9-day period preceding air exposure. Curves N 
F and G cover the final 11-day period of exposure to | 2 -.: NO. 3 
the salt water. On the final day of salt water ex- 
posure, an anodic polarization curve was also 
obtained on specimen No. 3, after which it was re- 
moved from the salt water, washed under running | | 
water and permitted to air dry. This was followed ‘ | 
by a spectrographic analysis of the coating, quanti- 
tative analysis for carbonates in the coating and 
finally, removal of the coating and reweighing. { 

The anodic polarization curves obtained on the 
final day that specimens Nos. 2 and 3 were exposed 
are shown in figure 4. The eventual superiority of 
the coating on specimen No. 3 as compared with that 
on specimen No. 2 is shown by the higher polarization 
rate of specimen No. 3. This will be discussed in on the final day of exposure corresponding with the time when 
more detail later in the paper. curves G (figs. 2 and 3) were obtained. — 
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FicgurRE 4. Anodic polarization of specimens Nos. 2 and 3 








In order to compare the polarization of specimens 
No. 2 and 3 having inorganic coatings with that of 
the unprotected specimen, No. 1, the curves shown 
in figure 5 were obtained. Curves A and B were 
run while the specimen was covered with the thick 
accumulation of rust, curve A before removing the 
rod for exposure to the atmosphere for 12 days, curve 
B after the rod had been back in the salt water for 
11 days following atmospheric exposure. Curves C 
through H were obtained on the same specimen 
during a 14-day exposure period after removing all 
the corrosion products and reexposing to the same 
solution. The specimen was weighed before and 
after the 14-day exposure period and anodic curves 
(fig. 6) were obtained to permit comparison of the 

calculated and actual we ight losses during the period. 


4. Interpretation of the Polarization Data 


The purpose of the discussion which follows is to 
show how the applied currents associated with 
changes-in-slope in polarization curves (usually refer- 
red to as breaks) are related to the actual polariza- 


tion over a specific range of the applied current 
(cathodic current here, because cathodic control 
predominated) and to show the relation between 
. . 2° e@ eee 
ee @€ €i6 hes. 
s al e e* @everve 


© + ce © © Peers eeee 
eee 
tee 
. 
+ 
%e 
+ ° ee ee Cees eeu, 
*e 
em, 
7 
o . © © eevee 
7 
fe. 
7 ° © © © weecs 








Figure 5. Cathodic polarization of a bare steel specimen, No. 


1, in a 3 percent sodium chloride solution. 


Curves A and B, after almost 4 years of exposure with 12 days of atmospheric 
exposure between polarization runs. Curves C through H, the same specimen 
solution for 
21; B, 16; C, 29; D, 40; E. 


after removing all corrosion products and exposing again to the same 
Polarization rates mv/ma are 
58; H, 59, 


14 days 
47; F, 


as follows: A, 
56; G, 
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Similar rela- 
tionships, with some variation, were previously 
shown for steel corroding in soils [10]. The signifi- 
cance of the slopes of polarization curves in relation 
to rates of corrosion has since been discussed by 
others, namely, Simmons [11], Skold and Larson [12], 
and Stern [13]. 

Curves A and B, figure 5 were the only cathodic 
polarization curves obtained on specimen No. 1 
(corroding freely) during the long exposure period, 
and while the specimen was covered with a thick 
adherent deposit of rust accumulated during almost 
4 yrs of continuous exposure. Figure 7 shows the 3 
specimens after 56 months of exposure to the salt 
water. The corrosion products on No. 1 are evident 
in the photograph. Reweighing of specimen No. 
after the final cleaning revealed a weight loss of 
11,835 mg during the intervening 1 days of 
exposure beginning after 12 months. Using the 
electrochemical equivalent (K=2.8938 x 10~* g/cou- 
lomb) and assuming a corrosion e flicienc y of 100 per- 
cent, the average corrosion current density was 
calculated to be 3.5 ma/ft®. It will be recalled that 
the average corrosion current density was 4.6 ma/ft? 
for the first 12 months. <A dropping off of the 
corrosion rate with the poeeseet of time is indicated. 
For this reason curves C through H, figure 5, and 
D through G, figure 6, were obtained to provide 
data on the corrosion rate of the clean steel rod during 
the first 14 days of exposure. The anodic curves, 
figure 6, are plotted on rectangular coordinate paper 
in order to permit a more accurate estimation of the 
current J,, than was possible on logarithmic coordi- 
nates. The data shown in table 1 were obtained 
from figures 5 and 6 and are used to calculate [10] 
the weight loss over the 14-day period of exposure. 
The corrosion current with time, being 
down to 6.8 ma/ft® after 10 days of exposure. The 
weight loss (299 mg) as calculated from the cathodic 
and anodic polarization curves compares reasonably 
well with the actual weight loss (280 mg). 


these data and the rates of corrosion. 
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Figure 6. 1 during the 


14 day exposure period mentioned in figure 6. 
¢ I Ju 


Anodic polarization of Spe cinmen No. 


Curves D, E, F, and Were obtained on the same day, 


D, E, F, and G (fig.). 


, 


respectively, as curves 











FiGureE 7. 


| ppea ance of steel spec~mens Nos. 1, 2, and 8 
before final cleaning 
Rust spots account for most of the we t loss on specimen No. 2. The dark 
appearance below the tape on No. 3 is not rust but the color of the inorganic 
ceating 
TABLE 1. Weight loss on specimen Vo. 1 caused by corrosion 
in 14 da Ss 
Calculations from polarization data (figs. 5 and 6 
Polar g current * Weight loss 
it break in curve Corrosion 
Exposure time current ame 
Cathod Anodic, Calculated*®| Actual 
I Iq cumulative 
Days 1a ma va mq mg 
l 1. 95 ¢ 19d 1.4 35 ical 
2 1.55 D 3.9 D 11 66 és og 
4 1.35 E 3.0 E 0. 93 116 ie 
7 1. 20 F 2.4 F . 80 180 . one 
10 1.05 G4 2 0G 68 234 sovessecnne 
l4 1.00 H 1.9 65 299 280 


* Measured values (density 0.1 Che capital letter following value of current 
designates the curve in fig. 5 or fig. 6 from which the value is obtained. 

» te p 14/1; 

© Weight loss Kti, where K=2.8938X10- g per coulomb, i=i,=average 
current (amperes) for the period (f in seconds) between successive readings. For 
the period before the initial polarization run, the initial value of i, was taken as 
average 

1 Based on the ratio 79/7, of the 2d day 

¢ Based on the ratio J,/J, of the 10th day, 


The high values of J, (curves A and B, fig. 5) are 
not consistent with the corrosion rate, which pre- 
sumably was at its lowest value when the data for 
these curves were obtained. Of course, there is no 
way of knowing what the corrosion current actually 
was at the time, except that it should be about 
equal to the average value (3.5 ma/ft®) previously 
calculated, based on the assumption that the cor- 
rosion rate was relatively stable for a long time. 
The answer does not lie in a change to anodic control, 
because an anodic polarization curve (not shown), 
obtained after curve B, showed that cathodic control 
predominated. It appears as though more work 
ought to be done to determine the effects of a heavy 
accumulation of corrosion products on polarization 
of the basis metal. 

After 56 mo exposure and before finally complet- 
ing the experiment it was decided, as previously 
mentioned, to leave specimen No. 2 without pro- 
tective current, observe the effect of scratches in the 
coating on polarization and possibly measure the 
actual corrosion which had occurred. The polariza- 
tion data for specimen No. 2 are shown in figures 
2 and 4 and the corrosion calculations based on these 
data in table 2. The only other time in almost 5 
yr of exposure that this specimen was subject to 
any corrosion was for a portion of the 3-day period 
after 42 mo of exposure. Based on subsequent 
measurements (table 2), the weight because 
of corrosion during this period was probably less than 
5mg. The weight loss calculations are based mainly 
on the cathodic polarization currents, J,, as the 
anodic curve (fig. 4) showed that the control was 
cathodic. Use of the ratio 7,=0.73 J, seems reason- 
able, in view of the data (table 1) obtained on 
specimen No. 1. Thus, for all practical purposes, in 
view of the agreement between actual and calculated 


loss 


weight losses, increases in J, were indicative of 
increases in corrosion rate. This was reflected in 


») 


curve D (fig. 2) as a result of a scratch in the coating 
purposely made on the preceding day. There was 
also some indication that exposure to the atmosphere 
for 12 days or the disturbance of the coating by a 
drop of dilute acid just prior to immersion increased 
the corrosion current (curve F). 

The cathodic polarization curves (fig. 3) for speci- 
men No. 3 cover three periods when there was no 
protective current. Although no corrosion was visi- 
ble at any time during these periods, and no attempt 
was made to scratch the coating, the potentials at 
zero current (substantially the same as those shown 
for the lowest current on the abscissa) are electro- 
positive to the protective potential —0.77 v and hence 
it was assumed that the possibility of corrosion, 
however slight, did exist. This does not imply that 
the potentials per se gave proof of corrosion, but 
that the potentials (fig. 3) were in the same range as 
were the potentials of specimen No. 2 (fig. 2), when 
corrosion was known to be taking place. Thus, No. 
3 was in the salt water for 29 days without protective 
current and for 21 days it is assumed, on the basis of 
the potential, that some corrosion could have oc- 
curred. During the almost 5 yr that this specimen 
was under cathodic protection, there was only one 
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TABLE 2. Weight loss on specimen No. 2 resulting from corrosion 
due to scratches in the coating 
Corrosion is calculated from the polarization data (figs. 2 and 4 


Polarizing current * Weight loss 





Exposure it break in curve Corrosion ¢ 
time current ? 
Calculated Actua 
Cathodie J, Anodic I, cumulative 
Days ma ma ma n I 
1 0. O72A 0. 053 1.3 
2 Om4B O61 7 
5 o9RCc O71 7% 
¢ 17D 12 10.0 
12 165 12 28 0) 
14 i 34.0 
Che specimen was removed from the salt water on the 14th day and put back 
iain on the 26th day Corrosion during this interim (12 days) was cons 
as negligible 
27 19} 14 37 
37 20G 0. 55 15 7H. ( 
® Measure V l rt pit { 
current designat ! f 2) fre which tt 1 btained 
in=Ip I n J,/1» ratio for the 37tl 1.731 
Same as footne 
Used 0.16 ma for J 
e From fig. 4 
occasion, over a weekend, when the current acci- 


dentally went off, as described in the previous section. 
The only anodic polarization curve on specimen No. 
3 was obtained on the last day of exposure and is 
shown in figure 4. Based on this curve and curve G 


(fig. 3), the ratio 1,/l,=approx. 3, and therefore, 
i,=0.75 1,. Using the value /,=0.03 ma (specimen 
area 0.1 ft?) from figure 3 as the average value 


for 21 days and calculating as previously, the caleu- 
lated weight loss was about 12 mg. The actual 
weight loss during the 5 vr was 35 mg. The only 
visual evidence of corrosion on the cleaned specimen 
was on the lower-half of the taped area and also 
above the tape, where salt had crystallized and 
penetrated the oil film. As the calculated corrosion 
would not include the area above the tape which was 
out of the solution, the calculated value seems to be 
quite reasonable. 

It is known from experience that there are times 
when it is difficult to observe significant breaks in 
polarization curves. The relative importance of 
cathodic and anodic polarization curves in estimating 
rates of corrosion was discussed several vears ago 
[10]. It was observed that, for steel in soils, a rela- 
tion existed between the rate of polarization (AV/A/ 
resulting from externally applied current) and the 
applied current J, (at the break in curve). It was 
also shown that both of these factors were related 
to the instantaneous rates of corrosion. The most 
useful values of AV seemed to be those potential 
changes which occurred in a range of applied currents 
of the order of J, or slightly larger. 

Simmons [11] found that inhibitors could be evalu- 
ated by AE/AT (called AV/A/ in this paper) from 
polarization curves. Small current densities caused 
changes in potential on exposed steel of 100 mv and 
more. In figures 2 and 3 of Simmon’s paper are 
shown cathodic curves for steel in salt water with 
and without the inhibitor Na, CrO,;. We noted that 
in both of Simmon’s figures the changes in potential 


the values at the breaks (corresponding to our /, 
values), are about equal. As these currents were of 
the relative magnitude of 10 to 1, Simmon’s con- 
ceeded this to be the relative magnitude of the 
corrosion currents which were either comparable in 
magnitude or less than the currents at the breaks 
in his curves. Thus, the ratios A//A/ for these data 
were in the reverse order, namely 1 to 10. 

Skold and Larson [12] showed that an empirical 
relation existed between the slopes (AE/Al) of 
polarization curves for low values of applied current 
and the corrosion of steel in fresh water 
determined by weight losses. They first tried cor- 
relating the breaks in polarization curves with cor- 
rosion rates but abandoned this in favor of the slope 
method. It is believed that the use of a null circuit 
is necessary for the detection of breaks in polarization 
curves obtained with a 3000 ohm-em solution such 
as was used by Skold and Larson. Nevertheless, 
their straight line relationship on logarithmic CO- 
ordinates between corrosion rate and A//A/ seems 
to have considerable merit. 


rates aus 


This has been further 
substantiated by the data reported in this paper. 
Stern [13] has published the theoretical basis for 
the method of Skold and Larson. 
the importance ol 


He has stressed 
applying polarizing currents which 
are smaller than the presumed corrosion currents in 
order to avoid marked polarization with resultant 
possible changes in the inherent rates of corrosion. 
Apparently, the inherent rates of corrosion observed 
in the present study were unaffected by the applied 
currents which at the end of a run were often 4 
the magnitude the corrosion currents. 
This is also true ol previous studies [7, 9] where 
controls subjected to reasonable amounts of polariza- 


times ot 


tion lost no more weight than those on which no 
polarization runs were made. Possibly the short 
polarizing time (10 to 15 min) is’ beneficial in 
negating this tendency. 

The polarization data (figs. 2 to 6 inel, were 


plotted from the strip chart records for additional 
interpretation and condensation. Furthermore, on 
the logarithmic current scale, the straight line portion 
of the eathodic curve, for applied eurrents larger 
than J,, becomes helpful in estimating AV/A/. The 
values of AV/A/ shown in figures 2, 3, and 5 were 
calculated from changes in potential between zero 
applied current and the currents indicated by the 
arrows (beginning of straight line) where all local 
action corrosion 1s presumed to have ceased. Figure 
8 shows the relation between AV/A/ and correspond- 
ing values of J, for the 3 specimens. It is interesting 
to note that the plot of these parameters on the 
logarithmic coordinates results in a slope close to the 
theoretical value of minus one mentioned by Stern 
[13], although he considered corrosion current instead 
of the current /,. Thus, figure 8 suggests, at least 
for steel in aqueous solutions, that a fairly constant 
ratio exists between the corrosion current and the 
current J,. The calculations of weight loss in tables 
1 and 2 are to some extent based on this premise. 
The small variations in the ratio 7,/7, may 


be at- 


of the d-c null curves, from the lowest currents to | tributed to errors in selecting the break current J, 
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POLARIZATION RATE, AV/AI, mv/(ma/ft? 

FiGguRE 8 Re lationship during cathodic polarization of steel 
specimens Nos. 1, 2, and 3 between the current, T,, at the 
break and the polarization AV/AI, on logarithmic 


coordinate Ss. 


Se. Specimen No 


rate, 


1 (fig. 5):@, No. 2 (fig. 2); ©, No. 3 (fig. 3). 


from the anodic polarization curves. The selection 
of J, is more subject to error than /, when the 
corrosion reaction is under cathodic control It is 
fortunate though that the accuracy with which the 
corrosion current can be caleulated depends chiefly 
on the more significant curve which also lends itself 
to the better interpretation. In a case of corrosion 
under anodic control, this would be the anodic curve. 

In order to further compare the data (fig. 8) with 
that of Skold and Larson, the actual weight losses 
of specimens 1 and 2 (tables 1 and 2) and the caleu- 
lated weight loss of specimen 3 were related to the 
polarization rates AV/A/ in their units of measure- 
ment (fig. 9 The ve weight losses, even though 
for only 3 specimens, carry considerable significance 
because of their range. The values are time-aver: iged 
and cover time periods of 14, 25, and 21 days, 
respectively, for specimens Nos. 1, 2, and 3. The 
calculated weight loss was used for No. 3 because, as 
previously stated, the actual value (35 mg) consisted 
chie fly of corrosion loss which occurred outside of 
the electrolyte and therefor could not be accounted 
for by the polarization measurements. Curves B 
and C were drawn from the data (fig. 1) of Skold and 
Larson [12]. Curve B covers approximately their 
range of data in the 3000 ohm-cem water and curve 
(is based on the mean value of each of two groups 
of their data. Curves B and C are subject to some 
error, B more than C, as mentioned by 
Skold and Larson and also by Stern [13] and therefore 
have greater negative slopes than has the NBS 
curve, A. 
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POLARIZATION RATE, AV/AI, v/(ma/dm? 
Figure 9. Logarithmic relationship between actual corrosion 
(weight loss) and cathodic polarization rate of steel. 
Curve A (NBS), @ Specimen No. 1, @ No. 2, No 
solution. Curves B and C [Skoid and Larson] in 
350-850 ohm-cm solution, respectively. 


3 all in about 35 ohm-cm 
3000 ohm-cm water and 


5. Nature of the Coatings 


After 57 months of exposure, approximately 50 
mg of coating was taken from specimens Nos. 2 and 3 
for analysis. The results of the spectrochemic ‘ 
analyses are shown in table 3. The coating of 
specimen No. 2 consisted mainly of carbonates 
(27% CO,) and silicates associated with calcium, 
magnesium, iron, and zinc. As the coating was 
almost white in appearance, the iron content was 
probably at the lower limit of the indicated range. 
The coating on specimen No. 3 was primarily in 
the form of carbonates (56% CO ;) combined with 
decreasing amounts of calcium, magnesium, and 
iron, in that order. Some silicate was also present. 
The zine and silicates apparently constitute the chief 
differences in coating composition of the two speci- 
mens. The coating of specimen No. 2 had both of 
these elements in considerable proportion, whereas 
the coating on specimen No. 3 did not, it belie 
chiefly calcium carbonate with minor amounts of 
magnesium carbonate and iron oxide. It was 
estimated that during the 57 months about 300 liters 


of city water, which accounted for the formation of 
the coatings, were poured into rach specimen jar 
to replace the water lost by evaporation. A com- 


posite mineral analysis of the city water (District 
of Columbia) for the year 1958 (avg) is shown in 


table 4. 
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The saline solution surrounding specimen No. 2 
probably had a zine ion concentration corresponding 
to the solubility product of zine carbonate or hy- 
droxide. _Carbonates were evident as a sediment in 
the container. A sediment of carbonates was also 
very noticeable in the jar containing specimen No. 3. 

The alkaline reactions that occur on the surface 
of steel during cathodic protection result in precipi- 
tates which either adhere to the metal surface as a 
coating or drop as a sediment. Some of the possible 


reactions are: 
Ca**+ +HCO;+OH-—CaCO,;+ H,O 
Me** + HCO; +OH---MeCO,+ HO 
Fe++ +HCO; +OH-->FeCO,+ HO 
Zn** +HCO;+0OH anCO, 4 H.O 
2 Fe***+6(OH -Fe,0;,+-3 HO 
Mg** +2 (OH)~—Mg(OH), 
Zn** 1/n(SiO,H,)n + 2(OH)~--ZnsSiO,+ 2 H.O 
Zn** +2(OH)~—Zn(OH), 


As a result of the foregoing reactions, the compounds 
formed can recombine and form other molecular 
compounds. The proportion of any one compound 
in the composition of the precipitate depends upon 
its solubility and the constituents of the solution. 
Based on table 3, the coating on specimen No. 2 
was, in addition to calcium carbonate, probably 
zine silicate. It is believed that the large amount of 
silicates in the coating of specimen No. 2 can be 
attributed to the zine in solution from the anode. 
This conclusion is reinforced by the general knowl- 
edge that zinc is effective in precipitating silicates 
out of solution. The relative absence of silicates 
in the coating on specimen No. 3 can only be 
attributed to the relative absence of zinc, as the 
solutions were similar, but a carbon rather than 
zine anode was used for ‘iteniie protection. 


TaBLeE 3. Results of spectrochemical analyses of coatings 
formed on specime ns nos. 2 and 3 
Specimen Specimen 
Element Element 
No. 2 No. 3 No, 2 No. 3 
Ag Vw I Li 1 1 
Al M W Mg Ss 8S 
B T Mn M W 
Ba W W Na M M 
Ca Vs Vs Ni W VW 
Cd W Pb M WwW 
Co Vw Pt Vw 
Cr W Vw Si VS M 
Cu M W Si M W 
Fe Ss Ss Ti VW 
K Vw Vw Zn Ss M 
Norte: In general, VS means greater than 10%; 8, 1 to 10%; M, 0.1 to 1.0%; W, 


0.01 to 0.1%; V W, 0.001 to 0.01%; T, 0.0001 to 0.001%: , hot detected. 


TABLE 4. 
Total solids Loss SiO Ca Cl Co I Fe 
ign 
15 37 4 31.7 10.0 2.1 1. 00 026 





Mineral analysis of District of Columbia water supply 
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The 


be a 
contributing factor in the formation of the coatings. 


mechanism of electrophoresis may 
Such a hypothesis would necessarily involve the 
attraction of positively-charged colloidal particles 
to the cathode. Although a highly conductive 
solution, such as 3 percent NaCl, is not a favorable 
environment for electrophoresis to take place, it 
is believed that over a long period of time the con- 
tribution of the mechanism to the formation of the 
coating, although perhaps taking place slowly, might 


be a significant factor. Thus, the formation of 
coatings on specimens Nos. 2 and 8 may then be 
attributed to two factors, first, the precipitation 
of relatively insoluble compounds, such as ear- 


bonates, silicates, and others in physical contact 
with the specimen surfaces and, second, the electro- 
phoretic deposition of colloidal particles in the pores 
of the coating. 

A cathodically protected surface promotes condi- 
tions favorable for the precipitation of carbonates, 
phosphates, hydroxides, and silicates. If the pro- 
tective current is excessive as the coating forms, the 
alkalinity around the uncoated areas will promote 
additional precipitation farther from the cathode 
surface because of the diffusion of hydroxide away 
from the cathode. On the other hand, a reduction 
of the protective current, consistent with the pro- 
tective potential requirement, would tend lessen 
this diffusion. For the formation of a compact 
coating, the volume and mass of insoluble colloidal 
particles might also be considered. It is believed 
that both volume and mass of the colloids increase 
with time during the precipitation process. Exces- 
sive alkalinity around the cathode areas, for example 
a result of too much current, will promote this 
increase because the precipitation begins further 
away from the cathodic areas. This will result: in 
two unfavorable effects, first, a slowing down of the 
migration of colloids toward the cathode and, second, 
the inability of the larger colloids which do arrive 
at the cathode to fill in the pores of the coating 
already in the formative stage. On the other hand, 
lower alkalinity in the porous or uncoated areas, 
as a result of reduced current, will result in smaller 
colloidal particles arriving at the surface, better 
penetration of the pores and hence a more compact 
coating. This condition is approached by reducing 
the protective current to the optimum value (con- 
sistent with protective potential) for cathodic 
protec ‘tion, 

The foregoing discussion might be evaluated by 
considering the properties of the inorganic coatings 
which were actually formed on specimens Nos. 2 and 
3. During the first 24 months of exposure, the coat- 
ing on No. 2 apparently was superior to that on No. 
3, as evidenced by the heavier coating on the former. 
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Hco Me Mn NO NO SO4 Alk Tot. hard C.H 
CaCO CaCO CaCO 
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It was also observed that the protective current 
density on No. 3 was several times greater than on 
No. 2. Current densities were about the same the 


following year and after 24 months the weight of 


coating formed on No. 2 was at least 4 times that on 
No. 3 (25mg). However, during the period between 
the 33d and 57th month of exposure the protective 
current to specimen No. 3 was progressively lowered 
qv reducing the applied voltage to that value which 
brought the potential of the specimen to —0.77 v. 
AKentually the current to No. 3 was smaller than that 
to No. 2, and the polarization of No. 3 was accom- 
plished at a smaller current density than the polari- 
zation of No. 2 (fig. 8). 

When the specimens were being cleaned for final 
reweighing, the coating on No. 3 was observed to be 
considerably more adherent than that formed on 
No. 2, even after subjecting the specimens to cath- 
odic cleaning (1 amp/ft®) for 2 hr. The coating on 
No. 2 was scraped off rather easily with a brass 
scraper but that formed on No. 3 had _ to first be 
hammered tapping. Unfortunately, no 
accurate comparison of coating weights can be given, 
but based on the weight of a portion (500 mg) of the 
coating removed from No. 3, it was conservatively 
estimated that the coating formed on No. 3 weighed 
about 1000 mg and its thickness was about 0.010 in. 
Thus, the thickness, adherence, and protective value 
of the coating apparently increased after lowering 
the protective current to the optimum value as 
determined by the potential of the specimen. 
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This paper discusses the errors from residuals in the Maxwell-Wien bridge and the 
effect of these on the measurement of inductors in a bridge not having a Wagner ground. 
Particular attention is given to the use of substitution methods for accurate measurements 
and especially to the ‘“‘equal-substitution’’ (comparison) method, which can yield excellent 
precision in the calibration of inductance standards 


1. Introduction 


Self inductors and mutual inductors can be con- 
structed to have inductance that is computable from 
their geometry and dimensions. Although the in- 
ductance of certain inductors having measurable 
mechanical dimensions and simple geometrical form 
can be computed with excellent accuracy, such 
inductors are not often used as reference standards 
for inductance measurements. In practice, induct- 
ance Measurements are more conveniently made with 
reference to noncomputable reactances in the form 
of compact and stable reference standards of induct- 
ance or capacitance. Accurate values of inductance 
can be determined by comparing the reactance of 
the inductor with the reactance of either standard 
capacitors or standard inductors by use of appro- 
priate bridge circuits. However, the most precise 
measurement of an unknown quantity (in this case 
inductance) is made by comparing it with a standard 
of like kind and magnitude, the small difference being 
measured by a corresponding small change in one 
element of the bridge circuit. When comparison 
methods are used, it is often possible to reduce the 
detrimental effects of residuals in the measuring 
circuit to such an extent that the precision of com- 
parison is much better than the accuracy with which 
the calibration standard is known. 

The Maxwell-Wien bridge circuit has long been 
used for the accurate measurement of inductance. 
For measurements of best accuracy with any alter- 
nating-current bridge circuit, even for comparison 
measurements, it is necessary to consider the effects 
of residuals, the ways of reducing these effects, and 
the handling of corrections to offset the net errors. 


2. Maxwell-Wien Bridge 


A bridge circuit originally developed by J. C. 
Maxwell [1]! for use with ballistic detectors was 
adapted by M. Wien [2] for a-c measurements. For 
many years the circuit has been used widely by 


*Contribution from the Radio Standards Laboratory, National Bureau of 
Standards, Boulder, Colo 
Figures in brackets indicate the literature references at the end of this paper. 


standardizing laboratories for measurements of in- 
ductance. Figure 1 shows a schematic of a well- 
shielded Maxwell-Wien bridge with Wagner arms for 
elimination of errors resulting from current diverted 
to ground through leakage impedances. The in- 
ductance, L, to be measured is placed in series with 
| a resistance, r,. If the components of the bridge 
arms shown in figure 1 were pure, the inductance 
measured when the bridge is balanced would be 


L ORpRs, (1) 


and to obtain the balance it is necessary that 


— 
an) 
— 
-~ 
~ 


PR ek (2) 









































on 











The Maxwell-Wien bridge circuit with complete 


| . 
| FicureE 1. 
| shielding and a Wagner ground. 
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In practical situations it is not possible to obtain 
completely pure components. The residuals in the 
components contribute to the systematic errors, and 
it is necessary to analyze the circuit sufficiently to 
correct for these errors or to determine the maximum 
uncertainty in the results if the errors are neglected. 
Figure 2 shows the components of figure 1 with the 
addition of appropriate residual components, but 
without the Wagner arms. For example, the series 
resistor, r, of figure 1, in the inductance arm is re pre- 

sented in figure 2 2 bya pure resistance, 7,, in series 
with an inductance, l,. The inductanc e, lr, is a 
lumped constant that may be either positive or nega- 
tive depending upon the relative magnitudes of 
inductive and capacitive effects. Although lumped 
constant residuals are not independent of frequency, 
the lumped constant concept provides a_ useful 
equivalent circuit at low frequencies. The other 
resistive components are similarly represented by a 
pure resistance and an effective series inductance. 
The variable capacitance, C, is presumed to contain 
an effective series resistance, r. The resistance 
R, in figure 2 is the equivalent series resistance of the 
inductor. 


3. Equation of Balance 
The impedances of the four bridge arms are 


. [(R+ ql) (1+-@?r?C?) 
Zane = “ - a one S —) (3) 
: + wC(wrC- 7) (R+ jal) 





a ae (4) 
Z; Rr ren qw(L T l,), (5) 

and 
ZB 4-deile. (6) 


The derivation of eq (3) is given in the appendix. 
When the bridge is balanced 


LacZ, Z pls. (7) 


After the multiplication indicated by eq (7) is per- 
formed, the equation of imaginary components can 
be written 

7 rzl = Rslp- 


(14+-orC?)[R(L+1,) + (Rk -Rpls] 


= (RRpRs4 T -w(RR pls +RRlp+RpRsl) orC 


—w?(Rl pls : Rplls - Rollp) — wll plswrC}. (S) 

Considerable simplification of eq (8) is possible if 
the residuals are expressed as time constants and 
dissipation factors, as follows. Let 


wrC=De, (9) 
l 
Rk 7*® (10) 
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Figure 2. The Mazxwell-Wien bridge circuit showing the resid- 
uals in the arms as lumped constants. 


lp . 

la a (11) 
and 

le 

te (12) 


With negligible error it will be found that 


L=CR pRs{1—w(1 prs +t et p+rets) —D, 
tr) De—w*t ptst RDe| 


> I RRs ‘ ‘ 
-(R, T 'y)/T rT R (T pt Ts) -l;, (13) 
t 


t+w(tpt+rs4t 


and a convenient simplification results if the small 
ve in brackets are replaced by the symbol K. 
Then the equation of balance becomes 


=a R a , 
L ( RpRs(1-4 K)- (PR, T ry) tr+ ct (Tp T Ts)—lz, 
(14) 
where 
K=—@(t pts ttetpttrts) —D, 
tw(tet+tstTr)De—@* 1 ptsteDe. (15) 


4. Effect of Stray Capacitance 


Before analyzing this equation in detail it is profit- 
able to consider the effect of stray capacitance across 
the terminals of the inductor. Such capacitance 

can be introduced by the leads connecting the 
inductor to the bridge (see fig. 3). For this analy sis, 
whatever capacitance exists within the case of the 
inductor is considered not a part of Cs, but rather 
part of the inductor. It should be understood that 
this internal shunt capacitance contributes to the 


wv gf 
Ps lr / L 
BRIDGE Cs= BRIDGE 
TERMINALS Rr TERMINALS Me 
“ea . 
(a) (b) 


Figure 3. The inductor under test, Ly, is connected to bridge 
terminals with leads hat ing stray capacitance, Cs, between 
them 


The bridge measures the equivalent inductance, L, which differs from Lr 
because of the stray capacitance 


effective series inductance measured at the terminals 
of the inductor, thereby causing the effective series 
inductance to vary with frequency to some extent. 

Figure 3a shows the equivalent series inductance 
under test, Lr, in series with the internal resistance 
of the inductor, Rr. The stray capacitance, Cs, 
external to the inductor is in parallel with Ly and 
Ry. The bridge measurement does not completely 
separate these elements but indicates the effective 
inductance, 1, as shown in figure 3b. The effective 
series inductance, Ly, at the terminals of the inductor 
is related to the measured value, 1, at the ends of 
the leads by the approximation 


L, L( l an S 4? R2( 's)° (16) 
T / 


The derivation of eq (16) is given in the appendix. 
The last three terms in parentheses constitute the 
small correction terms resulting from stray capaci- 
tance. The first and the last of these vary as the 
square of the frequency. The last term is usually 
much the smaller of the two in practical situations. 
The second correction term is not dependent upon 
frequency and must never be overlooked as a source 
of error in measurements of the highest accuracy. 
For example, in the measurement of a 10 h inductor 
having an internal resistance of 10,000 ohms, at the 
end of a coaxial cable having a capacitance of 100 pf, 
the error from neglecting the second correction term 
is 0.01 h, equivale nt to 0.1 percent of the quantity 
measured. Often the first correction term is sig- 
nificant while the other terms are negligible, and the 
effective series inductance then can be written 
L,;~=L(1—eL 7Cs) (17) 
However, in this paper both the first and second 
correction terms will be carried in the analysis to 


follow. 


5. Substitution Methods 


The existence of the residuals and the consequent 
correction terms in eqs (14) and (16) are detrimental 
to good accuracy for residuals are often very difficult 
to measure or estimate to better than an order of 
magnitude. To eliminate the effect on measure- 
ments of some of these residuals, substitution 














methods are used. For example, if a standard 
inductor of equal nominal value is substituted for the 
test inductor and the difference of inductance is 
measured by changing C and r, (or F) slightly, then 
the residuals in the resistance arms of the bridge have 
negligible effect on the measurement. Also, the 
effect. of the stray capacitance across the inductor 
under test or across the standard inductor is much 
reduced. The standard inductor must be very 
accurately known, however. 

If accurate measurements of large inductors are 
contemplated and a shorting plug ‘of negligible, or 
relatively small, inductance is substituted for the 
test inductor, Cs; must be reduced to a negligible 
value by complete shielding, and a Wagner ground 
should be used. The effect of the residual, /,, in the 
resistor, 7, is eliminated by this method only if /; is 
the same for the two conditions of balance. 

The latter method, in which a relatively 
inductance replaces the inductor under test, 
designated the ‘“zero-substitution” method to dis- 
tinguish it from the ‘“equal-substitution” method, 
in which the test inductor is replaced by a standard 
of equal nominal inductance. 

A general formula for inductance measurements by 
substitution methods with the Maxwell-Wien bridge 
may be derived from the two balance equations, 
the first with the inductor, Z7, in the circuit, and the 
second with a standard, Ly, in its place. If the 
small correction terms represented by A in eq (14) 
are constant for both conditions of balance and if 
(sis changed only slightly (to Csv) when the standard 
inductor or shorting plug is connected, then with 
certain justifiable approximations 


small 
can be 


Lin Lon + (O—Cy)RpRS{I Ln)Cg : Kk] 


—w*(Lr t 
— wo Lrby(C g —( SN ) 


(R2.— R%)Cs—(L,—lin). (18) 
The derivation of eq (18) is given in the appendix. 
This equation is useful in showing the effect of un- 
certainties in the magnitude of circuit components 
on the determination of inductance. 

If good quality components are used and the bridge 
is well designed, the time constants that contribute 
to K will be small. However, K cannot be deter- 
mined with certainty and is somewhat variable 
(dependent upon actual circuit constants). The 
uncertainty in K sets an ultimate practical limit to 
the accuracy of measurement of inductance by the 


Maxwell-Wien bridge. 


6. Zero-Substitution Method 


When the zero-substitution method is employed 
the small inductance, Ly, need not be accurately 
known, but the capacitor, C, and the resistance 
product, RpRs, mene be known with better accuracy 
than that desired for the inductance, Ly. It is 
possible to construct very small inductances that are 
computable from geometry and dimensions. A short- 
ing bar between two terminals has a finite inductance 
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that can often be calculated or estimated well enough 
to serve as the “‘zero-reference”’, Ly, for bridge meas- 
urements of inductors of much larger magnitude. If 
the measured inductance is defined as the increase 
in inductance when a shorting-link or switch associ- 
ated with the inductor is manipulated in a specified 
manner, then the actual inductance of the short 
circuit, Ly, is immaterial. The term (/,—J/,x) ean 
be minimized if the resistance, rz, is inductance com- 
pensated; i.e., if designed so that the effective series 
inductance remains constant as the resistance 
varied. The uncertainties in this term may limit 
the accuracy of measurement of small inductors. If 
Cs is not extremely small, it is evident that the 
uncertainty in Cs can have a significant effect on 
the accuracy of measurement by this method. In 
eq (18) the term in brackets containing Cs can 
become objectionably large when large values of 
inductance are being measured. For example, if the 
inductance being measured is approximately 10 h, 
if Cs=100 pf+10 pf, and if o=10,000 radians per 
second, then wLC,=0.10+0.01. Thus the correction 
term is 10 percent of the measured inductance, and 
the measurement uncertainty is 1 percent from this 
cause alone, disregarding the probable lack of validity 
of the assumptions (when the correction terms are 
so large) that were used to derive eq (18). The effect 
of capacitance Cs thus places a definite limitation 
on the accuracy which can be obtained by the zero- 
substitution method with a grounded bridge. In 
order to obtain better accuracy when «ZL is large 
it is necessary to use the “‘equal-substitution method”’ 
(described below) or to resort to separate shielding 
on the inductor leads with an ungrounded bridge 
using Wagner arms. This latter possibility involves 
considerations beyond the scope of this paper. 


is 


7. Equal-Substitution Method 


If an accurately known standard inductance, Ly, 
is substituted for the test inductor, Ly, the capaci- 
tance, C, and the resistance product, pls, need not 
be so accurately calibrated. Between the two con- 
ditions of balance the only circuit components 
changed are C and r,, both by small amounts. The 
residuals De and /, are essentially unchanged in this 
vase, and all the other residuals are identical for both 
conditions of balance. For the equal-substitution 
method the general eq (18) can be rewritten with 
negligible additional error as 


Le T (C’- Cy) RpRs( l 


wLlF(Cs— (hz 


2w?LeCs+K) 


feng) RR)Cs—(l,—lin). (19) 

In eqs (18) and (19) the stray capacitances, (s 
and Csy, are not considered to be exactly equal; 
however, it is expected that they will be maintained 
nearly equal. The importance of this, even in equal- 
substitution measurements, is evident from eq (19). 
In the comparison of 10 h inductors at an angular 
frequency of 10,000 radians per second, if the differ- 
ence Cs—Cgy is only 1 pf, then wL}(Cs—Csgy) is 
0.01 h. This error of 0.1 percent shows the impor- 





tance of maintaining (s nearly constant during 
the substitution. 

The term (/,—J/,y) can be reduced to a minimum 
by employing an adjustable resistor, r,, having in- 
ductance compensation. Uncertainties in this term 
may limit accuracy if the resistances of the test and 
standard inductors differ greatly, and the inductance 
compensation of the resistor, 7,, is inadequate. The 
adjustable resistor is in series with the inductor and 
is connected to the grounded corner of the bridge 
as shown in figure 2. If, instead, this resistor were 
connected to the ungrounded end of the inductance 
arm of the bridge, the capacitance to ground from 
the several decades of this resistor would shunt the 
inductor, and as the resistance is varied to accom- 
modate inductors having different internal resistance 
the change in stray capacitance would cause an error 
that would be difficult to correct. If one end of the 
resistor is grounded, and if care is taken to keep 
the resistance as low as possible, the capacitance to 
ground within the resistor merely shunts the resist- 
ance, contributing to the effective residual in- 
ductance, /;, which is generally insignificant. 

The term (R}?—%)Cs deserves special attention 
because it is directly dependent upon the effective 
series resistances of the inductors. As an example 
of the effect of this term, consider that, in the meas- 


urement of a 10 h inductor, Rp=10,000 ohms, 
Ry=8,000 ohms, and Cs=100 pf. Then the term 
(R7— RX) Cs=0.0036 h, nearly 0.04 percent of the 


quantity measured. This correction, being inde- 
pendent of frequency (if skin effect may be neglected) 
exists even at low frequencies, at which most of the 
other correction terms are negligible. 

If the time constants, Tp, Ts, and Tr, are 1 psec or 
less and De is 0.001 or less (these are conservative 
but reasonable estimates for commercially available 
components), and the angular frequency is 10,000 
radians per second, the magnitude of K is less than 
0.0003. For a 1 percent difference in inductors, if 
this term were neglected the resulting error would 
be only 3 ppm of the measured inductance. 

Thus, the equal-substitution method can be em- 
ployed without the necessity of using a Wagner 
ground providing care is taken to keep the stray 
capacitance, Cs, reasonably small and nearly con- 
stant and to apply corrections to the measured 
values to offset errors resulting from the difference 
in the internal resistances of the inductors. 


8. Effect of Frequency Difference 


It has been assumed that the frequency of the 
a-c supply is stable. If the frequency, wy, when 
the standard inductor is in the circuit is not equal to 
the frequency, w, when the test inductor is connected, 
there can be an error resulting from the stray capaci- 
tance, Cs. <A difference of frequency between the 
test balance and the standard balance will necessitate 
the addition of the term 


CL 7 Ly(w?— w%) 
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to the right side of eq (18), and since w~wy in any 
practical case, this can be factored to produce 
2wl ‘sh rby(w Wy). 


Since Ly~Ly for the equal-substitution method, 
the additional term to be added to eq (19) is 


If w=10,000 radians per second, wy=10,010 
radians per second (a 0.1% difference), and Cs=100 
pf, for the measurement of a 10 h inductor by the 
equal-substitution method the magnitude of the 
correction term is 0.002 h, an error of 0.02 percent. 
It is unportant to realize that this analysis of the 





effect of frequency variation does not account for 


the change of effective series inductance with fre- 
quency resulting from eddy currents, skin effect, or 


distributed capacitance within the inductors. 


9. Calibration of Inductors 


In recent years standard inductors having good 
stability have become-available in a wide range of 
nominal values. These have made accurate meas- 
urements by the equal-substitution method feasible 
and convenient. Several complete sets of standards 
are maintained in the laboratories of NBS both at 
Washington and at Boulder, and are used regularly 
for calibrating similar standards submitted for certi- 
fication. The values assigned to NBS working 
standards ultimately depend upon a computable 
inductor or capacitor. 

The equal-substitution method can be utilized 
with a variety of a-c bridge circuits. The essential 
requirements are that the bridge have good short- 
time stability, adequate resolution, and means for 
externally equalizing the storage factor (Q) of the 
inductors being compared. An advantage of the 
equal-substitution method is the lessened need for 
accuracy of adjustment of the bridge components 
because the bridge is used merely to measure small 
differences. A number of commercially available 
bridges embody the Maxwell-Wien bridge circuit to 
which the analysis given in this paper is primarily 
devoted. It is almost invariably necessary, however, 
to improve their resolution by adding a calibrated 
variable capacitor in parallel with those built into 
the bridge. 

The resistor, r,, is partially imductance com- 
pensated and is connected in series with the cable 
connecting the inductors to the bridge. It has been 
determined that the existing variations of inductance 


are negligible relative to the inductance being 
measured by this method. This resistor may be 
regarded as serving the purpose of externally 


equalizing the Q of the standard and test inductors. 

It is seldom necessary to apply a correction for the 
term (2}—R)Cs in eq (19). On rare occasions 
when this term is significant, a crude measurement of 
Ry, and Ry is adequate. 





The resistance product (range) in the bridge is set 

that the difference between the test and the 
standard inductors can be accommodated by adjust- 
ment of the externally connected variable capacitor 
and resistor, the other bridge controls being left 
unchanged. 

The analysis of errors in the Maxwell-Wien bridge 
circuit described in this paper was carried out as part 
of the investigation of the feasibility of adopting the 
equal-substitution method for the rapid and conven- 
lent measurement of inductance. 


sO 


The use of the equal-substitution method for the 
precise comparison of inductors at the Electronic 
Calibration Center, NBS, Boulder, Colo., was 
instigated by Chester Peterson, NBS, Washington, 
D.C. The author is indebted to Mr. Peterson for 
helpful comments and suggestions pertaining to this 
paper. 


10. Appendix: Derivation of Equations 


Equation (3) is derived from the equation of the 
impedance of Ze and Zz in parallel. 


= l wl! 

Z re R T jool orl j 
: l wl (orl ty) 
R-+-jol wrC?+ 1 
Lor? C?+-wC (or C- ete ~jol) 

(R- qul) (1. torr) 
Therefore 

F (R+-jol) (1+???) 

Lx = wa : Y Y : ) : = 
L+w??*C* + aw (ar lex + qul) 


Equation (8) is obtained by substituting eqs (3) 


. . -= . . ? 
(4), (5), and (6) into eq (7), giving 


(1-0?) (R+-jol)[(Ry+rr)+jo(L+lz)] 


[1+-wr?C0?+-w0 (wrC+-7) (R +-jal)| 


(Rp+julp)(Rs 


+-juls), 


from which the equation of imaginary components is 


(1+-w*r?C?)[R(L+1,)+(R_+r,)l—Rslp—R pls] 
O[RR pRs+o(RR pls+RRslp+R pR sl)orC 
—w(Rlpls+R pllst+R sll p)—o'll pl wr C}. 
Equation (16) describes the effect of stray capac- 
itance. Figure 3a represents an inductor, L7, having 


an internal resistance, Rr, connected to the bridge 
terminals by a cable having shunt capacitance, C3. 
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The impedance of this circuit, Z, is derived as follows | the two conditions of balance; hence, the assump- 
tion that A is constant is justifiable. In the above 


1 1 ee equations since 
ZR; +jol, 1 CR pRs~Ly 
and 
Z (Rrt+jol 7) (1—w*L rl 's—joR rf 's) CyR pRs~Ly 
s— — _ 5 > - - : - “ , 
(lw LC 54+ joR Cs) (1—w*® LC s—joR Cs) and 


(Ry+rz,) ~(Rint+rry) 





grt qo Lr(1 wT rl 's) —R7C 's] 


(1l- wl rfl 's)?+4 wR %, OR pR sR 7 's L T * R7z( 's and 
Bs hd — CyR pR sR XCsy/ Ly ~ RXC sw. 
In the equivalent circuit shown in figure 3b vR pR sRRCsy/Iy = RRC gy 
‘ ie fi Therefore 
Z R-4 joL, : 


; me ; Lr—Ly=|(C—Cy)R pRs—R pR go*(CL 7 Cs 
and from the two preceding equations it is evident rte iii sie or 


that 





—CyLyCsy)4 R7Cs RC snl + Ky) (dz, lin). 
L-(1 


L a: 





The factor (CLrCs—CyLyCsy) can be expanded as 
follows: 


and, with negligible error, (CL Cs—CyLyCen) =(C—Cy)(LrCe+LyCen) 


Oh g—CLiCen. 


~ 


9 y RFCs 
Lr ( ] wr —— — ) 
od . 5 ; : 
(1 — 2? LCs +o? R2C2) | If Cs and Csgy can be reduced to zero, this factor 
becomes zero, and the equation is greatly simplified. 
However, if Cs is not zero, but if Cs approximately 
equals Csy, 


It should be noted that the discarded term, 
w'L2C2, in the denominator is always much smaller 
than the term, 2w°Z7Cs, which is usually much OL C.—CeL (C—-C\ LOC 
smaller than 1 in practical inductance measurements a) ee ees eee 
at low audiofrequencies. Since the terms in paren- | LOvLeCs—CL Cen 


theses, other than 1, are small, this equation can be | and 
solved for Ly and simplified by neglecting higher 
order terms. Lr—Ly={(C—Cy)RpRs{1 —w?(L 7+ Ly) Cs] 


R20. a ee ee ee ee) ee 
Lr=L ( l—w* L7C'st+ 7 = w? R208 ). ltphisw*(CyLrCs—CLyxCsy) 
4T | } (R? R2)Cs (1-4 k)—(l,—l,y), 
Equation (18) is obtained from the following two pia ; sl 
equations which are derived from eq (14) modified and with negligible error the substitution CyRplis 
3a . Ly and CRpR;s=Lyr can be made in the small 
according to eq (16) ; ae : 
.; correction terms. The second term in braces can be 


R20 reduced to —wLrly(Cs—Cey). If the terms under 
‘ OR PRs (1 FCs aie 1 kK’) | consideration are not relatively small, the approxi- 
Lr mations are not valid. 
RR | Neglecting second order terms, the above equation 
, Jiplis 


(Rr+rz)tr+ can be simplified, giving eq (18): 


R (tp+Ts)—l, 


and 


] Is Ln t(C’' Cn \RpRf1 w*( Le T DCs T Kk] 
twLien , 
Ly=CyR pRs (1—«*LyCsy+—F™ ) (1+-K) 


- 0 Lp Ly(Cs— Cay) + (Ri — Ri) Ca— (1 —lzw). 


R pk s 


(Rr T Trn)TrR+ R (Tp TTs 
i 
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precise determination of field strength of 


equipment, which is battery- 


operated, contains provisions for the reception of VLF signals and the generation of standard 


fields to an 
The 
antennas positioned coaxially at a spacing 
technique was developed for use at VLF, i 
loop antennas, 
or similar 


accuracy of 5 percent. 


generators, 
applications. 


and voltmeters, 


1. Introduction 


With the increased use of the very low frequency 
portion of the radio spectrum, an urgent need exists 
for the accurate calibration of receiving antennas 
operating in this frequency range. Conventional 
methods of calibrating loop antennas require that 
they be placed in a radiation field. In other words, 
the separation between the test source and the 
receiving antenna is large compared to the wave- 
length. At VLE this requirement becomes imprac- 
tical; for example, the wavelength at 10 ke/s is 30 
km. <A test source placed at distances greater than 
this may require several kilowatts of power to pro- 
vide a signal that will compete with the atmospheric 
noise level. 

It is suggested by Greene [1951] that a free space 
radiation field could be simulated adequately at the 
receiving antenna if the test source were also a loop 
antenna and placed coaxially with, and close to, the 
receiving loop antenna which is to be calibrated. 
The magnitude of the equivalent radiation field at 
the receiving loop antenna is then calculated in 
terms of the geometry and configuration of the two 
loops. Such an equation has been stated by Greene 
[1951] without proof. 

It is the purpose of this paper to investigate the 
general applicability of Greene’s equation, to outline 
a technique used for the precise and accurate deter- 
mination of field strength using this equation, and 
finally to describe the equipment used to perform 
these measurements. 


2. Theoretical Basis 


The necessary equations for calibrating a loop 
receiving antenna with a coaxial test transmitting 
loop are first considered. The distance of separation 
is arbitrary, and the dimensions of the loops are 
assumed small compared to the wavelength. The 
situation is illustrated in figure 1. 


! Contribution from Central Radio Propagation Laboratory and Radio Stand- 
ards Laboratory, National Bureau of Standards, Boulder, Colo. 
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Both the receiving end transmitting antennas are loops. 
: field strength is determined in terms of a quasi-static magnetic field with the two loop 
f approximately 
it can be used at higher frequencies for calibrating 
and for det 


two Although the 


meters. 


ermining effective heights of antennas, 
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Simplified view of two coaxial loops. 


The actual value of the voltage v induced in the 
receiving loop may be expressed 1 : terms of the max- 
imum equivalent electric field, H,, that would exist 
in a free-space radiation field. The relationship for 
an electrically small loop is simply 


9 
. rd aT / 
v Bark, where p=>? (1) 
but v=JZ,,, where Z, is the mutual impedance 
between the loops. Where all dimensions are elec- 
trically small and currents in the loop are constant, 


we have exactly 
PY... ; 
~5) K®— EW | (2) 


Zm=(2ar2) due : y ( 1 
(“xe (er Tee...), 


_s ar 2 
JONG 
(3) 


where k | aris . 
(7; tT>)*-+¢ 
mitting and receiving loop radii, respectively, and d 
is their center-to-center separation. A(k) and E(k) 
are elliptic integrals as tabulated by Dwight [1947]. 
The above equation for Z,, may be obtained from a 
number of different sources, including Smythe [1950] 
who gives a formula for the ¢-component of the 


1 
p and 7; and r, are the trans- 
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vector potential, Ay, of a loop carrying a constant 
current J. The mutual impedance Z,, is then 
calculated by integrating L,——jwA, around the 
receiving loop. The length of this path is 277». 


Thus, 
io 
xf 1 +5h 


Neglecting terms in the series above k’, it 
that this expression can be reduced as follows 


60ar2] 
|d?+-(r; +15 )2}8 2 


E,= 


is seen 


@,% 


3 allo wane ; rly . 
‘a Ra 3-2 term: in ( p ) 
|d? TT; 


| 1 2 tite. 8. TRIS 
1g rile] 
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a ay : 2/8 x[14 terms in (> ) | (5) 
+ri+rsy d 


The last expression is the denominator in Greene’s 
equation [1951]. If both loops have radii compa- 
rable to d, the quantity & is not small compared to 

unity, and higher terms must be retained in the | 
series. At higher frequencies where d/\ cannot be 

neglected compared to unity, a correction factor 

may be applied which is the ratio of the coaxial 

magnetic field of a small loop at wavelength \ to the 

infinite wavelength case. This is known 
1+-72rd/d and has a magnitude 


k -( ie 


as stated by Greene [1951]. It is strictly valid for 
d ‘7, but it is only a minor second order correction 
at medium frequencies and is completely negligible 
at VLF. 

Equation (4) is valid for determining the equiva- 
lent free-space electric field strength if 7, r2., and d 
are electrically small compared to \. Since this 
equation is strictly valid only in free-space, the 
effects of neighboring objects have been verified 
experimentally by moving metallic objects of various 
sizes near the loop antennas and observing the 
effects on the system. If the nearest sizable objects 
are separated from the loops by approximately two 
to three times the loop spacing, d, the resulting 
error will be less than about 1 percent. The meas- 
urements should not be performed in a_ shielded 
room, and the measurement site should be free of 
large metallic objects such as overhead power lines 
or similar objects. 


to be 


2ard 


3. Calibration Procedure 


Two basic methods are applicable to the measure- 
nent of field intensity [IRE Standards, 1942]; one 
is the standard-antenna method, and the other is 
the standard-field method. In the former a standard 
receiving antenna is used so that the voltage pro- 
duced in it by a field of given intensity and polariza- 
tion may be computed. The usual voltage measured 
is not the open-circuit antenna voltage but is related 


to it through a voltage transfer ratio. It is some- 
times difficult to determine this voltage transfer 
ratio to a high degree of accuracy. 

In the standard-field method, a standard field is 
set up by a local transmitting arrangement, and the 
field at the receiving antenna is computed from the 
dimensions of the receiving and transmitting an- 
tennas, the separation distance between the antennas, 
and the current in the transmitting loop. In this 
technique the equipment of a suitable 
receiver equipped with a meter to indicate the rela- 
tive strength of the received field. The strength of 
the received field is determined by comparing the 
receiver output voltage produced by the field under 
measurement with the voltage produced by a stand- 
ard field. 

Each of the basic methods offers certain advantages 
and disadvantages dependent upon the application. 
In the case under consideration, the standard-field 
method was adopted because of the straightforward 
calibration method which can be evolved using co- 
axial loop antenna arrangements as discussed above. 

The calibration equipment is represented in figure 
2; that to the left is used to receive and indicate the 
relative amplitudes of fields, while that to the right 
is used to produce VLF fields of appropriate fre- 
quency and amplitude. In calibrating, the deflec- 
tion of the receiver output meter in response to the 
unknown VLF field is obtained with the receiving 
loop in the plane of propagation. The receiving 
loop is then rotated 90° to minimize the receiver re- 
sponse to this field. The transmitting loop 
positioned coaxially with the receiving antenna at : 
separation distance of 1 or 2m. Without changing 
the receiver tuning or sensitivity, the frequency and 
output of the signal generator are adjusted to repro- 
duce the same receiver output meter indication that 
was obtained in response to the unknown VLF field. 
The equivalent value of the radiation field incident 
at the receiving antenna is calculated from eq (4) 
above in terms of the antenna radii, their separation 
distance, and the RF current in the transmitting an- 
tenna. The strength of the unknown field is equal 
to the equivalent radiation component of the locally 
generated field. 

In figure 2, the “Metering Equipment” block is 
svmbolic of the method used to determine the RF 
current in the transmitting loop. The technique is 


consists 


Is 


represented in detail in figure 3. Balanced RF 
voltages, obtained from the RF signal generator and 
Balun balaneced-to-un balanced transformer), are 
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FIGURE 2. diagram of loop antenna calibration 
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FiaureE 3. Calibration of the transmitting loop antenna, 
connected to the single-turn unshielded transmitting 
loop. A vacuum thermoelement is mounted in the 
top-center of the transmitting loop as indicated. 
The RF current introduced into the loop flows 
through the thermocouple heater element. The d-c 
output voltage of the thermoelement is indicated by 
the d-c millivoltmeter. The RF current then 
disconnected from the loop antenna and replaced by 
a direct current. The direct current is adjusted to 
produce the same thermoelement d-c output voltage 
that was attained in response to the RF current. 
Since equal values of rms RF current and direct 
current produce equal thermoelement output volt- 
ages, the value of direct current read from the milli- 
ammeter, figure 3, is used in eq (4) in calculating the 
field strength at the receiving loop antenna. 

An REF micropotentiometer? [Selby, 1953] was 
incorporated into the calibration equipment to pro- 
vide an accurate of RF voltages. The 
arrangement used for calibrating the micropotenti- 
ometer and establishing known RF voltages is shown 
in figure 4. The thermoelement of the micropo- 
tentiometer is the same type used in the transmitting 
loop. The switching arrangement indicated permits 
either RE or direct current to flow through the 
thermoelement and the micropotentiometer resistor 
element. To establish an accurately known RF 
voltage across the micropotentiometer, the current 
required to produce the desired voltage is calculated. 
A direct current of this magnitude (as indicated by 
the milliammeter) to flow through the 
thermocouple and resistor element of the micropo- 
tentiometer. The d-c output of the thermoelement 
in response to the direct current is observed. An 
RF current of appropriate frequency is then sub- 
stituted for the direct current and adjusted in mag- 
nitude to produce the same thermoelement output 
voltage that was obtained in response to the direct 
current. The required RF voltage then exists 
across the resistive element of the micropotentiom- 
eter. The micropotentiometer, having an internal 
impedance of a few milliohms, is usually used to 
calibrate equipment having input impedances rang- 
ing from 50 ohms upwards. In such the 
standard voltage across the micropotentiometer is 
not significantly altered when the equipment under 
calibration is connected. 


is 


source 


Is caused 


cases, 


2 Thisis simply a low-impedance source of RF voltage adjustable from 1 to 105nv 
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Ficure 4. Calibration of the micropotentiometer. 


4. Description of Equipment 


A photograph of the prototype of the field intensity 
and calibration equipment is given in figure 5. As 
viewed from left to right, the individual items are: 
The receiving loop antenna shown connected by a 
cable to a box containing a balanced-to-unbalanced 
electrostatically shielded transformer. The third 
unit from the left is a transistorized VLF receiver. 
The next item is a commercially made transistorized 
VLF signal generator. The transmitting loop an- 
tenna is shown mounted on a supporting pedestal 
with the input terminals near the bottom and the 
thermoelement and its output at the top. The small 
box in front and to the right of the transmitting loop 
contains a transistorized final amplifier for the signal 
generator. The large box at the extreme right is the 
calibration unit containing the laboratory type d-c 
milliammeter and d-c millivoltmeter, the micropoten- 
tiometer, switching circuits, etc. A photograph of a 
collapsible frame to support the transmitting loop 
antenna (to the left) and the receiving antenna (to the 
right) is shown in figure 6. The spacing between the 
loop antennas in this photograph is 2m. The frame- 
work is made of short sections of tubing which may 
be disassembled into a “bundle of sticks” 
transport. 


for easy 








FiGcuReE 5. 


Field intensity and calibration equipment, 











Fiacure 6. 


Loop antennas mounted on a collapsible frame for 
field use. 


Although all of the equipment shown is battery- 
operated and has been used in field work, no serious 
effort was made to minimize the size or weight of 
this prototype equipment. Obvious improvements 
could be made, such as reducing the diameter of the 
receiving loop and the input transformer, positioning 
the final amplifier within the signal generator, minia- 
turizing the calibration unit, ete. 

The receiving loop antenna is electrostatically 
shielded and balanced with respect to ground. 
The input transformer is a balanced-to- unbalanced 
unit with an electrostatic shield betweeen the 
primary and secondary windings. A high degree of 
balance was achieved in the loop input circuit to 
minimize undesirable unbalanced currents. 

The transistorized, battery-operated VLF receiver 
shown in the photograph is an early model developed 
at NBS. An improved model has since been devel- 
oped which has a sensitivity of about 100 »v/m and 
a bandwidth adjustable from 40 to 400 c/s. The 
receiver is tunable from 14.5 ke/s to 22 ke/s and i 
provided with a voltmeter to indicate the antes 
amplitude of received signals. In the experiment 
cited, voltmeter readings were made during intervals 
when station NLK was transmitting a ste ady carrier 
signal. It is intended to describe this receiver in a 
future publication. 
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FiGuRE 7. Schematic of calibration unit. 
A schematic diagram of the calibration unit is 


shown in figure 7. The apparatus for the generation 
of standard fields (represented in fig. 3) and for the 
calibration of the micropotentiometer (represented 
in fig. 4) are combined in this unit. The d-c milliam- 
meter and millivoltmeter are high quality laboratory 
instruments with accuracies of 0.5 percent. Thermo- 
elements of 5- to 100-ma capacity were used in the 
micropotentiometers, which enabled standard volt- 
ages to be generated within a range of 100 to 50,000 
uv with an accuracy of better than 3 percent. This 
type of thermoelement when calibrated with direct 
current can be used as a current-measuring device 
up to 300 Me/s to an accuracy within 1 percent. 

The transmitting loop consists of a single-turn 
unshielded loop with a 100 ma_ thermoelement 
mounted in the top-center. The balun transformer 
has a voltage unblance that is less than 2 percent of 
the balanced output voltage, and an insertion loss 
that is less than 0.3 db. This system can be used to 
establish standard fields with excellent precision and 
an accuracy that is better than 5-percent. 

The transmitting and receiving loop antennas, 
with radii of approximately 0.1 and 0.35 m, respective- 
ly, were spaced 2.0 m from each other. Under these 
conditions the bracketed series in eq (4) yields a 
correction factor of approximately 1.028. This repre- 
sents an increase of approximately 3 percent over 
the value of the field if the correction is neglected. 

The equipment described in this paper was de- 
signed to measure field strengths in the range from 
230 to 23,000 uv/m at frequencies between 15 and 
20 ke/s. The range of frequency coverage can be 
extended by proper selection of the balun trans- 
former, the signal generator, and the receiver. The 
micropotentiometer can be used over a voltage range 
of 1 uv to 200 mv at frequencies up to 300 Me/s. 
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The ranges of operation of this equipment can be 
extended in frequency and amplitude for use in other 
applications. Excellent precision is obtained with 
this system with an accuracy not before attainable. 
All calibrations are performed with direct current. 
This eliminates many sources of error, and the d-c 
measuring instruments can readily be checked against 
reliable d-c working standards. This unit was 
recalibrated after an extensive field trip, and the 
calibration of the instruments did not vary more than 
0.2 percent. 
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Location of the Plane of Best Average Definition With Low 
Contrast Resolution Patterns 


Francis E. Washer and William P. Tayman 


(April 28, 1961) 


The plane of best average definition is located for each of several airplane-camera lenses 


using two types of low contrast test pattern and two emulsions. 


A low contrast pattern 


composed of dark lines on a light background and the reverse pattern consisting of light 


lines on a dark background are used, 


The results of measurement indicate that the position 


of the plane of best focus and the numerical magnitudes of the root mean product mean 
VR; 7, Value of the resolving power are not significantly affected by this reversal of contrast 
In addition, the results obtained using low contrast targets are compared with those obtained 


with high contrast targets 


Hence the position of the selected focal plane remains invariant 


although the values of the measured resolving power are substantially higher for the high 


contrast targets. 


1. Introduction 


In specifications dealing 
cameras, great emphasis is 
graphic resolving power of the lenses used in the 
camera. The method of measurement to be used 
in evaluating the resolving power is usually described 
in careful terms, with the type of test chart, the 
characteristics of the registering emulsion, the 
physical conditions of test, and the manner of report- 
ing results all given in detail. A 1950 military 
specification, MIL-STD-—150 [1]! attempted to 
establish standard practices in the evaluation of lens 
performance. A recent revision of this compre- 
hensive specification, MIL-STD-—150A [2] was issued 
in 1959. At the Washington meeting of the Inter- 
national Society of Photogrammetry in 1952, a 
specification dealing with the calibration of photo- 
grammetric cameras was adopted for trial and 
discussion [3]. 

At the time the International Specification for 
the calibration of photogrammetric cameras was 
drawn, it was not possible to standardize on a single 
type of resolving power chart. Consequently several 
were listed as suitable for use; these include the 
three-line chart used by the U.S. Air Force [2], the 
Cobb two-line chart used in Great Britain and the 
annulus chart used in Canada [4]. These charts 
were of low contrast except that the three-line 
chart of the Air Force was provided in both high 
and low contrast versions. The three-line chart of 
the National Bureau of Standards [5] was not included 
as it has been available for only a few months prior 
to the drawing of the specification and its properties 
were not widely known. 


with airplane mapping 
placed upon the photo- 





1 Figures in brackets indicate the literature references at the end of this paper, 





In the eight years that have passed since the 
adoption of this specification, agreement on a single 
test chart has not yet developed. A few publications 
have appeared [6, 7] that permit an estimate of the 
differences likely to be found in comparing various 
types of charts. 

Because of the differences that exist among the 
various test charts used in various laboratories, a 
study was initiated at the National Bureau of 
Standards in which the values of resolving power 
for the same group of lenses were determined for 
au variety of test charts, target contrasts, and 
emulsions. The results of this study should be 
useful in estimating the probable values of resolving 
power that would be obtained for a lens with a 
given type of chart when results are available for 
the same lens with a different type of chart. 

This is the concluding portion of a three part 
investigation. The results for the first part dealing 
with various indices and methods used in locating 
the plane of best average definition have been 
reported in Photogrammetric Engineering [8]; the 
results for the second part which comprises results 
for lone-line, short-line, and annular charts are 
reported in J. Research NBS [9]. The results of 
the third part which deals with variations arising 
from differences in target contrast for two emulsions 
are reported in the present paper. 

Measurements of photographic resolving power 
at the National Bureau of Standards are customarily 
made using a long-line three line high contrast test 
chart having dark lines on a light background with 
the images registered on a high contrast emulsion 
which has a higher resolving power than the emulsion 
commonly used in aerial photography. In many 
other calibration laboratories, the tendency is to 
use test charts having light lines on a dark background 
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with the images registered on a medium contrast 
emulsion with resolving power comparable to that 
being currently used in ‘aerial photography. 

In the present study, values are obtained for two 
emulsions and three targets of differing contrast. 
Results are given for four lenses of types commonly 
used in photogrammetric cameras. 


2. Method of Measurement 
2.1. Test Camera 


The National Bureau of Standards precision lens 
testing camera [10] was used in making the negatives 
from which the values of resolving power are deter- 
mined. Its construction and mode of operation 
simplify the process of evaluating the resolving 
power of a lens in a series of focal planes spaced 
at definite intervals along the optical axis. The 
precision lens testing camera is equipped with 10 
collimators spaced at 5° intervals and spanning the 
range from 0° to 45° 
placed in the focal plane of the collimators and 
umaged by the lens under test permit the study of 
quality of imagery for a series of focal planes of the 
test lens. By using one type of chart in 5 of the 10 
collimators (designated the’ odd-numbered  colli- 


Resolving power test charts | 





mators) spaced at 10° intervals and another type | 


in the remaining five collimators (designated 
even-numbered collimators) it is possible by proper 
manipulation of the various controls to record the 
imagery of a lens under test on a single photographic 
plate for each of two types of chart in the same series 
of focal planes for the same angular separations, 
8, from the axis of the lens under test. 


2.2. Test Charts 


Three types of test chart are used in the present 
investigation. All three are based upon the NBS 
resolution chart of 1952 but differ in contrast between 
the lines and spaces. Each chart has the same 
range of sizes of the test object with the size pro- 
ceeding in a geometric progression with ~/2 (or 
1.1892) as the common ratio. The same size test 
chart is used in the various collimators so that it is 
necessary to apply the ‘cos’ and “cos”’ corrections 
to images formed off-axis. The various types are 
described in the following section. 


a. High Contrast Test Chart (HDL) 


This is a composite chart based upon the NBS 
resolution test chart of 1952 which is described in 
detail in NBS Circular 533 [5]. It is a long-line 
three line chart with dark lines on a light background. 
The difference of log luminance between the dark 
lines and the light background exceeds 2.0. 


b. Low Contrast Test Chart (LDL) 


This chart is identical in all respects to that 
described in (a) except that the contrast is low 
between the dark lines and the light background. 
The difference of log luminance between the dark 
lines and the light background is 0.20 +0.02. 


the | 


c. Low Contrast Test Chart (LLD) 

This chart is identical to that described in (b) 
except that the contrast is reversed, that is, the 
difference in log luminance between the light lines 
and the dark background is 0.20 +0.02. 


2.3. Test Negative 

In making the test negative, the lens under test is 
initially so alined that its optical axis is parallel to 
and approximately coincident with the axis of the 
collimated beam emergent from the first collimator. 
The lens is adjusted along the bench to a location 
such that the collimated beam at 40° from the axis 
fills the front aperture of the lens under test as viewed 
through the lens at an inclination of 40° from the axis. 
The plate holder is adjusted to a position such that 
the front surface of the emulsion is in the plane of 
best visual axial focus for the central row of images 
to be registered on the plate. The plate holder is 
then moved to a position 1.05 mm nearer to the lens, 
where an exposure made by illuminating the 
reticles in odd-numbered collimators, which thus 
records the imagery on the plate at 10° intervals for 
the range of angles from B=0° to B=40°. This 
process is repeated with the plate moved 0.15 mm 
farther from the lens until 15 exposures have been 
made with the last for the plane 1.05 mm farther from 
the lens than the position of best visual axial focus. 
Between each exposure, the plate holder is also moved 
downward by an amount sufficient to avoid super- 
position of successive rows of images. The foregoing 
operation registers the imagery for the long-line high 
contrast patterns on the plate. The plate holder 
is then returned to its initial position for which the 
emulsion surface is 1.05 mm nearer to the lens than 
the plane of best visual axial focus; the plate is dis- 
placed sidewise in its holder approximately 12 mm; 
and the entire bench on which the camera is mounted 
is rotated on its pivot by 5° so that the axis of the 
lens is now parallel to and nearly coincident with the 
axis of the collimated beam emergent from the second 
collimator. The foregoing procedure is then re- 
peated with the exception that exposures are made by 
illuminating the reticles in the even-numbered col- 
limators, which again registers imagery on the plate 
at 10° intervals for the range of angles from B=0° 
to B=40°. This process registers the imagery for 
the long-line low contrast patterns on the plate. 
The exposed plate is then processed to form the 
finished negative from which values of resolving 
power for tangential lines (73) and radial lines (Rs) 
are determined. 

All exposures are made with the reticles illuminated 
by light from tungsten lamps with Wratten K- 
filters between the light source and test charts. 
Neutral filters are used to adjust the intensity of the 
light reaching the plate so that the final resulting 
optical density of each image on the negative is 
approximately the same for all values of 8. The 
two types of photographic plates used in this work 
were Eastman Spectroscopic VF which has a fine 
grained panchromatic emulsion and Eastman Super 
Panchro Press, T ype C. The plates were developed 
for 3 min. in D-19 de »veloper at 68 °F. 
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2.4. Reading the Negative 


The negative images were examined with a micro- 
scope using powers ranging from 30 to 50x. The 
criteria for conside ‘ring a particular line pattern to be 
the finest resolved were that all coarser patterns were 
resolved and that the number of lines in a given 
pattern was the same as that of the corresponding 
pattern in the object. 


3. Results of Measurements 


The results of measurement on four wide-angle 
lenses are reported and analyzed in this study. 


Two lenses, designated Nos. 1 and 5, are essentially 
distortion-free; the other two, designated 7 and 8, 
have moderate amounts of distortion. Measure- 


ments of resolving power were made at 10° intervals 
from 0° to 40 

Negatives were made with both VF 
sions for each of three target contrasts for each of 
the four The measured equivalent focal 
length for the plane of best visual axial focus was 
the same for each emulsion and each type of target 
pattern. The maximum range of change in measured 
equivalent focal length did not exces ed + 0.02 mm 
for a given type of lens. 

The results of measurement are shown graphically 
in figures l and 2. The 
resolving power y Rs Ts are shown plotted against the 
separation, Af, from the plane of best visual axial 
focus. Positive values of Af denote positions of the 
image plane farther from the lens than the plane of 
best visual axial focus. In figure 1, values of ¥R, T3 
obtained with VF emulsion are shown for values of 8 
ranging from 0° to 40°. The lowest frame in each 
column shows values the root mean product 


and SP emul- 


lenses. 


values of the geometric mean 


of 


mean [8] V/s 7's versus Af. The curves, numbered 1, 
show the results obtained with the high contrast 
target having dark lines on a light background. The 
curves, marked 2, show the results for the low con- 
trast target having dark lines on a light background. 
The curves, marked 3, show results for the low con- 
trast targets having light lines on a dark background. 
In figure 2, comparable results are shown for the 
same four lenses using SP emulsion. 

Values of ys 73 are shown in these graphs as this 
quantity of primary interest in evaluating the 
image forming qualities of a lens intended for use in 
an airplane mapping camera. In two earlier papers 
[8, 9], values of Rs and 7's are shown for lens No. 


is 


4. Location of the Plane of Best Average 
Definition 

The plane of best average definition was located 
for each set of conditions by the maximum value of 
the index viz 73. The plane so located is likely to 
be slightly farther from the lens than that located by 
graphical analysis [8]. For the purposes of the pres- 
ent paper, the use of the maximum of V/?, 7’; is satis- 
factory. The values of 


the index Vs 73 for the 


| nearer to the lens for 





range of Af extending from 0.45 mm nearer to the 
| lens to 0.30 mm farther from the lens than the plane 
| of best visual axial definition is shown in table 1 for 
the four lenses for all of the various conditions. The 
values listed are the averages for several runs ranging 
from 1 to 6. The total number of runs for the high 
contrast target is equal to the sum of the runs for 
the two low contrast targets. This occurs because a 
single negative always contains a record of the im- 
agery for the high contrast target and for one of the 
two low contrast targets. In table 1 the underlined 
value of the index in each column is the maximum 
for the indicated set of conditions. 


5. Comparison of Results Obtained With 
the Three Charts 


The principal points to be considered in the 
analysis of the results of measurement are possible 
differences in the location of the focal plane for the 
various conditions, the reduction in average resolving 
power that may result from such differences, relative 
magnitudes of average resolving power for the two 
emulsions and effects of contrast. 
cussed in the present section. 


These are dis- 


5.1. Relative Positions of Focal Plane 


The relative positions of the focal plane for the 
various conditions are indicated for the four lenses 
in table The displacements for each of the con- 
trasts and emulsions from the focal plane of best 
average definition for the high contrast target using 
VF emulsion of the corresponding planes are shown 
in table 2. It is clear from this table, that differences 
in contrast of target have little or no effect upon the 
location of the focal plane of best average definition 
for either emulsion. Onlv in the case of the low 
contrast target having light lines on a dark back- 
ground (LLD) used with SP emulsion is there any 
indication of a possible effect on the location of the 
focal plane and even in this instance the apparent 
shift may be a result of random error. 

On the average, the plane of best average definition 
for the various contrasts appear to be 0.15 mm 
the SP emulsion than is the 
corresponding plane for VF emulsion. This difference 
is small and may be real. However, it would be 
unlikely to present serious cause for concern in the 
locating of the plane of best average definition. It 
is likely that the average plane selected using VF 
emulsion will be satisfactory if this plane be used for 
photography with SP emulsion. Likewise the plane 
selected using SP emulsion would be reasonably 
satisfactory for photography with VF emulsion. 


5.2. Reduction in Resolving Power From the 
Maximum 


If the differences in the relative location of the 
focal plane are regarded as genuine, it is instructive 
to determine the magnitude of the effect of these 
differences on resolving power under the various 
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TaBLe 1. Variation of v R; Rp Ts. with separation, Af, from the 


plane of best visual azial focus for four lenses 


— 

Values of ¥ Rg Ts in lines per millimeter are given for two emulsions and three 
target contrasts. The columns, designated HDL, show values obtained with a 
high contrast target having dark lines on a light background. The columns, 
designated LDL, show values for a low contrast target having dark lines on a 
light background. The columns, designated LLD, show values obtained with 
a low contrast target having light lines on a dark background. The maximum 


values of V Rs Ts are underlined. The number of test negatives used in each 
set of determinations is shown in the rew marked n 


Values of ¥ Rg 7's for lens No. 


1 obtained with emulsion 





VF SP 
n 2 1 1 3 l 2 
Af HDL LDL LLD HDL LDL LLD 
mm 
—0.45- _ 39.3 27.2 25. 2 20. 0 14.6 13. 5 
— .30 42.0 28.9 27.8 22. 2 15.3 13.8 
—.15 46.4 31.3 28. 3 22.9 15.4 13. 7 
-00 48. 6 32.4 28.7 21.8 14.9 13. 5 
15, 48.4 299.3 6. 6 19.4 12.7 11.8 
30. 38. 6 23. 2 3. 2 15.6 11.1 10.4 
Values of ¥ Rg 7 for lens No. 5 obtained with emulsion 
VI SP 
n 3 2 1 6 3 3 
Af HDI DL LLD HDI DI LLD 
mm 
—(0.45 30.1 22. 2 19.8 14.4 13.1 
30 34.6 6. 0 21.8 15.1 14.5 
15 39.8 27.6 23. 4 15.6 14. 1 
.00 42 28. 8 22. 8 15.2 13.8 
15 42.1 26. 8 21.4 13. 5 12.6 
.30 36. € 23. 0 17.6 11.7 10.8 
Values of V Rs T3 for lens No. 7 obtained with emulsion 
VI SP 
n 4 2 2 2 1 1 
Af HDI LDL LLD HDI LDL LLD 
mm 
0.45__. 42.6 31.6 9. 6 24.6 16. 1 14.3 
— .30 46.5 33. 0 29. 8 24.8 16. 4 14.3 
—.15 » a 48. 1 33.0 30.0 24.5 16. 2 14.0 
00 46.8 32.3 29. 2 23. 6 16. 2 13.9 
15 43.4 27.2 27.4 21.5 13.9 12.4 
30 36.5 24. 2 23. 4 17.9 12.8 11.1 
Values of ¥ Rz T3 for lens No. 8 obtained with emulsion 
VF SP 
n 3 l 2 2 1 1 
As HDI LDL LLD HDI LDL LLD 
mm 
—0.45 _— 34.5 28. 0 24.2 21.4 13.3 13. 2 
— .30_. . . 36. 28.3 26. 6 21.5 13.3 13.2 
—.15. ‘ = 37.4 28.8 26. 4 21.7 13.7 12.9 
2 ————— 38.9 27.8 25. 9 21.0 13. 4 12.5 
+ ass seiee SME 37.4 25.5 24. 4 18.5 11.7 11.6 
Naa 33. 2 23.2 21.6 17.2 10.8 10.9 








TaBLE 2. Location of the focal plane of best average definition 
with respect to that obtained for a high contrast target using 
VF emulsion for four lenses 


Results are given for VF and SP emulsion and for two types of low contrast 
target. The symbols in the colmn heading are defined in table 1. A negative 
value of A/indicates a position of a selected focal plane nearer to the lens than that 
selected for the high contrast target using VF emulsion. The focal plane of best 


average definition is determined by the maximum value of V Rs Ts 


Separations, Af, in millimeters from plane of best average 
definition obtained with high contrast target for emulsion 


Lens 
Vi SP 
HDL LDL LLD HDL LDL LLD 
1 0.00 0. 00 0. 00 0.15 0.15 0. 30 
5 v0 OO 15 15 15 30 
7 00 00 oO 15 15 15 
x 00 15 30 15 15 30 
A verage 0. 00 0.04 0. 04 0.15 0.15 0. 26 
TABLE 3 Effect on average resol ing power, Rs Ts for various 


contrast targets produced by using focal plane of best average 
definition obtained with VF the 
target 


emulsion for high contrast 


Values of the reduction from the maximum are given in percent for four lenses 


using VF and SP emulsions for one high and two low contrast targets. The 
meanings of the symbols heading the columns are as given in the caption of table 1. 
Reduc nin ©, of ¥ Rg7 from its maximum for 
Lens V F emulsion SP emulsion 
HDI LDI LLD HDL LDL LLD 
l 0.0 0.0 0.0 4.8 3.2 2.2 
5 0 0 3.4 2.6 2.6 1s 
7 0 0 0.0 1.2 1.2 2. 1 
8 0 3. § 2.6 3.2 2.2 


conditions. For example, one may select the focal 
plane of best average definition for high contrast 
targets with VF emulsion and then use this focal 
plane for photography under all conditions of con- 
trast for both emulsions. This has been done and 
the effects on resolving power for the other conditions 
are shown in table 3. From this table it is clear that 
for 3 of the 4 lenses, there is no reduction in resolving 
power when using VF emulsion for the low contrast 
target having dark lines on a light background 
(LDL) while there is a 3.5 percent loss in the case of 
lens No. 4. For the low contrast target having light 
lines on a dark background, the change is zero for 
two lenses and reductions of 3.9 and 2.6 percent are 
shown for the other two. When SP emulsion is used, 
reductions in resolving power ranging from 1.2 to 
5. oy percent occur, with the reduction under 3 percent 

1 7 of 12 instances. One can expect a 3 percent 


variation in the magnitude of y 2373 on the basis of 
observational error, hence in only five of the twelve 
cases is the reduction in resolving power in excess of 
the probable observational error. 

If the plane of best average definition for the high 
contrast target using SP emulsion is taken as a re ef- 
erence plane, values of the reductions in 
power can be determined for the other 


resolving 
condition. 
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The values of these reductions are shown in table 4 
For this condition, the reduction is zero for the low 
contrast target (LDL) with SP emulsion. For the 
VF emulsion, the reduction is under 3 percent in 5 of 
12 cases and does not exceed 5.6 percent for any of 
the remainder. 


5.3. Relative Magnitude of the Resolving Power 


From figures 1 and 2, it is apparent that the re- 
solving power for VF emulsion is higher than that 
for SP emulsion for all values of 8 ranging from 0° to 
40°. Itis moreover clear that this is true for all three 
conditions of contrast. It is not possible to determine 
a numerical value of the ratio of the resolving powers 
that would be invariant for all values of 6 and all 
valuesofAf. Accordingly, average values of the ratios 


based upon the values of V/s 73 shown in table 1 
have been determined for each lens and are shown in 
table 5. In table 5, the ratios of average resolving 
power for the two emulsions and three contrasts with 
respect to the average resolving power for the high 
contrast target obtained with SP emulsions for four 
lenses. From the average values of the ratios given 
in this table, an estimate can be made of the probable 
value of other ratios that are of interest. For ex- 
ample, values of resolving power obtained with the 
low contrast target are approximately two-thirds 
of that obtained with the high contrast target for 
either of the two emulsions. 


TasBLe 4. Effect on average resolving power, ¥ Rg Ts for various 
contrast targets produced by using focal plane of best average 
definition obtained with SP emulsion for the high contrast 
target 


Values of the reduction from the maximum are given in percent for four lenses 


using VF and SP emulsions for one high and two low contrast targets. The 
meanings of the symbols heading the columns are as given in the caption of table 1 
Reduction in % of V Rs T's from its maximum for 
Len VF emulsion SP emulsion 
HDL LDL LLD HDL LDL LLD 
1 4.5 3 1.4 0.0 0.0 0.7 
5 5. 6 1.2 3.9 0 .0 2.8 
7 3.3 0.0 0.7 0 0 0.0 
x 3.8 0 8 0 0 2.3 
TABLE 5. Ratios of average resolving powers for two emulsions 


and three targets of various contrasts 


Ratios are given showing the relative magnitude of the average resolving power 
obtained under specified conditions to that obtained with a high contrast target 
on SP emulsion. Values are given for four lenses, using two emulsions and three 
contrasts. Values of the averages are given together with the average probable 
error PE determination. The meani 
is given in the caption of table 1. 


, of a single 


gs of the symbols heading the 
columns are 


Ratios of average resolving powers for 


Lens V F emulsion SP emulsion 

HDI LDL LLD HDL LDL LLD 
Raia 2.17 1. 42 1, 32 1. 00 0. 69 0. 63 
tics 1.77 22 1, 21 1. 00 . 67 . 62 
a a 1. 93 1. 32 1. 24 1. 00 . 67 . 59 
8... 1. 80 1. 33 1. 23 | 1. 00 63 . 63 
Avert iv esues 1.92 | 1. 32 1. 25 | 1.00 | 0. 66 | 0. 62 

+4.0 | ot 5 ore +2. 6 | +5 








5.4. Effect of Reversed Contrast 

In many specifications that require measurements 
of resolving power of lenses with low contrast targets, 
it is stipulated that the test chari shall consist of 
light lines on a dark background. This is so speci- 
fied because of the belief that the use of a target 
having dark lines on a light background will result 
in lower measured values of resolving power than 
would be found with a target of opposite contrast 
because of the effect of deopeig glare [3]. In figures 
1 and 2, the curves marked 2 show values of Vs Ts 
as a function of Af for the low contrast target basin 
dark lines on a light background; the curves marked 
3 show values of ys; Ts as a fraction of Af for a low 
contrast target having light lines on a dark back- 
ground. It is noteworthy that in most instances, 
the curve marked 2 lies slightly above the curve 
marked 3. Even if one ascribes the difference to 
experimental error, it is improbable that such error 
could reduce all of the observed values of curve 2 


below those for curve 3. From these results, it is 
clear that the values of resolving power obtained 


using targets having dark lines on a light background 
are not lower than those obtained using targets 
having light lines on a dark background when the 
difference in log luminance between line and back- 
ground is 0.2 

It seems aaa therefore, that when resolving 
power tests are made using transparent charts 
located in the focal plane of collimators, the use of 
light lines on a dark background will not, in general, 
vield higher values of the resolving power than will 
the use of dark lines on a light background. In this 
case, the values of resolving power are not noticeably 
reduced by the effect of veiling glare. Although 
veiling glare may affect the values of resolving power 
on the axis which may be the case for lens No. 5 in 
figure 2, it is not likely to do so for extra-axial imagery. 
The effect of veiling glare possibly would be more 
pronounced in using resolution targets on a range 
rather than in collimators which may account for 
the prevalence of the specification of targets having 
light lines on a dark background. 


6. Conclusion 


In this study, the effect of contrast upon choice of 
a focal plane has been investigated using four lenses, 
two target contrasts, and two emulsions. Analysis 
of the results of measurement leads to the following 
conclusions. 

(1) For either of the two emulsions used, the choice 
of focal plane of best average definition is not affected 
by differences in contrast between lines and back- 
ground in the target. 

When a lens images targets where the differ- 
ence in log luminance between lines and background 
equals 0.20, the values of measured resolving power 
are not significantly affected by reversal of target 
contrast. In other words, low contrast targets with 
either dark lines on a light background or light lines 
on a dark background will yield substantially the 
same measured values of aver age resolving power for 
a given lens. 
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FiGcurE 1. Average resolving powers versus position of the image plane for four lenses with VF emulsion. 
{ 9 fT } J ge } Z J 


Values of the average resolving power V Rg 7's are shown at 10° intervals from 0° to 40°. The lowest box in each column shows values of / Rs Ts . Curvel 
shows the result obtained with a high contrast target having dark lines on alight background with contrast between lines and background greater than 2.0 on a densit 
scale. Curve 2 shows the results for a low contrast target (0.2 difference of log luminance between lines and background), having dark lines on a light background, 
Curve 3 shows results obtained for a low contrast target (0.2 difference of log luminance between lines and background) but having light lines on a dark background. 
The zero of abscissas marks the position of best visual axial focus, and positive values of Af indicate positions farther from the Jens 
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FIGURE 2. Average resolving power versus position of the image plane for four lenses with SP emulsion. 
Values of the average resolving power V Rs T's are shown at 10° intervals from 0° to 40°. The lowest box in each column shows values of VR Ts . Curvel 


shows the results obtained with a high contrast target having dark lines on a light background with contrast between lines and background greater than 2.0 on a 


density scale 


Curve 2 shows the results for a low contrast target (0.2 difference of log luminance between lines and background), having dark lines on a light back 


zround. Curve 3 shows results obtained for a low contrast target (0.2 difference of log luminance between lines and background) but having light lines on a dark 


yackground 


The zero of abscissas marks the position of best visual axial focus, and positive values of A f indicate positions farther from the lens. 
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At relative humidities close to a hundred percent, bare film packets are unsatisfactory 


for personnel dosimetry, regardless of temperature. 
siderable protection from excessive humidity. 


1. Introduction 


Fading characteristics of the photographic latent 
image are known to depend on the type of emulsion, 
the temperature, humidity, and chemical contamina- 
tion of the atmosphere, as well as on the type of 
radiation exposure [1].!. Because of the variety of 
conditions under which photographic personnel 
dosimeters are presently used, and because of the 
long periods of time over which they are carried, a 
knowledge of the fading characteristics of the photo- 
graphic latent image is of considerable importance 
in personnel dosimetry. 

Karly fading studies performed at the National 
Bureau of Standards on the latent image formed with 
X-rays were carried out over a period of only one 
week; all films were exposed at laboratory tempera- 
atures to low-energy X-radiation and then stored 
for the prescribed length of time at about 25 °C 
and 45 percent relative humidity [2]. More recently, 
Ziegler and Chleck [3] carried out fading studies 
over a period of 90 days, on films exposed at labora- 
tory temperatures and stored at temperatures 
between 27 and 34 °C, either in air without humidity 
control, or in a dry atmosphere of argon. Tomoda 
et al. [4,5] did latent-image fading studies on 
several types of X-ray film exposed at laboratory 
temperatures and stored up to 30 days at 30° C at 
relative humidities between 40 and 75 percent. 

The aim of the present study was to determine the 
response of one type of the more widely used dosim- 
eter films under conditions simulating as closely as 
possible the temperatures and relative humidities 
that may be encountered during its use in personnel 
dosimetry. The films were used in their original 
packets, which had been sealed commercially into 
polyethylene For the irradiation, one set 
of packets was removed from the polvethylene bags, 
the other set remained protected. All film packets 
were exposed to Co® gamima radiation of low in- 
tensity for a period of about 1 month under 1 of 12 
different combinations of atmospheric temperature 
and relative humidity. After completion of the 
exposures, all film packets were enclosed in poly- 
ethvlene and placed in a refrigerator. About 5 
days later they were developed according to con- 
ventional procedure, along with control films exposed 
at room temperature. 


bags. 


iby the U.S 
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Sealed polyethylene bags afford con- 


Tomoda found that different film types behaved 
similarly under his experimental conditions [5]. 
The experimental conditions employed here differ 
from those of Ziegler and Chleck and also from those 
of Tomoda, who exposed at room temperature and 
heated the films after exposure. However, it is 
reasonable to assume that different types of film 
will also behave similarly under the conditions of 
the present experiment. Therefore, the results 
obtained here for one dosimeter film should be 
applicable to other film types as well. 


2. Experimental Technique 


A sealed Co® source with an activity of about 
50 we was placed in the center of a spiral array of 
lead-enclosed wooden supports, each holding a 
sealed glass jar containing a small amount of satu- 
rated salt solution that produced a known relative 
humidity within the jar. Four different relative 
humidities were obtained by using lithium chloride, 
magnesium chloride, sodium chloride, and potas- 
sium sulfate, listed in the order of increasing values 
of relative humidities. The exact humidity value 
for each of the three temperatures used was taken 
from the data of Wexler and Hasegawa [6]. At the 
temperatures employed for the experiments, the 
amount of solutes present in the gaseous phase was 
negligible. 

Each glass jar contained two film packets, one bare 
and one sealed in a polyethylene bag. The packets 
were supported by small clamps glued to the inside 
of the jar lids. The walls of the jars were thick 
enough to establish electronic equilibrium for Co® 
radiation at the inner wall surface. They absorbed 
about 2 percent of the incident Co® gamma rays. 

A total of 24 jars was used. The source-to-jar 
distances were chosen in such a way as to yield 5 
different film exposures at each of the 4 different 
relative humidities; also 1 jar at each relative 
humidity was shielded from the radiation. 

The entire setup was placed in a room whose 
temperature was held as closely as possible at a 
constant value for 1 month. About 5 days after 
completion of the exposure, the films were developed 
in Kodak Liquid X-ray Developer for 5 min at 20 °C, 
along with ‘control films”. The control films were 
exposed at room temperature to the radiation from 
a stronger Co® source about 4 days prior to develop- 
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ment;-all control packets were enclosed in poly- 
ethylene and were supported from the lids of dry 
jars during the exposures which lasted for 1 to 30 min. 
The same procedure was followed for 3 consecutive 
months, during which the exposure room 
maintained at 3 different temperatures. 


3. Results 


The table shows the net densities (i.e., the densities 
above base and fog) obtained at the 4 different 
exposure levels and for each of the 12 exposure con- 
ditions; the listed values were taken from curves of 


was 





the plotted experimental data. The uncertainty in 
the exposures is estimated to be less than 10 percent, 
while that in the density readings is less than 0.02 
density units. The tabulation enables the reader to 
plot, for purposes of comparison, one characteristic 
curve for each exposure condition. However, com- 
parisons should be made only between the densities 
obtained at one temperature, since only the films 
exposed at the same temperature were developed 
simultaneously. 

The table also lists the photographic sensitivities 
relative to the sensitivity of the controls for all 
exposure conditions.2, At 5.5 °C all sensitivities 
but those obtained at 98.3 percent relative humidity 
were within 10 percent of the values for the controls. 
However, at 98.3 percent relative humidity, the 
sensitivity was considerably lower. At this high 
humidity, the film exposed in a polyethylene bag to 
0.25 r reached the same density as the film exposed 
to 1.0 rin a bare packet. The 12-percent difference 
between the sensitivity of the controls and that 
obtained at 98.3 percent relative humidity with the 
film packets inside the polyethylene bags may be 
real, inasmuch as polyethylene is not entirely 
impermeable to water vapor. 





2 Photographic sensitivity is here defined as the reciprocal of the exposure in r 
required for a given density. The sensitivity values used for the tables are 
averages of the sensitivities determined at the density levels 0.2 and 0.5. 


At 19.8 °C, the influence of high humidity was 
even more pronounced. Here, 0.57 given to a film in 
a bare packet in an atmosphere of 97.2 percent rela- 
tive humidity produced no measurable density above 
background, while at lower relative humidities the 
bare films showed a tendency to fog. The sensitivity 
of the films exposed in polyethylene-enclosed packets 
was not significantly different from that of the 
control. 

At 32.5 °C, there was a significant increase in 
the sensitivity of all films, as much as 30 percent 
for films protected from atmospheric changes. 
This was to be expected on the basis of the kineties 
of the photographic process [7]. However, super- 
imposed on this effect was, in the case of the bare 
packets, the effect of the varying atmospheric con- 
ditions. As a result, the films exposed in bare 
packets at relative humidities up to 75.6 percent 
showed an additional sensitivity increase of up to 
15 percent. At a relative humidity of 96.5 percent, 
the films were tightly stuck to adjacent films con- 
tained in the same packets. During the process 
of separation, they were damaged considerably. 
The undamaged portions of the exposed films were 
fogged to a degree that made evaluation impossible. 
The density of the corresponding unexposed film, 
although higher than that of the film in the pro- 
tected packet, was measurable. The density of the 
heated, unexposed film from the packet protected 
by polyethylene was not appreciably higher than 
that of the unexposed control film, which had been 
kept at room temperature. 


4. Conclusions 


Regardless of temperature, bare film packets 
cannot be used for satisfactory personnel dosimetry 
in areas where prolonged exposure to relative 
humidities close to a hundred percent are expected. 
Such high relative humidities cause physical damage 
to the film packets; moreover, at low temperatures, 


Change in film response with temperature and relative humidity 


Net density for exposure of the jars to 


= = Relative sensitivity 
Relative 
Temperature humidity 0.10 r 0.25 r 0.50 r 1OOr 
Bare Plastic Bare Plastic Bare Plastic sare Plastic Bare Plastic 
C t 

543.5 13. 9+1 0. 02 0. 04 0. 08 0.09 0. 22 0. 22 0. 46 0. 43 1.00 1. 05 
34. 542 03 03 10 10 23 23 44 44 9S YS 
75. 1+0.6 03 03 09 Oy 21 21 45 45 06 06 
98. 34 7 a() 02 a() wo 2 02 19 09 40 very low RS 
03 09 22 4 46 ¢1.00 
19.8-4-0.5 12. 4+0.5 0. 03 0. 03 0.10 0. OR 0.19 0.16 0. 42 0. 37 25 0. OS 
33.5% .5 . 03 03 OS 07 16 15 37 35 1. 00 93 
75.54 .6 03 03 10 OS 21 18 44 39 1.16 1. 02 
97.24 .5 a) 03 e() 09 a() 17 (2 38 very low 1. 02 
©. 03 07 17 38 ©1.00 
32.541 11. 740.5 0. 06 0. 03 0. 12 0.10 0. 27 0. 23 0. 54 0. 49 1.45 1. 29 
32.64 .5 O4 04 12 12 21 21 44 44 1. 22 1.22 
75.64 .7 06 03 13 O8 26 20 49 1) 1. 38 1.07 
96. 5 5 21.70 06 21.68 12 21.75 22 21.85 44 ds 4 1. 22 
02 O8 Is 37 ©1.00 


* The films stuck to adjacent films during exposure and had to be forced apart with a knife 


> Numerical evaluation impossible. 
Controls, exposed at laboratory temperature in dry jars. 
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they produce excessive fading of the latent image 
and, at high temperatures, are the cause of emulsion 
fogging. 

Personnel badge calibrations are usually made 
at normal laboratory temperatures, regardless of 
the ambient temperatures at which the personnel 
badges are to be used. When the temperatures 
are higher during use, one expects the personnel 
exposures to be overestimated. This error is less 
when the film packets that are to be used in the 
personnel badges are sealed in polyethylene bags. 
Thus, although polyethylene not completely 
impermeable to water vapor, it does provide a 
considerable measure of protection to photographic 
film exposed under conditions of excessive atmos- 
pheric humidity. 


Is 
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Evaluation of the NBS Unit 
computable capacitor, R. D. 
65A (Phys. and Chem.) No. 3, 147 ( Vay-—J une 
An evaluation of the unit of resistance maintained at the 
National Bureau of Standards, based on the prototype 
standards of length and time, is described. The evaluation 
is based on a nominally one-picofarad capacitor whose value 
may be calculated from its mechanical dimensions to high 
accuracy. This capacitor is used to calibrate an 0.01- 
microfarad capacitor. A frequency-dependent bridge involv- 
ing this capacitor establishes the value of a 104-ohm resistor. 
Comparison of that resistor with the bank of onae-ohm resistors 
maintaining the NBS unit of resistance establishes that this 
unit is 


based on a 
Research NBS 


1961) 75 cents. 


of resistance 
Cutkosky, J. 


ev = 1.000002; ohms + 2.1 ppm. 


The indicated uncertainty is an estimated 50 percent error 
of the reported value based on the statistical uncertainty of 
the measurements and allowing for known sources of possible 
systematic errors other than in the speed of light, and assum- 
ing that the speed of light is c= 2.997925 « 10'cm/see. 


Digital methods for the extraction of phase and amplitude 
information from a modulated signal, R. 8S. Lawrence, 
J. L. Jespersen, and R. C. Lamb, J. Research NBS 65D 
(Radio Prop.) No. 4, 351 (July-August 1961) 765 cents. 

A description is given of three digital methods which have 
been used to recover amplitude and phase information from 
a modulated sinusoidal signal sampled at equal intervals of 
not more than one-sixth of a period. The first method, the 
“zero-crossing’’ method, is economical of computer time and, 
for modulation which is not too deep and does not contain 
frequencies near the carrier frequency, accurately 
the phase and amplitude modulation. The second method 
the ‘filter’? method, is more laborious but it gives better 
accuracy and will operate with deeper and more rapid modula- 
tion. The third method, a statistical approach, will work 
with severly overmodulated signals, but it yields only a 
statistical summary of the modulation. The methods were 
designed specifically for analysis of radio-star scintillation 
records but they may be applied to many other modulated 
signals, 


recovers 


Development of high-temperature strain gages, J. 
and D. G. Moore, NBS Mono. 26 (1961) 20 cents. 
{ summary is presented of a research program aimed at the 
improvement of high-temperature strain gages of the electrical 
resistance type. Potential ceramic and metal components 
were evaluated and a gage was devised that was based on 


W. Pitts 


these evaluations. This gage (NBS 5B) was flexible and 
easy to install; however, it lacked resistance stability at 
higher temperatures. In an attempt to minimize this 


deficiency, ceramic cements were developed that showed 
greater electrical resistivity than had been previously observed 
in the range 800 to 1,800 °F; also, a technique was devised 
for increasing the resistance to ground by applying a fired-on 
ceramic coating to the grid of a specially developed unbacked 
gage 

A study was made of the cause of the erratic response of 
cemented gages that had not been preheated prior to use. 
There were strong indications that the erratic response was 
caused mostly by the rapid decrease in resistance that 
accompanied structural changes in the cement. 





Precision measurement and calibration, 5S. F. Booth: 
Electricity and electronics, NBS Handb. 77, Vol. I (1961) 
$6.00; Heat and mechanics, Vol. JJ (1961) $6.75; Opties, 
metrology, and radiation, Vol. J] J (1961) $7.00. 

This Handbook is a compilation of the more important 


National Bureau of Standards publications over a period of 
years, dealing with precision measurement and the calibration 
of standards. The publications, originally issued as circulars, 
research papers, chapters of books, and articles in scientific 
periodicals, were selected for inclusion in the Handbook on 
the basis of their having, in the past, best served the needs 
of scientists in specialized fields of measurement. It is 
intended to serve as a quick reference source for workers in 
the field of standards, and also as a textbook and aid to 
scientists and engineers in standards laboratories. The list 
of specific titles included in each volume is given in the Table 
of Contents of that volume. Because the three volumes are 
sold separately, the subject index and authors’ index are 
complete and identical in each volume. 


Units of weight and measure (United States customary and 
metric) definitions and tables of equivalents, L. V. Judson, 
NBS Misc. Publ., 233 (1960) 40 cents. 

The units of length, area, volume, capacity, and mass in the 
United States are defined. Tables of interrelation and tables 
of equivalents for these units in the metric system and in the 
U.S. customary system are given. All values in the text 
and in the tables have been revised to be in accord with the 
Federal Register announcement of July 1, 1959, an announce- 
ment entitled ‘Refinement of Values for the Yard and 
Pound” in accordance with an agreement among the directors 
of National Standards Laboratories of English-speaking 
nations to obtain uniformity in precise measurements involv- 
ing the yard and the pound. The long tables were all recom- 
puted on an automatic computer and printer. Material on 
the status of the metric system in the United States was 
deleted as not being germane to the main purpose of the 
publication. The tables showing interrelation (1 to 999 
units) betweer bushels and hectoliters were deleted. A 
more complete table (0.01 to 19.99 units) showing equivalents 
of inches in millimeters has been added. This is a revision 
of Miscellaneous Publications 214. 


Report of the 45th national conference on weights and 
measures 1960, NBS Misc. Publ. 235 (1960) 765 cents. 

A report of the proceedings of the forty-fifth National 
Conference on Weights and Measures, held in Washington, 
D.C., June 6, 7, 8, 9, and 10, 1960, and attended by state, 
county, and city weights and measures officials. 159 pages. 


The sample, the procedure, and the laboratory, W. J. Youden, 
Anal. Chem. 32, No. 13, 23A-—37A (Dec. 1960). 

Statistics is accorded a place in the activities of analytical 
chemists because analysts are interested in tracking down the 
sources of error in their determinations. Uniformity of the 
material samples, the inherent capacity of a procedure to give 
precise results free from systematic errors and the careful 
adherence of the analyst to the prescribed procedure are all 
involved in the final result. When different laboratories 
analyse the same materials the results usually show much 
greater discrepancies than that found for determinations 
made by one laboratory. 

This paper describes (1) a logical breakdown of the analytical 
error into meaningful components: and (2) a graphical method 
for presenting the results that clearly differentiates between 
random and systematic errors. 
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Evaluation of resistance strain gages at elevated tempera- 
tures, R. L. Bloss, Am. Soc. Testing Materials Proc. 1, No. 1, 
9-15 (Jan. 1961). 

To satisfy the need for strain measurements at elevated 
temperatures a number of types of “‘high temperature”’ strain 
gages are available. A comprehensive evaluation of these 
gages is needed because of the lack of available information 
on their performance and because of difficulty which has been 
encountered in their use. A facility for conducting such an 
evaluation is described. Typical results obtained during 
evaluation tests are given to show the capability of the equip- 
ment, to point out some of the problems encountered in the 
use of these gages, and to illustrate the need for comprehensive 
gage evaluation prior to use. The results shown are only 
indicative of the performance of the particular gage tested. 
Significantly different results might be obtained from other 
gage types. 


Techniques in calorimetry. I. A noble-metal thermocouple for 
differential use, E. D. West, Rev. Sci. Inst. 31, No. 8, 896- 
897 (Aug. 1960). 

The application of a noble-metal thermocouple to problems 
of the observation and control of temperature differences in 
calorimetry is described. Fora couple of gold-40%-palladium 
and platinum-10°%-rhodium dE/dT is 37 to 57 microvolts/°C 
from 25° to 1100°C. The uncertainty in the temperature 
difference measured by a 20-junction thermopile in a nearly 
isothermal region is less than 0.4 millidegrees. 


Evaluation of micrometer and microscopical methods for 
measuring thickness of floor coverings, 1. Horowitz, J. 
Mandel, R. J. Capott, and T. H. Boone, ASTM Bull. 1, No. 2, 
99-102 (Feb. 1961 

The dial micrometer and microscopical (optical) methods for 
measuring the thickness of floor coverings were evaluated in 
an interlaboratory study of four samples having wearing 
layer thicknesses in the range 5.5 to 22 mils. The data are 
analyzed statistically in terms of three parameters: (1) The 
overall average, wu; (2) the slope of the straight line obtained 
by plotting the average results of a given participant versus 
the corresponding grand average of all the participants, 8; 
and (3) the scatter of the individual results about the straight 
line, V(»).. With experienced personnel, there is no significant 
difference in the results obtained by the two methods; the 
grand averages agree to within 0.04 mils. Using four replicate 
measurements the standard deviations for the micrometer 
and optical methods are 1.15 and 1.19 mils, respectively. 
The measurement of thickness does not appear to be affected 
by different diameter pressure feet between the limits of 0.10 
and 0.25 inches when the load on the pressure foot is main- 
tained between 18.6 and 22.4 ]b/in?. 


Standards and testing, the key to quality, B. £. Foster, Eng: 
Bull., Purdue Univ. XLIV, No. 6, 54-62 (Nov. 1960). 
Measurement standards, dimensional standards, quality 
standards, and construction standards are briefly discussed. 
The advantages of quality standards in furnishing guidance to 
the selection of materials, in promoting uniformity of materials 
purchased, and from the contractural standpoint are brought 
out. The importance of testing is emphasized. 


An analysis of the accumulated error in a hierarchy of cali- 
brations, E. L. Crow, 7RE Trans. Instrumentation I-9, No. 2, 
105-114 (Sept. 1960). 

Calibrations of many types are performed in a hierarchy of 
-alibration laboratories fanning out from a national standard. 
Often the statement is made that the accuracy of each echelon 
of the hierarchy should be 10 times the accuracy of the im- 
mediately following echelon. The validity of such statements 
is examined by deriving formulas for the total error accumu- 
lated over the entire sequence when systematic and random 
errors may occur in each echelon, and by determining how a 
given total error may be achieved at minimum total cost under 
reasonable assumptions for the form of the cost-error functions. 


A recommended standard resistor-noise test system, G. T. 
Conrad, Jr., M. Newman, and A. P. Stansbury, 7RE Trans. 
Component Parts, CP-7, No. 3, 71-88 (Sept. 1960). 

This paper describes a recommended standard resistor noise 
test system, with the theroy of operation given in considerable 








detail. The system provides for the measurement of a defin- 
itive resistor noise quality Index. This Index, the ‘“‘micro- 
volts-per-volt’’ Index, is a modified form of the “conversion 
gain’’ index formerly used at NBS. The Index is a decibel 
expression for the number of microvolts of noise per volt of 
applied de voltage transmitted in a frequency decade. A 
practical example is given to illustrate how the Index may be 
used to calculate the magnitude of current noise generated 
by fixed resistors in a given application. To complete the 
example, the combined effects of thermal noise, tube noise, 
and current noise are then considered. 

A description is also given of an existing test set which is based 
upon specifications associated with the recommended measur- 
ing system. The full specifications for this set are given and 
distinctive design features are discussed. The simple opera- 
tion of the test set and the computation of the Index are 
carefully explained. 


A transfer instrument for the intercomparison of microwave 
power meters, C. F. Engen, JRE Trans. Instrumentation I-9, 
No. 2, 202-208 (Sept. 1960). 

In the intercomparison of two microwave power standards, 
or in the calibration of a microwave power meter by means of 
a second or ‘‘standard’’ power meter, the measurement of 
microwave impedance has played a major role. 

Through an extension of the reflectometer concept it is now 
possible to devise a four arm junction which, when properly 
adjusted, makes possible the intercomparison or calibration 
of such power meters with little or no regard for their im- 
pedance characteristics. In addition, the method is substan- 
tially independent of the impedance discontinuity which may 


be present at the input flange or connector. This latter 
result is of particular value in coaxial systems. 
Propagation of error in a chain of standards, A. G. McNish 


and J. M. Cameron, 7RE Trans 
101-104 (Sept. 1960). 

Methods for dealing with errors developed in individual lab- 
oratories in a chain of standardizing laboratories must be 
considered with respect to the end result. It is shown that if 
a large number of laboratories is involved in the chain the 
most satisfactory results are obtained by regarding the total 
error of the chain as the root of the sum of the squares of the 
errors of the individual laboratories. If only a few laboratories 
are involved it is safer, and satisfactory too, to add the indi- 
vidual errors. 


Instrumentation I-9, No. 2, 


Other NBS Publications 


Journal of Research 65A 

April 1961) 75 cents. 

Mass spectra of some deuteroethanes. 
L. Mohler. 

Heats of hydrolysis and formation of potassium borohydride. 
W. H. Johnson, R. H. Schumm, I. H. Wilson, and E. J. 
Prosen. 

Heat of combustion of borazine B;N;Hg. 
H. Johnson, and E. J. Prosen. 

Thermodynamic properties of thorium dioxide from 298 to 
1,200 °K. A.C. Victor and T. B. Douglas. 

Calculated energy dissipation distribution in air by fast 
electrons from a gun source. J. E. Crew. 

Vitrons as flow units in alkali silicate binary glasses. L. 
Tilton. 

Tetragermanates of strontium, lead, and barium of formula 
type AB,Oy. C. R. Robbins and E. M. Levin. 


Phys. and Chem.) No. 2 (March- 


E. I. Quinn and F, 


M. V. Kilday, W. 


W. 


Journal of Research 65A (Phys. and Chem.) No. 3 (May- 

June 1961) 75 cents. 

International practical temperature scale of 1948, 
vision of 1960. H. F. Stimson. 

Evaluation of the NBS unit of resistance based on a com- 
putable capacitor. R.D.Cutkosky. (See above abstract.) 

Wavelengths and intensities in the first spectrum of bromine, 

2000 to 13000 A. J. L. Tech and C. H. Corliss. 

Torsional resonance vibrations of uniform bars of square cross 
section. W. E. Tefft and 8. Spinner. 


Text re- 
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Infrared studies of aragonite, 
structures in the borates, carbonates, 
Weir and E. R. Lippincott. 

Dielectric properties of polyamides: polyhexamethylene adip- 
amide and polyhexamethylene sebacamide. A. J. Curtis. 

Heat of formation of calcium aluminate monocarbonate at 
25°C. H. A. Berman and E. 8. Newman, 

Thermodynamic constants for association of isomeric chloro- 
benzoic and toluic acids with 1,3-diphenylguanidine in 
benzene. M. M. Davis and H. B. Hetzer. 

Heats of combustion and formation of trimethylborane, 
ethylborane, and tri-n-butylborane. W. H. Johnson, 
V. Kilday, and E. J. Prosen. 

Pyrolysis of linear copolymers of ethylene and propylene. 
S. Straus and L. A. Wall. 

Pyrolysis of fluorocarbon 
Straus. 

Preparation of fluoro- and bromofluoroaryl compounds by 
copyrolysis of bromofluoroalkanes. L. A. Wall, o: & 
Fearn, W. J. Pummer, and R. E. Lowry, 

Thermal stability of polydivinylbenzene and of copolymers 
of styrene with divinylbenzene and with tribinylbenzene. 
S. Straus and 8. L. Madorsky. 

Conformations of the pyranoid sugars. 
tion spectra 7 some fully 
Tipson and H. 8. Isbell. 

A standard for the measurement of the pH of blood and other 
physiological media. V. E. Bower, M. Paabo, and R. G., 
Bates, 


calcite, and vaterite type 


and nitrates. C. E. 


tri- 
M. 


polymers. L. A. Wall and 8. 


LY, 


acetylated 


Infrared absorp- 
pyranoses, R. S. 


Journal of Research 65B (Math. and Math. Phys.) No. 2 

(April-June 1961) 70 cents. 

Optimal approximation for functions prescribed at equally 
spaced points. H. F. Weinberger. 


Truncations in the method of intermediate problems for 


lower bounds to eigenvalues. N. W. Bazley and D. W. 
Fox. 

Comparison — sms for symmetric functions of character- 
istic roots By tga 

Some prope sae of the empirical distribution function of a 


random process. M. M. Siddiqui. 

Another extension of Heinz’s inequality. M. Marcus. 

Mean motions in conditionally periodic separable systems. 
J. P. Vinti. 

Some boundary value problems involving plasma media. J. 
R. Wait. 

A new decomposition formula in the 
J. H. Bramble and L. E. Payne. 

Pointwise bounds in the Cauchy problem of elastic plates. 
L. E. Payne. 


theory of elasticity. 


Journal of Research 65D (Radio Prop.) No. 4 
1961) 75 cents. 

Almost fifty years of URSI. J. H. Dellinger. 
Power density requirements for airglow excitation by gyro- 
waves. V. A. Bailey. 
On the validity of some 


(July-August 


approximations to the 
Hartree formula. K. Davies and G. A. M. King. 

Amplitude and angular scintillations of the radio source 
Cygnus-A observed at Boulder, Colorado. R.S. Lawrence, 
J. L. Jespersen, and R. C. Lamb. 

Digital methods for the extraction of phase and amplitude 
information from a modulated signal. R. 8S. Lawrence, 
J. L. Jespersen, and R. C. Lamb. (See above abstract.) 

Comparison between mode theory and ray theory of VLF 
propagation. H. Volland. 

Antenna coupling error in direction finders. W. Harrison, Jr. 

The electrically short antenna as a probe for measuring free 
electron densities and collision frequencies in an ionized 
region. R. W. P. King, C. W. Harrison, Jr., and D. H. 
Denton, Jr. 

Effect of multiple atmospheric inversions on tropospheric 
radio propagation. F. H. Northover. 

A few observations of the perturbations in the phase of the 
low-frequency ground wave. J. M. Ross and J. EK. Kirch. 
Smooth earth diffraction calculations for horizontal polariza- 

tion. L. E. Vogler. 

On the theory of mixed-path ground-wave propagation on 
aspherical earth. J. R. Wait. 


Appleton- 








Atomic energy levels in crystals, J. I 
19 (1961) 60 cents. 

Climatic charts and data of the radio refractive index for the 
United States and the world, B. R. Bean, J. D. Horn, and 
A. M. Ozanich, Jr., NBS Mono. 22 (1960) $2.00. 

A spectrophotometric atlas of the spectrum of CH from 3000A 
to 5000A, A. M. Bass and H. P. Broida, NBS Mono. 24 
(1961) 20 cents. 

Measurement of absorbed dose of neutrons, and of mixtures 
of neutrons and gamma rays, NBS Handb. H75 (1961) 
35 cents. 

Medical X-ray protection up 
Handb. H76 (1961) 25 cents. 

Report of the International Commission on Radiological 
Units and Measurements (ICRU) 1959, NBS Handb, H78 
(1961) 65 cents. 

Quarterly radio noise data—June, July, 
Crichlow, R. T. Disney, and M. A. 
(PB151377-7) (1960) $1.75. 

Quarterly radio noise data, September, 


1. Prather, NBS Mono. 


to three million volts, NBS 


August 1960, W. Q. 
Jenkins, NBS TN18-7 


October, November 


1960, W. Q. Crichlow, R. T. Disney, and M. A. Jenkins, 
NBS TN18-8 (PB151377-8) (1961) $1.75. 
VHF radio propagation data for Cedar Rapids-Sterling, 


Anchorage-Barrow, and rie € Churchill test paths April 

1951 through June 1958, R. Sugar and K. W. Sullivan, 
NBS TN79 (PB161580) C960) $4.00. 

Bibliography of tropospheric radio wave 
Abbott, NBS TN80 (PB161581) 

A survey of spread-f, F. M. Glover, 
(1960) $1.75. 

The NBS meteor-burst propagation project 
port, C. E. Hornback, L. D. Brayfogle, and G. R. Sugar, 
NBS TN86 (PB161587) (1960) $1.25. 

A theoretical aah of spor: adic-E structure in the light of radio 
measurements, K. Tao, NBS TN87 (PB161588) (1961) 
$1.25. 


scattering, R. L. 
(1960) $2.25. 


NBS TN82 (PB161583) 


a progress re- 


Prolonged space-wave fadeouts in tropospheric propagation, 
\. P. Barsis and M. E. Johnson, NBS TN88 (PB161589) 
(1961) $2.00. 

Some ceramic dielectrics with a very low temperature coeffi- 
cient of capacitance, 8S. Marzullo and E. N. Bunting, J. Am. 
Ceram. Soc. 43, No. 11, 609, (Nov. 1960). 

Vibration-rotation bands of H.O, E. D. Tidwell and E. K. 


Plyler, J. Opt. Soe. Am. 50, No. 12, 717-720 (Dee. 1960). 
Nuclear magnetic resonance in tantalum metal, J. I. Budniek 
and L. N. Bennett, J. Phys. Chem. Solids 16, No. 1/2, 
37-38 (Jan. 12, 1960). 
The relation of hmarF2 to M(3000)F2 and H,F2, J. W. 
Wright and R. E. McDuffie, J. Radio Research Lab. 7, 
No. 32, 409-420 (July 1960). 


S-quinolinol precipitation of the elements, J. I. 
Chem. Anal. 49, No. 4, 126 (Dee. 1960). 

Dipole moments of hydrocarbons, D. R. Lide, Jr., J. 
Phys. 33, No. 6, 1879 (Dee. 1960). 

Collection of ions produced by alpha particles in air, Z. 
Bay and N. H. Seliger, Phys. Rev. 120, No. 1, 141-143 
(Oct. 1, 1960). 

A test of a procedure for easy estimation of representative 

monthly electron density profiles for the ionosphere, J. 
W. Wright, J. Geophys. Research 65, No. 10, 3215-3217 
(Oct. 1960). 

A method for efficiently providing low temperature liquids on 
a large scale to an accelerated experimental program, A. 
K. Stober, Proc. 10th Intern. Congress of Refrigeration, 
Copenhagen, Denmark 1, No. 1, 17-18 (Pergamon Press 
Ltd., London, England, 1959). 

Seasonal variations in the twilight enhancement of [OI] 
L. R. Megill, P. M. Jamnick and J. E. Cruz, J. 
and Terrest. Phys. 18, 309-314 (Aug. 1960). 

Synthetic mica, A. V. Valkenburg, Book, Encyclopedia of 
Chemical Technology, p. 480-487 (Interscience Encyclo- 
pedia, Inc., New York, N.Y., 1960). 

Nonresonant microwave absorption and relaxation frequency 
at elevated pressures, > A. Maryott and G. Birnbaum, J. 
Phys. Chem. 64, No. 11, 1778-1780 (Nov. 1960). 

Precipitation of the > oeca on addition of aqueous ammonia 


Hoffman, 


Chem. 


cee 


5577, 
Atmospheric 


to their clear solution, J. I. Hoffman, Chem.-Anal. 49, 
No. 3, 94 (Sept. 1960). 


Low scatter high current gas target for D-D neutrons, A. 


209 








C. B. Richardson, Rev. Sci. Inst. 31, No. 11, 1202-1203 
(Noy. 1960). 

The heat of combustion of dicyanoacetylene, G. T. Armstrong 
and 8S. Marantz, J. Phys. Chem. 64, 1776-1777 (1960). 

Errors in dielectric measurements due to a sample insertion 
hole in a cavity, A. J. Estine, and H. E. Bussey, IRE 
Trans. Microwave Theory and Tech. MTT-8, No.6, 
650-653 (Nov. 1960). 

Ensemble method in the theory of irreversibility, R. Zwanzig, 
J. Chem. Phys. 33, No. 5, 1338-1341 (Nov. 1960). 

Determination of pentosans. Interlaboratory comparison of 
the aniline acetate, orcinol, and bromination methods, W. 

<. Wilson and J. Mandel, Tappi 43, No. 12, 998-1004 
(Dec. 1960). 

A research for geomagnetic singular days, C. 
Pohrte, and H. MacDonald, J. Geophys. 
No. 9, 3013-3015 (Sept. 1960). 

Letter to editor of Revue Des Matériaux (first), E. S. Newman 
and H. A. Berman, Revue Des Matériaux C, No. 540, 231 
(Dec. 1, 1960). 

Geomagnetic storms and the space around the earth, S. 
Chapman, Nature 187, No. 4740, 824-827 (Sept. 3, 1960). 

On the theory of the critical point of a simple fluid, M. 8. 
Green, J. Chem. Phys. 33, No. 5, 1403-1409 (Nov. 1960). 

Some implications of slant-E,, E. K. Smith and R. W. Knecht, 
Polar Atmosphere Symp., Part 2, Ionospheric Section, 
Oslo, Norway, July 2—8, 1956, p. 195-204 (Pergamon Press, 
Inc., New York, N.Y., 1957). 

Mismatch errors in microwave phase shift measurements, G. 
E. Schafer, IRE Trans. Microwave Theory and Tech. 
MTT-8, No. 617-622 (Nov. 1960). 

Residual are spectra of seventy elements diluted in copper, 
C. H. Corliss, W. F. Meggers and B. F. Scribner, Book, 
Colloquium Spectroscopicum Internationale VIII, 119-121 
(1959). 

The rotational constants of hydrogen chloride, E. 
and E. D. Tidwell, Z. Elektrochem. 64, No. 
(Feb. 1960). 

The minima of cyclic sums, K. Goldberg, J. 
Soc. 35, 262-264 (1960). 

Normal approximation to the chi-square and non-central F 
probability functions, N. C. Severo and M. Zelen, Bio- 
metrika 47, No. 3 & 4, 411-416 (1960). 

Microtechnique for the infrared study of solids, diamonds, and 
sapphires as cell materials, E. R. Lippincott, F. E. Welsh, 
and C. E. Weir, Anal. Chem. 33, 137-143 (Jan. 1961). 

Definitions relating to metals and metalworking, M. R. 
Meyerson and 8. J. Rosenberg, Am. Soc. Metals Handb. 1, 
1—41 (1961 

Stress-rupture tests at 1350° F on type 304 stainless steel, 
W. D. Jenkins, W. A. Willard and W. J. Youden, ASTM 
Bull. 1, No. 2, 104-108 (Feb. 1961). 

Paratellurite, a new mineral from Mexico, G. Switzer and 
H. E. Swanson, Am. Mineralogist 45, 1272 (Nov.—Dec. 
1960). 

The nature, cause, and effect of porosity of electrodeposits 


Warwick, T. 
Research 65, 


K. Plyler 
5, 717-720 


London Math. 


Vi Note on a magnetic method of detecting corrosion 
currents, F. Ogburn and W. H. Roberts, Plating 48, No. 
2, 168-169 (Feb. 1961 


The height of maximum luminosity in an auroral are, F. E. 
Roach, J. G. Moore, E. C. Bruner, Jr., H. Cronin, and S. M. 
Silverman, J. Geophys. Research 64, No. 11, 3575-3580 
(Nov. 1960). 

Physical metallurgy and mechanical properties of materials: 
Ductility and the strength of metallic structures, J. M. 
Frankland, J. Eng. Mech. Div. Proc. Am. Soc. Civil Eng. 
86, No. EM 6, 45-52 (Dec. 1960). 

Cryogenic impurity adsorption from hydrogen, M. J. 
Chem. Engr. Progress 55, No. 10, 68-71 (Oct. 1960). 

Nomenclature for standards of radioactivity, A. G. McNish, 
Intern. J. Appl. Radiation and Isotopes 8, 145-146 (Jan. 
1960). 

Lithium, M. R. Meyerson, Am. Soc. for Metals Handbook 1, 
1213 (1961). 

Tensile cryostat for the temperature range 4° to 300 
R. M. McClintock and K. A. Warren, ASTM 
No. 2, 95-98 (Feb. 1961 


Hiza, 


Kelvin, 


Bull. 1, 


A radio-frequency permittimeter, R. C. Powell and A. L. 
Rasmussen, IRE Trans. Instrumentation I-9, No. 2, 


179-184 (Sept. 1960). 








Subgroups of the modular group and sums of squares, H. 
Newman, Am. J. Math. 82, No. 4, 761-778 (Oct. 1960). 
Phase equilibria in the system cadmium oxide-niobium oxide, 

R. 8. Roth, J. Am. Ceram. Soc. 44, No. 1, 49-50 (Jan. 1961). 

Redetermination of the chromium and nickel solvuses in the 
chromium-nickel system, C. J. Bechtoldt and H. C. Vacher, 
Trans. Metallurgical Soc., Am. Inst. Mining Engrs. 221, 
14-18 (Feb. 1961 

A modulated subcarrier technique of measuring microwave 
phase shifts, G. E. Schafer, IRE Trans. Instrumentation 
I-9, No. 2, 217-219 (Sept. 1960). 

Hydrogen-bonding in calcium-deficient hydroxyapatites, 
A. 5. Posner, J. M. Stutman, and E. R. Lippincott, 
Nature 188, No. 4749, 486-487 (Nov. 1960). 

Hydrogen sulfide precipitation of the elements from 0.2-0.5 
normal hydrochloric acid, J. I. Hoffman, Chemist Analyst 
50, No. 1, 30 (Mar. 1961). 

Dielectric constant and dielectric loss of TiO. (Rutile) at 


low frequencies, R. A. Parker and J. H. Wasilik, Phys. 
Rev. 120, No. 5, 1631-1637 (Dec. 1960). 
Generating functions for formal power series in non-com- 


muting variables, K. Goldberg, Proc. Am. Math. Soe. 11, 
No. 6, 988-991 (Dec. 1960). 

Tests for regression coefficients when errors are correlated, 
M. M. Siddiqui, Ann. Math. Stat. 31, No. 4, 929-938 
(Dec. 1960). 

Application of the method of polarized orbitals to the scatter- 
ing of electrons from hydrogen, A. Temkin and J. C. 
Lamkin, Phys. Rev. 121, No. 3, 788-794 (Feb. 1961). 

On the absorption spectrum of CF» and its vibrational analysis, 
D. E. Mann and B. A. Thrush, J. Chem. Phys. 33, No. 6 
1732-1734 (Dec. 1960). 

Standards and measurements of microwave surface imped- 
ance, skin depth, conductivity and Q, H. E. Bussey, IRE 
Trans. Instrumentation I-9, No. 2, 171-175 (Sept. 1960). 

Plating gun bores, V. A. Lamb and J. P. Young, Ordnance 
XLV, No. 245, 725-727 (Mar.—Apr. 1961). 

Integrated starlight over the sky, F. E. Roach and L. 
Megill, Astrophys. J. 133, No. 1, 228-242 (Jan. 1961). 

Low temperature static seals using elastomers and plastics, 


> 


R. 


D. H. Weitzel, R. F. Robbins, G. R. Bopp, and W. R. 
sjorklund, Rey. Sci. Inst. 31, No. 12, 1350-1351 (Dec. 
1960). 

Electrodeposition of metals from nonaqueous media, A. 


Brenner, Encyclopedia of Chemical Technology, 2nd Suppl., 
315-324 (Interscience Encyclopedia 1960). 

FM and SSB radiotelephone tests on a VHF ionospheric 
scatter link during multipath conditions, J. W. Koch, 
W. B. Harding, and R. J. Jansen, IRE Trans. Commun. 
Systems CS-8, No. 3, 183-186 (Sept. 1960). 

Absolute measurement of temperatures of microwave noise 


A. J. Estin, C. L. Trembath, J. S. Wells, and 


sources, 


W. C. Daywitt, IRE Trans. Instrumentation I-9, No. 2, 
209-213 (Sept. 1960). 

Electron scattering in high magnetic field, A. H. Kahn, 
Phys. Rev. 119, No. 4, 1189-1192 (Aug. 1960). 


Measurement of reflections and losses of waveguide joints 
and connectors using microwave reflectometer techniques, 
R. W. Beatty, G. F. Engen, and W. J. Anson, IRE Trans. 
Instrumentation I-9, No. 2, 219-226 (Sept. 1960). 

Microwave spectrum of cis-difluoroethylene, V. W. 
J. Chem. Phys. 34, No. 1 (Jan. 1961 

Surface roughness of gold castings, J. J. Barone, R. L. Huff, 
and G. Dickson, Dental Progr. 1, No. 2, 78-84 (Jan. 1961). 

The sensitivity of photographic film to 3-mev neutrons and 

to thermal neutrons, M. Ehrlich, Health Phys. 4, 113 

128 (1960). 

prototype 

Carpenter, E. 

Stone, IRE 

1960). 

Low even configurations in the first spectrum of thorium 
(Th 1), R. E. Trees, Physica 26, 353-360 (1960). 

Comparison measurements with intensity standards for high 
energy bremsstrahlung, J. 8. Pruitt and W. Pohlit, Z. 
Naturforschung 15b, No. 9, 617-619 (1960). 

High-dispersion spectra of Jupiter, C. C. Kiess, C. H. Corliss 
and H. K. Kiess, Astrophys. J. 182, No. 1, 221-231 (July 
1960). 


Laurie, 


A rubidium vapor frequency standard, R. J. 
C. Beaty, P. L. Bender, 8. Saito, and R. O. 


Trans. Instrumentation I-9, No. 2 (Sept. 


210 


Propagation of electromagnetic waves along a thin plasma 
sheet, J. R. Wait, Can. J. Phys. 38, 1586-1594 (1960). 

Preparation and properties of aromatic fluorocarbons, W. J. 
Pummer and L. A. Wall, J. Chem. Engrs. Data 6, No. 1, 
76-78 (Jan. 1961). 

Stability of thermoset plastics at high temperatures, 8. L. 
Madorsky and 8. Straus, Modern Plastics 38, 134-140 
(Feb. 1961). 

Immiscibility and the system lanthanum oxide boric-oxide, 
E. M. Levin, C. R. Robbins, and J. L. Waring, J. Am. 
Ceram. Soc. 44, No. 2, 87-91 (Feb. 1961). 

On the properties of the vapor pressure curve, E. H. Brown, 
Cryogenics 1, No. 1, 37-40 (Sept. 1960). 





O 


Absorption and scattering of photons by holmiun and erbium, 
E. G. Fuller and E. Hayward, Proc. Intern. Conf. Nuclear 
Structure, 761—766 (Kingston, Canada, 1960). 

Spectroscopy of fluorine flames. I. Hydrogen-fluorine flame 
and the vibration-rotation emission spectrum of HF, 
D. E. Mann, B. A. Thrush, D. R. Lide, J. J. Ball, and 
N. Acquista, J. Chem. Phys. 34, No. 2, 420-431 (Feb. 1961). 


* Publications for which a price is indicated (except for NBS Technical Notes) are 
available only from the Superintendent of Documents, U.S. Government Printing 
Office, Washington 25, D.C. (foreign postage, one-fourth additional). Technical 
Notes are available only from the Office of Technical Services, U.S. Department of 
Commerce, Washington 25, D.C. (Order by PB number). Reprints from outside 
journals and the NBS Journal of Research may often be obtained directly from 
the authors. 


211 





oa a Be PENALTY FOR PRIVATE USE TE 
UNITED STATES a coe a PAYMENT OF POSTAGE, { 
GOVERNMENT PRINTING OFFICE — vat . 
S°IBION OF PUBLIC OCU MENTS 
WASHINGTON 115, D.C 


OFFICIAL BUSINESS 
Oo eae nanan mapanaaa 
JOURMAL CW RASEARCH 6b dhe Souci! Parecn of Sbade Vol, 65C, No. 3 
C. Engineering aad Instrumentation 





Contents 


Precliction of symptoms of cavitation. Robert B, Jacobs. 


Heating atid cocling of air flowing through an underground tunnel. 
B. A, Peavy. 


Stress-corrosion cracking of the AZ31B magnesium alloy. Hugh L. 
Logan. . 


Coatings formed on steel by cathodic protection and their evaluation 
by polarization measurements. W. J. Schwerdtfeger and Raul J. 
Manuele. 


Calibration of inductance standards in the Maxwell-Wien bridge 
circuit. Thomas L. Zapf. 


Calibration of loop antennas at VLF. A. G. Jean, H. E. Taggart, and 
J. R. Wait. 


Location of the plane of best average définition with low contrast resolu- 
tion patterns. Francis E. Washer and William P. Tayman. 


Influence of temperature and relative humidity on the photographic 
response to Co® gamma radiation. Margarete < rlich. 


Publications of the National Bureau of Standards. 








For salle by the Superintaucen ot Documents, U.S. Goverument Printing Otfice, Washington 25, D.C. Price 75 cents { 
Subscription priest $2.25 a year; 50 cents additional for foreign mailing. 


UNITED STATES GOVERNMENT ——_ OFFICE, WASHINGTON 1961 














4 : : ' : : 5 hone 








